收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

李海森等:多波束合成孔径声呐技术研究进展

[复制链接]
4 v( S2 ], q; {. e

来源:《测绘学报》

% j) L8 I4 B; S7 P# L

转自: 溪流之海洋人生

8 |) z3 x0 H3 ]( }

作者:李海森、魏波、杜伟东

3 i! ~ a, Z3 c. C2 G

近年来,随着现代水声信号处理技术和水声换能器技术的大幅度进步,水下目标精细探测和成像声呐技术已然成为了国内外研究的热点,在民用和军用领域都有着其他声呐不可替代的作用。在民用方面,成像声呐技术可用于海洋资源开发、海底地质勘探、海底地形地貌测绘、水下物体探测等海洋工程领域;在军事上,高隐蔽性水下军事小目标(如军用无人潜器、鱼雷、水雷、蛙人等)的探测与识别、港口锚地和舰艇的安全防范、地形匹配导航等领域上也迫切要求应用高分辨的水下目标精细探测和成像声呐技术。目前国内外已有多种先进的成像声呐技术,主流的主要包括干涉侧扫声呐技术、多波束测深声呐技术及合成孔径声呐技术等。

5 V/ x T5 [- p% h

干涉侧扫声呐一般需搭载水下拖体进行工作,其设备安装简单、目标横向分辨率较高,可以借助阴影对目标进行识别判断。但是,由于其探测机理制约不容易获得精确海底深度,并且测量垂底区域存在缝隙,需要单独的声呐设备或者方法进行补隙。多波束测深声呐能较精确地测量出海底深度并获得水体成像(watercolumn),能得到直观的、精确定位的全覆盖三维海底地形图,然而多波束测深声呐波束脚印随着深度增加而扩大,对远距离情况下的目标探测分辨率较低,对小目标的探测更为困难。合成孔径声呐(syntheticaperturesonar,SAS)使用小孔径的声呐换能器阵,通过运动形成虚拟大孔径的方法,来获取更高的航迹向分辨率。相比于实孔径声呐,SAS最突出的优势是航迹向分辨率与作用距离、信号的频率无关。然而,现阶段对于合成孔径声呐技术的研究主要集中在侧扫式合成孔径上,因此同样存在测深精度不佳和垂底探测缝隙等局限性。

$ P1 `( f8 `9 H$ u

/ p) }7 M3 y9 M: c) G

综上,迫切需要一种水下目标精细探测和成像声呐技术以满足对水下小目标探测能力的需求。为此本文将多波束测深声呐技术与合成孔径声呐成像技术相结合,在新的发射和接收基阵结构基础上,提出一种多波束合成孔径声呐(multibeamsyntheticaperturesonar,MbSAS)机理,理论分析和试验皆证明其可以获得与目标作用距离及发射信号频率无关的航迹向高分辨力,且可以精确测深和垂底区域没有缝隙。

1 多波束合成孔径技术研究现状# n5 G% S9 @8 s

多波束合成孔径技术是一种将多波束测深技术和合成孔径技术相结合的新型水下目标成像技术,通过载体运动在航迹向上虚拟合成较大的基阵孔径,既可以在航迹向上获取较高的分辨率,用于对地形地貌的全覆盖测量,还可以在距离向上通过波束形成确定目标所处方位,最终可以精确地测量出目标的深度信息,对目标进行三维成像。多波束合成孔径技术的发展,紧随着多波束测深技术和合成孔径技术的发展趋势,结合二者技术优势,实现水下目标的精细探测。

( \5 E: S2 G' l$ z5 E7 X' A

⒈ 多波束测深声呐技术现状与发展趋势

% u9 F; {: L3 a! j

多波束测深技术是随着现代水声、电子、计算机、信号处理技术的进步而发展起来的,至今多波束测深技术的研究已经经历了半个多世纪的发展,逐渐地形成了各种功能的实用化商业声呐产品。声呐系统供应商根据不同测量水深范围发展系列化的测深仪器,分为浅水、中水、深水多波束3类;按照搭载常规测量船只、水面无人船、水下AUV等不同设备载体研发抗压性、密闭性不同的换能器基阵,分为船载式、无人式和潜用式;根据不同的客户需求研发便携式悬挂基阵、V型组合基阵、内嵌式壳体基阵等不同的适装类型,极大地拓展了多波束测深系统的应用领域;从探测对象不同可以分为水面探测、水体探测和水底探测型;按照采用的测深信号处理算法不同又分为幅度检测法和相位检测法。概括起来,现阶段多波束测深技术的主要发展趋势是朝向超宽覆盖、小水深测量、运动姿态稳定、精细化测量等方向发展。

$ x! H8 n! l6 e; ^0 \

⑴超宽覆盖

( e9 D: a1 w0 K2 ^

限制多波束测深系统覆盖宽度的主要问题在于小掠射角情况下,外侧波束回波信噪比较低,波束展宽严重且容易受到中央波束“隧道效应”干扰,限制了外侧波束回波到达时间检测的有效性。针对此问题,国内外研究者主要从换能器基阵阵型设计和信号处理方法两个方面展开研究。

1 D" \7 ^+ l8 J# S5 Z

在换能器设计方面,通过阵型设计可以提高发射换能器外侧角度的发射响应,或者改善接收基阵外侧接收灵敏度,使换能器基阵对外侧回波的响应得到改善。也可以通过增加接收基阵阵元数目,减小接收波束宽度,改善外侧波束的测深精度。在信号处理方面,主要的研究趋势是研究分辨率更高的算法,提升深度测量精度,主要有3种技术途径:第1种是利用信号子空间类高分辨方法代替常规波束形成方法,使系统目标DOA分辨能力超过瑞利限,如多重信号特征法、子空间旋转法、解卷积类方法以及子空间拟合类方法等。第2种是利用相位法代替幅度法的波达时间估计方法,如多子阵幅度G相位联合检测法等。第3种是基于常规波束形成输出的拟合法算法,如BDI算法等。

- x# G0 S8 A+ j# k7 ?# e

⑵小水深测量

1 t( p9 [$ [- U5 I( w9 p5 ?

多波束测深系统不但需要对远处目标进行探测,同时还需要对小水深情况下的目标进行精细化探测,这就需要对近场环境下的目标回波按照球面波假设进行波束形成。虽然近场聚焦波束形成算法的基本原理比较简单,但其运算过程非常复杂,实时实现难度较高,浅水多波束实时动态聚焦方面的研究是小水深测量的基本技术保障。

3 ?3 v: I& W ~1 J' ]

⑶运动姿态稳定与补偿

- I1 D& C: [; R' V

多波束测深系统的载体在航行过程中不可避免地受到风浪的影响,因此对于载体的运动姿态稳定研究是提高多波束测量精度的另一个热门方向。载体运动过程中需要通过姿态传感器设备实时记录载体运动姿态,通过算法进行姿态补偿,从而得到高精细度的测量图像。需要通过波束形成技术控制发射波束和接收波束所对应的波束角度,包括接收横摇补偿、发射纵摇补偿、航行艏向补偿等技术。

! v, H }* p% x- l7 q( D) t

⑷精细化测量

4 E0 {8 c8 ]& i" A! c3 ]1 d. ~2 F

随着多波束测深技术的不断发展,研究者们希望通过多波束测深系统得到更为精细的测量结果,因此研究的方向主要集中于距离向精细测量、水平向精细测量、航迹向精确测量3个方面。距离向的精细度主要取决于系统的采样频率,采样频率越高则对回波到达时间的估计越精细,同时LFM信号的匹配滤波技术也能够提升信号的处理增益和时间分辨能力。水平向的精细测量主要取决于波束密度和波束宽度,更多的波束数目、更小的波束角度能够带来更精细的测绘条带。航迹向精确测量的局限性在于航速与帧率的制约以及多波束系统的航迹向波束脚印较宽,目标分辨能力不够,需要一种新的探测机理有效地提升系统的航迹向分辨率。结合了多波束测深技术和合成孔径技术的MbSAS技术,这种有效的途径近年来逐渐受到研究者们的关注。

% D% I, v8 X! u7 A

⒉ 合成孔径声呐研究进展

6 {' i, O! D8 b0 _- r# Q Y

合成孔经声呐技术的发展最早可以追溯到1967年美国Raython公司的Walsh等人,他们从1967年到1969年分别发表文章阐述他们把合成孔径技术应用到对海底小目标如锚雷等进行高分辨成像的研究结果。近些年来,合成孔径技术的发展已经由实验室走到了外场,更多的理论验证样机和海洋试验出现在学术界的视野内。目前主流的合成孔径声呐一般采用侧扫式合成孔径方法,国内外学者和声呐厂商纷纷推出各自的研究成果并推向实际应用。但是这些研究都没有很好地解决垂底区域存在缝隙问题,普遍需要单独使用多波束测深声呐或者成像声呐进行补隙,数据拼合效果有待提升。现阶段合成孔径声呐的研究热点主要集中在目标回波模拟、合成孔径成像算法、载体运动姿态补偿等方面。

$ J* T( F. A% z& D

⑴目标回波仿真

* H- ]# y( b7 V( a; W6 B

由于水下目标探测外场试验条件复杂,不可控因素多成本高,需要进行大量理论仿真研究,比如目标回波模拟以代替部分外场试验,然而目标回波模拟是一项相当复杂的工作。目前国内外很多专业机构已经展开了相关的研究并取得了相应的进展,例如北约水下研究中心的SIGMAS软件仿真系统以及新西兰Cantbury大学开展的掩埋目标回波研究等。

9 C8 R$ Z Y5 b5 v' Q

⑵合成孔径成像算法

$ d% F8 }. {- t5 \* F. w* r0 I

合成孔径成像算法的基本原理就是利用接收到的回波信号的时延信息求解出目标与收发换能器之间的距离,进而推导出目标的所在位置。常见的算法有:时域延时求和算法、距离多普勒算法、ChirpGScaling算法、波数域算法等。根据所使用基阵的阵型推导出各阵元与目标之间的时延差,并提出实用的成像算法是合成孔径技术的研究热点。

6 G) Y+ Q. o) ?. v6 O5 Y0 ^, p

⑶载体运动姿态补偿

4 W/ \$ K* q" q. s) t

如果想获得航迹向虚拟大孔径的分辨能力就需要非常准确的航迹向航行轨迹,而实际上载体航向的偏移等运动误差形式是一直存在的,这种载体的运动误差会造成图像的散焦,所以在合成孔径技术的研究中运动误差的估计与补偿是其实用化的最大瓶颈。相位梯度自聚焦算法(PGA)算法利用回波信号相位上存在的冗余度,理论上实现对任意误差的校正。多子阵SAS系统中可以采用冗余相位中心(DPCA)算法,通过重叠目标的相关处理获得相位误差信息。寻找有效并且价格相对低廉的载体多自由度运动误差估计和补偿方法是目前合成孔径技术的研究热点之一。

5 d& ?2 N1 L$ [+ O, d \

⒊ 多波束合成孔径声呐新技术

$ F) A" M, j' a3 y$ A

相比于侧扫合成孔径声呐而言多波束合成孔径声呐的研究起步较晚,最先见于文献的是2001年日本的研究人员在SeaBeam2000多波束测深声呐的基础上使用了合成孔径的算法,得到了很好的探测效果。2002年美国研究者向美国专利局申请了多波束合成孔径声呐的发明专利申请,在国际上首次提出了多波束合成孔径声呐的初步设想,然而其后,国际上未有该机构研究者利用多波束测深声呐进行合成孔径算法深入研究的文章公开发表。

& w) U$ p6 F( Z, t3 x/ @% {" i) N, c

2015年,Kongsberg公司首次利用该公司EM2040C浅水多波束测深系统数据,进行合成孔径算法处理,并将结果与多波束测深声呐结果进行对比。对比结果表明,经合成孔径算法处理后,能够得到更为精细的水下地形图像,该公司将这套系统称为HISAS 2040,这也是国外目前为止见到的最新利用多波束声呐数据进行合成孔径算法处理的实例。

( w2 r1 i) s1 v& \9 Z. B- u

哈尔滨工程大学通过对侧扫合成孔径声呐的研究后在国内率先提出多波束合成孔径声呐的概念,并独立开展了利用现有基于单线阵的国产多波束测深系统进行试验,证明了多波束合成孔径声呐的可行性,相较于传统多波束测深系统分辨率具有显著提高,并且能够一次测绘得到全覆盖测绘的结果,对目标的深度信息、航迹向坐标信息等有良好的成像效果,可以在保证与侧扫合成孔径声呐相同航迹向分辨率的前提下有效地提高合成孔径声呐的距离向分辨率并完成正下方无缝隙测绘,成像效果如图1所示。目前正在上述研究的基础上,开展基于二维面阵的多波束合成孔径探测机理研究,期望能够获得更好的航迹向分辨率和有效的提升系统探测效率。

! i l8 M- E( d: J

& ]8 z4 v' m8 T

图1单线阵多波束合成孔径声呐成像效果图

2 多波束合成孔径声呐基本原理

⒈ 多波束合成孔径声呐基本模型

. k# s6 ^) i1 U+ N9 M( P4 ? Y0 _

多波束测深声呐的基阵排布方式一般为接收阵元沿距离向依次直线排布,合成孔径声呐的收发阵元一般为沿着航迹向直线排列。为了解决侧扫式合成孔径声呐的不足,笔者融合了合成孔径声呐和多波束测深声呐的基本模型提出了一种多波束合成孔径声呐测量模型,能够一次性地完成测绘区的全覆盖测绘,不需要额外进行补隙,同时多波束合成孔径声呐能够通过距离向的波束形成,得到目标回波方向,从而解算出目标的深度,形成一种三维成像声呐,基本模型如图2所示。多波束合成孔径声呐与多波束测深声呐的最大区别是前者的发射波束沿航迹向的开角很大,这样在航迹向的不同位置波束会多次照射到目标,从而可以通过合成孔径提高航迹向的分辨能力。

0 z* A: U3 m5 [: R7 U2 l+ Q

- c. w @- |# E" M

图2多波束SAS的基本模型

1 \/ G5 Q8 v$ ^8 f! Z* c

⒉ 二维多子阵多波束合成孔径声呐阵列

( N3 g8 J; L# r

多波束测深声呐多采用Mill’s交叉的“T”型换能器结构,但是当多个目标的斜距相同时,栅瓣会使成像模糊,尤其在对大面积水底地形进行测量时影响显著。在实际的测量中,受到探测机理的限制,侧扫式合成孔径声呐测量效率将非常低。在SAS系统中,常采用的方法是在航向使用多个接收阵列,即多子阵SAS,可以有效地提高测量效率,如图3所示。因此,根据多波束合成孔径的原理和多子阵SAS结构,多波束合成孔径声呐换能器阵型选为平面阵结构,其距离向接收单元能够完成垂直于航行方向剖面内的波束形成,而其航迹向接收单元能够保证距离向波束合成处理时的栅瓣抑制。

. Q3 ]/ i. D8 x0 [+ B

|, S8 J& M0 V: h

图3二维多子阵多波束合成孔径声呐阵列

3 多波束合成孔径声呐关键技术

现阶段多波束测深技术的主要发展趋势是朝向小水深测量、运动姿态稳定、精细化测量等方向发展,而合成孔径声呐的研究热点主要集中在目标回波模拟、合成孔径成像算法、载体运动姿态补偿等方面。通过对比研究可以发现,二者研究的共同热点方向在于:利用复杂的探测信号形式,得到更好的回波信噪比和更精确的波达时间分辨力;利用更好的聚焦波束形成技术,实现对近距离目标的精细探测;对载体运动姿态进行有效的补偿,得到质量更高的声呐图像;利用新的探测机理,获得更高的目标分辨力。因此,针对多波束合成孔径声呐新机理开展了如下关键技术的研究。

. G8 J3 N8 U6 i$ h2 O& d. {

⒈ 合成孔径声呐目标仿真

$ a4 A( y8 }+ U9 ?7 g% q

目标回波仿真是多波束合成孔径声呐技术研究的基础,基于声呐运动模型的目标回波模型仿真结果对于多波束合成孔径声呐成像算法以及运动补偿算法的研究会有很大的促进作用。目标的三维仿真模型包括声呐的运动模型、目标的回波模型、目标的阴影模型等,利用抗干扰能力较强、距离分辨力高的线性调频信号进行目标探测研究。

8 R, L+ s6 a( W3 Q3 G

⑴声呐的运动模型

: A6 F; I8 Z, G; r3 ?+ X

多波束合成孔径声呐的基本运动模型如图2所示,接收基阵为由无指向性的阵元组成的换能器线阵。在工作过程中发射换能器以一定的重复间隔向海底发射脉冲信号,记载体运动方向为y轴方向,接收线阵所在的横向为x轴方向。计算得到声呐接收阵在不同时刻接收到的信号的时延,对于多个目标点可以分别得到各个点的回波,然后将回波进行叠加,即可得到目标信号仿真数据。

# A* ^% X6 W: `3 _7 l- _- `3 ~

⑵目标回波仿真

$ i8 U' l ~2 r

假设探测面目标或者体目标时,采用点目标重构法,将目标分解为一个个单独的点目标。以正方体为例,首先对目标体以固定间距进行切线分解,将目标分解成众多小块,然后利用体表面各条切线的交点来构造点目标,完成整个目标分解。将体目标分解为点目标,求出各个点目标的回波并进行叠加,可将叠加后的目标回波视为面目标或者体目标的回波,目标模型如图4所示。

' |8 P' E4 F5 U& K' M& U( ~

- |8 U% P4 x% Y" e$ l) w- k1 l

图4目标模型分解示意图

& x$ S$ B1 z0 H8 A0 t/ L5 K

⑶阴影区域仿真

6 j3 ], X7 S5 W$ k! o D

在声呐基阵沿着航迹向运动时,发射的声波束照射到物体,由于物体的遮挡会在物体的后方形成阴影,在接收基阵与仿真目标之间进行连线,延长线与水底相交,目标的着底点与延长线交点所围成的区域即为声波阴影区,目标模型如图5所示。

0 R2 W; z; p a4 g- {5 X

$ e- ^# ]7 `; `/ N( ~: h/ \

图5立方体目标的阴影区域

/ m3 e- i" X' f" ]' ~# L5 z

⒉ 多波束合成孔径声呐成像算法

6 ~. ^ _3 L* c$ @

多波束合成孔径声呐成像算法原理是将SAS逐点成像算法与波束形成算法相结合,经合成孔径技术处理后可得到目标所在的航迹向坐标和斜距两个物理量,波束形成技术在空间预成多个波束,将空间划分为不同的波束角度,根据目标的斜距和所在波束角度可计算出目标的深度信息,从而对目标实现三维成像。根据二维面阵结构推导目标到各接收阵元的距离,从而计算出各接收阵元接收到信号的时延差,阵元时延结构如图6所示。

' B8 T; k3 v3 y* Q5 @6 _) |' S% U$ @) w

6 B7 x& G- r" T

图6多波束合成孔径声呐各位置时延示意图

: h' F: n- L, r: R8 t

线性调频信号可以提高时间带宽积,通过对载频线性调制的方法使其频谱展宽。线性调频信号可以使声呐同时获得较大的作用距离和距离分辨率,同时其具有较大的抗干扰能力,因此多波束合成孔径声呐采用线性调频信号进行合成孔径声呐的探测,处理流程由正交变换、脉冲压缩、航迹向的合成孔径处理、距离向的波束形成4个部分组成。

4 c5 }; K8 x2 M& g* B8 C

⒊ 多波束合成孔径声呐联合运动误差估计

3 }$ ]3 }: p1 r9 o

多波束合成孔径声呐的载体在航行时受到风浪的影响,不可避免地会发生运动轨迹的偏移以及载体自身的摇摆,载体的运动失配将会造成图像的散焦。所以在多波束SAS的研究中运动误差的估计与补偿是其实用化的最大瓶颈,因此需要一种有效并且成本相对低廉的载体运动误差估计和补偿方法。

, V& t/ q& _2 E- R0 J9 `! U

多波束合成孔径声呐载体的运动可由6个自由度分别表示,各种单自由度运动估计算法都是根据相关函数估计出信号的时延,从而对各自由度的运动误差分别估算。但是目标到载体的距离是6个自由度共同作用的结果,如果各个自由度的误差分别做估计,将给运动估计带来较大的误差。因此,提出一种根据多个强点目标的回波数据的6个自由度的联合估计方法,同时对6个自由度的运动偏差做出估计,信号模型如图7所示。

6 C; E: m0 T8 L. }' |; \+ A5 O

" F. c" D2 s) Z+ V( Y2 l, z

图7多波束合成孔径载体运动失配

. T; l: ?3 M5 A' L4 R" b6 D

选择基阵的3个顶点上的阵元,记为A0(x0,y0,z0),B0(x0,y0,z0),C0(x0,y0,z0)。然后根据阵元到目标点的距离计算出O、P、Q三点坐标,根据基阵在S1位置接收到的回波和求解出的O、P、Q三个点的坐标来计算阵元A1,B1,C1的位置,即可得出此时载体的运动误差。对估计出的运动姿态误差进行算法补偿,即可得到更为清晰的声呐图像。

4 多波束合成孔径声呐试验研究

为了验证多波束合成孔径算法的有效性,开展了基于二维面阵的多波束合成孔径声呐试验研究,通过水池航车走航试验,进行了不同目标的探测试验并与常规多波束测深系统成像结果作出对比,试验系统结构图如图8所示。

9 Y2 \3 Z1 k! C

; Y1 o% w, B( ^5 o8 S2 e

图8多波束合成孔径声呐试验系统示意图

, c8 N/ Q+ ?0 V( H* B

首先进行了边长30cm方块目标的探测试验,观察目标成像结果可以发现,直接使用常规多波束成像算法时,由于受到波束脚印扩展的影响,目标尺寸发生了明显的增大,成像结果不能反映出探测目标的真实尺寸。经多波束合成孔径算法处理后,可以观察到探测目标的分辨率较常规多波束测深声呐有明显提高,并且目标回波强度也获得了显著的增强。其后又进行了直径13cm的双球目标探测实验,经多波束合成孔径算法处理后可以发现,双球目标能够明显的被区分开,可见该算法对相邻小目标也具有较好的分辨能力。由以上试验可以证明,多波束合成孔径算法能够有效地增强目标回波强度,提升目标的分辨能力,方块目标成像对比图如图9所示,双球目标成像效果如图10所示。

' y, T7 {; s; D' P, s

/ U4 M8 B: Q+ B/ l) B9 y% `* @

图9方块目标成像效果对比图

" K( p: f6 r# H! j# d r

: h0 l& M! F: n5 l. W3 V$ ?

图10双球目标成像效果图

5 结论与展望

多波束合成孔径声呐技术结合了多波束测深技术与合成孔径技术的优势,能够通过载体运动在航迹向上虚拟合成较大的基阵孔径,既可以在航迹向上获取较高的分辨率,用于对地形地貌的全覆盖测量,还可以在距离向上通过波束形成确定目标所处方位,并可以精确地测量出目标的深度信息,对目标进行三维成像。通过对多波束合成孔径声呐关键技术的研究和水池试验的多种目标探测试验,初步验证了多波束合成孔径探测机理的有效性。通过与常规多波束成像算法的对比,证明了多波束合成孔径技术具有提升水下目标分辨能力的潜力。

1 O/ m% P/ ]- u+ ~ y; a

多波束合成孔径技术的发展,紧随着多波束测深技术和合成孔径技术的发展趋势,代表了多波束测深声呐未来的一个重要发展方向,将会在水下小目标探测领域发挥重要的作用,是一种具有广阔应用前景的新颖探测技术。对于多波束合成孔径声呐技术的研究还有很多的工作需要开展,更快速的成像算法是多波束合成孔径声呐系统实用化的必要保障,寻找消除运动姿态估计累积误差的滤波方法能够有效地提高成像质量,复杂环境下的目标成像试验能够更真实地反映算法的有效性,这些都将是未来研究工作的重点方向。

; B* D4 J5 x8 V9 y

回复关键字,获取相关主题精选文章关键字:名家言论 | 政策法规 | 产业快讯 | 无人机 | 人工智能 | 王凤岭专栏 | 新材料 | 科技前沿应用 | 3D打印 | 创新科技 | 人才教育其他主题文章陆续整理中,敬请期待······

& R; `- D n' u: `1 O' p" u/ s

蓝海长青系列自媒体

& w' z1 \# b" }. R

% S7 a" Y) m( O" A

" @8 y, O1 t7 l

9 S% V" [5 g G+ h8 G0 K( ~ 1 D% k S ^: y- b. y/ G) g* C 5 l1 j& ~+ T$ W4 D6 P3 A/ V: g7 N8 V% Q; G2 N- i, `& p- w $ F" N% \( i2 T) l
回复

举报 使用道具

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
超级爸爸
活跃在2026-4-4
快速回复 返回顶部 返回列表