海洋技术▏高精度和高分辨率水下地形地貌探测技术综述

[复制链接]
* a3 i; C" I+ i$ O% E% U

《测绘学报》

6 F9 y, {# ~& V( q

构建与学术的桥梁 拉近与权威的距离

8 f- ~' F1 \; f$ A n% S

海洋蕴藏丰富的资源,对海洋的科学认知和高效开发利用有利于国家维护海洋权益和实现可持续发展,而这需要详细和精确的海底地形地貌信息作为保障。此外,河道安全、淤泥挖掘、水库监测和水下考古等也须获取水下地形地貌信息。因此,对高精度和高分辨率水下地形地貌信息的需求日益加大,相关探测技术的价值不断显现。

. o# {) o7 y( U) d3 q% c

随着声学尤其是多波束探测技术以及卫星遥感和计算机等技术的进步,彻底改变水下地形地貌探测方式,即采取数字化测量、自动化操作、智能化处理和可视化显示等综合技术,以各种平台为载体,在不同深度水域,快速获取高精度和高分辨率水下地形地貌信息。其中,探测平台包括卫星、飞机、船舶、水下机器人和水下观测网,探测深度包括浅层(水深小于200m)、浅中层(水深200~2000m)和深层(水深大于2000m)。

, Z" r8 J) R: ^$ _; H7 `

一、探测技术

' q1 l. Q0 l# d/ h/ c

⒈机载激光探测技术

! e- ~0 T& _/ J. H

与传统的船载声学探测技术相比,机载激光探测技术具有精度高、范围大、密度大、周期短、成本低和机动性高等特点,且具有浅层水域探测优势,探测深度最小可达0.15m,有效填补近海探测的空白。此外,最大探测深度为90m,精度通常小于0.3,测点密度最大可达0.12m×0.12m。

5 ?& f6 k" ?$ b& Y* j; P

自20世纪60年代起,中国、美国、加拿大、俄罗斯、法国、瑞典、荷兰和澳大利亚等国家相继开展机载激光探测技术研究,常见的机载激光探测系统如表1所示。

9 U- s+ ^- z' w) w9 |; J

表1 机载激光探测系统的主要参数

; ]) {6 E9 y `( `+ Q+ Y% s; I6 O

注:计算探测效率时,作业航速为70m/s,重叠探测覆盖率为200%。

$ b7 c" V H) I3 `% ^

机载激光探测技术在获取浅层水域尤其是近海的高精度和高分辨率水下地形地貌信息方面具有广阔的应用前景,目前我国在该领域仍未实现产业化,但已研发多代样机并开展多次探测实验,实验结果符合预期,发展潜力较大;但仍须突破若干关键技术,进一步提高探测精度,满足IHO S-44一类标准的基本要求和探测任务的实际需求。

+ J" q$ z" Y' }9 a

⑴定位和定向技术。机载激光探测技术对探测点的精确定位和定向通常通过高精度航空遥感运动补偿姿态系统(POS)实现,目前我国相关研究已取得突破性进展,如基于滤波器建模的时延补偿法的对称8点减振惯性敏感组件结构、采用伪距和伪距率的紧耦合导航系统以及适用于PPP/INS紧密组合的前后向平滑算法,都可有效提高POS的精度和稳定性。

+ \; N! r; D2 F1 M! u

⑵波浪和潮汐修正技术。黄谟涛等和胡善江等研究并总结无修正法、滤波法和惯导辅助修正法3种波浪和潮汐修正方法,利用平均海平面的高程计算探测点的海底高程,解决探测瞬时海面与深度基准面的相互转换问题,从而提高探测精度。

/ ?7 t- U# \9 [, N% m1 n

⑶探测精度评估技术。对探测数据误差的分析和校正是获取高精度和高分辨率水下地形地貌信息的关键,黄谟涛等和欧阳永忠等对深度归算技术的研究,为评估探测精度提供参考。

& {9 l, K: } k8 \, K

⑷多源数据融合技术。在处理探测数据的过程中,须对航线边缘的覆盖探测、不同时间的分段探测和多种设备的组合探测等的多源数据进行精度补偿,从而提高整体探测精度。

. t' U/ h* u0 w; d7 Z

⒉多波束探测技术

7 J! k/ x2 Y4 V5 Z! l ]- {" ]

多波束探测技术采用较窄波束的条带式探测和宽扇面的覆盖探测、高密度的探测点云、集成化的组合系统、结合传感器的姿态修正以及声速剖面的声线修正,具有探测范围大、速度快、效率高和精度高等显著优势,保障高精度和高分辨率水下地形地貌信息的获取。

' N. n/ \& ~; e/ J

根据波束形成的原理,多波束探测系统可分为电子多波束探测系统和相干多波束探测系统。与电子多波束探测系统接收多个波束相比,相干多波束探测系统通过回波角度的相位获取数据,具有数据快速采集和处理的特点。适用于不同水域深度的多波束探测系统如表2所示。

8 ~8 W* j- k2 e. H9 |% M

表2 多波束探测系统的主要参数

6 M" n; \& \. R! e

注:计算探测效率时,作业航速为5m/s,条带有效覆盖宽度分别为水深的7.5倍(浅层水域)、6倍(浅中层水域)和4倍(深层水域),重叠探测覆盖率为200%。

. m9 |2 `/ ?3 P- J* U- ^

多波束探测技术是获取高精度和高分辨率水下地形地貌信息的最重要的声学探测技术之一,基于应用需求,须从硬件设备和软件数据2个方面着手,取得突破性的研究进展和技术创新。其中,对采集数据的后期处理是至关重要的环节,涉及的关键技术包括声速剖面处理技术、异常数据检测和处理技术以及条带数据拼接处理技术等。

0 _- R$ o" J. i0 |

⑴声速剖面处理技术。声速剖面处理是多波束探测数据处理中的重要环节之一。由于海水具有非均质性,声线在海水不同介质层间发生折射而弯曲,不仅影响探测深度(声线的传播距离),而且影响探测结果,且越往条带边缘处声线的弯曲越明显,影响也越大。赵建虎等和丁继胜等提出采用实测声速、温度和盐度的剖面数据,建立声线跟踪模型和计算公式,从而减小声线弯曲造成的误差,提高探测精度和分辨率。

- Y" C. |* t& I# b! K( l8 i

⑵异常数据检测和处理技术。受海风和海浪等海洋环境、船身震动、环境噪声和气泡等因素的影响,采集的原始数据难免存在异常,对异常数据的检测和处理决定探测结果的真实性和可靠性。

" O+ `, I5 Z4 S5 h5 B/ J

此外,在数据量极大的情况下,人-机交互的数据处理方式可靠性较高但效率较低。目前CUBE算法取得极大的研究进展,并在异常数据检测和处理中得到广泛应用。在实际应用中,采用人-机交互和智能算法相结合的方式可获得更好的效果。

4 Q8 r$ l: @% k. O5 w. o

⑶条带数据拼接处理技术。为实现对水下地形地貌的全覆盖探测,相邻条带的边缘一般有所重叠。受随机误差和环境干扰等因素的影响,条带边缘重叠区的数据容易产生明显误差。因此,须建立相应的精度补偿模型,基于多源数据融合技术,采用误差验后补偿理论和两步平差计算方法,实现条带数据的合理拼接处理。

7 L% H+ s; c8 i

⒊侧扫声呐探测技术

( _% K) F7 ~% Z. q8 `" U! d

侧扫声呐探测技术根据回声探测仪的工作原理,利用声线的反向散射信号获取水下地形地貌信息,并生成直观反映水下微地形地貌的图像,同时还可根据海底回波的强度定性分析海底底质成分。

* O, T( Q2 O& f: c! Y) v. g) I

侧扫声呐可分为船载式和拖曳式,其中船载式侧扫声呐的探测幅度宽和效率高,但分辨率较低;拖曳式侧扫声呐可根据探测任务的需求调节拖体距水底的高度,以获取高精度和高分辨率水下地形地貌信息。

2 @2 x. T; G# G6 s! N9 r* z

传统的电子型侧扫声呐的正下方存在探测盲区,而相干型侧扫声呐结合条带探测手段,可通过声线回波的相位精确定位探测点,并可提高探测精度。但水深大于200m时,相干型侧扫声呐的探测效果不如电子型侧扫声呐。

$ [# p5 A6 S" y+ @ h/ s( F/ t' d8 X9 j

侧扫声呐探测的最大有效作用距离取决于声呐设备的工作频率,探测分辨率则与脉冲宽度(距离分辨率)和水平波束开角(水平分辨率)有关。常见的侧扫声呐探测系统如表3所示(不同探测频率条件下)。

6 v, b' n" v' d# X

表3 侧扫声呐探测系统的主要参数

1 ~( K/ }, z1 h6 r' h
/ V3 `6 g# g, p# \* r7 K2 ^3 l
6 x& g% t: L. b1 ~; B* z( b

系统

, m- k, A5 i8 R, s: w2 L

EdgeTech 4125

9 J, z. {0 o3 U8 `

GeoSwathPlus

k1 Z: T- U/ @6 Z# b/ S6 @& W

声呐类型

# ]! c2 Q3 X2 l9 E) z+ s

电子型

) l5 Z3 ` h; r5 i7 m2 f Z w- ?

相干型

: ^- C! d) L1 p2 i* A9 h* E

探测频率/KHz

# I: F4 O- a2 b! E! V! {; v1 J

400/1600

1 R0 F. g, p" i1 J2 K6 K$ c1 t

125/500

: M9 n: F k% V. e& r

最大探测深度/m

( `) U$ \" Q1 I, O

200

' e' k3 {* K- ^* i: I4 q

200/50

5 ^) @: d: |( l1 @# E; ?* h) q

最大条带宽度/m

8 @" F8 G+ Y3 k# N; j8 N

150-35

( ^1 }2 S" X9 L% R/ i% O3 Q5 a

600/150

# `$ ~8 y S' Z, E

探测精度/cm

- V; |5 ]; J# k4 v7 O

2.3-0.6

& s$ b2 i& v5 k" G) e% b

2.4-0.6

1 l7 |, P) j9 ?1 I/ t

探测效率/km2·h-1

3 T: D- |0 Z/ m( [

2.70-0.64

" |3 m6 D; v) A- V. _

10.08-2.7

6 h9 I+ g9 X1 I
* N& `* M3 B1 ]: {( F6 r9 ]
$ w* c0 d+ q5 t0 f/ h" a

注:计算探测效率时,作业航速为5m/s,全覆盖探测时只须测线采用2倍有效作用距离,无须设计重叠带

2 K/ ~% J: |4 f" ^1 E, I7 S

侧扫声呐探测技术具有造价较低、分辨率高和可连续获取二维海底图像等优势,得到广泛应用。拖曳式侧扫声呐的拖体须尽可能地贴近海底,因而探测效率较低;同时,为保障拖体的安全,作业航速不能过快,进一步限制探测速度。

: b& X& U. j8 Y# V5 G

⒋浅地层剖面探测技术

( _4 f) s5 {4 T/ Y& u

浅地层剖面探测技术基于水声学原理,利用声线反射,连续走航式探测海底地形地貌、浅层沉积地层和基底、特殊目标以及浅表层灾害地质体,具有效率高、分辨率高和成本低等优势,且探测所得的声学记录剖面在形态上与真实地质剖面极其相近。

4 N$ J% L& P, l

随着非线性声线传播技术和实时海底追踪时延记录技术等的发展,浅地层剖面探测技术实现高精度和高分辨率水下地形地貌探测,并应用到深层水域。常见的浅地层剖面探测系统如表4所示。

8 i) o8 U7 Q$ {9 i* q2 U

表4 浅地层剖面探测系统的主要参数

0 V0 E/ v! d. }- E5 j4 ^6 p

海底底质、环境噪声和船身震动等都会干扰浅地层剖面探测,为保障获取高精度和高分辨率水下地形地貌信息,须减小甚至消除这些因素的影响,可采用直达波、多次波和侧反射等方法进行分析和校正,并改变声源和接收基阵的指向等。此外,根据应用需求,浅地层剖面探测技术也应有所突破,如:采用超宽频技术,解决分辨率与地层穿透深度之间的冲突,提高在深层水域的探测能力;与多波束探测技术和侧扫声呐探测技术相结合,实现三维探测;进一步提高探测精度和效率。

. g7 R& g$ N' C; v

⒌双频识别声呐探测技术

a T; Y$ n6 x1 ~8 T

双频识别声呐具有较高的分辨率,利用声学透镜折射率大和聚焦短的特点形成较窄波束,最大限度地减小传播损失,可在能见度极低的水下生成接近光学照片的高质量图像,相当于高性能水下摄影机。常见的双频识别声呐探测系统如表5所示。

" M9 f6 \5 l+ Y# U% t

表5 双频识别声呐探测系统的主要参数

5 _: u G! @1 s2 M; a) j$ f/ F
4 {8 e. G+ _2 h" Q4 z5 u& A- X$ _
4 i: H0 Z6 W, \6 [6 _0 c6 n

系统

! g- @' d" q8 A2 a

DIDSON标准型

7 V+ J9 f0 q. H" ?2 [+ Z f

ARIS 3000

% M$ t7 L% K# P2 Z* d- b' O

低频

/ T4 k3 X; X7 i" W# J

高频

7 g( N# _& ?0 e

低频

p: |( y9 m7 t( C+ s; W) m: `

高频

! Y' C7 I. L5 p

探测频率/MHz

7 Y1 `. k# G$ ^1 i8 G& q

1.1

+ V/ }+ _1 S4 @& ]8 b! h+ L

1.8

$ z3 V5 y9 s9 ?3 g% Q2 r

1.8

% C# g/ K- u4 b! o# `

3.0

6 y# A. u6 T; W- b" R9 O

最大探测深度/m

. m, s+ L1 R) I d6 n

300

- f T) ]6 _9 G% H5 g, H0 s

300

! Q' V8 J1 {5 U, k$ F4 K+ B* q

300

/ H+ y$ o5 s+ Y9 [8 p

4000

9 U0 ^4 i/ u0 N* {

探测距离/m

, M1 P9 c; N+ j6 Y Y g, ^8 F! T

4.5,9,18,36

$ l. ^: Q6 t( [

1.13,2.25,4.5,9

/ Q" Z: L8 [) c4 ~

15

# q* {% t5 o- K+ |: I8 N; v! [

5

7 W8 f8 S8 s2 e5 a! s

取样间隔/mm

3 w; | a0 k4 _& u

8,17,35,70

! s( U$ M/ z2 o5 |0 B

2.2,4.4,8,18

4 m, r4 P( u. V2 Q9 v" s c

3

0 W# ~" N8 a0 @" _7 _7 X8 g0 N

3

. P1 d! }) ~; R6 a' B

水平视角/°

$ n1 F. A% H0 I

29

" C+ E) b) N, U: q, G! `

29

9 F' M5 u' z" v

30

8 z9 @7 c) W3 v) d

30

/ J% ?* m/ |+ J$ ^/ W

最大探测效率/km2·h-1

. R9 z4 V: h8 g+ R. }8 E4 f

0.16

0 U8 b: V y: W" k+ S) J l, k

0.04

! x! L7 D& f C

0.06

$ X5 Q; b1 {* ^( g5 G! V( P1 q6 b

0.02

0 I9 P p7 |! E) O' ~4 w1 p
" j4 T; \6 k" \3 o0 E# b- l' H
9 Z& n! _! h) V* l9 A, {

双频识别声呐探测系统在获取高精度和高分辨率水下地形地貌信息的过程中,对声呐信号的处理和分析尤为重要,其影响因素主要包括海水介质、海底底质、海面影响、环境噪声和设备条件等。除尽可能地消除这些干扰外,对声呐数据的后期处理也是关键环节,可采用特征图像配准法实现对图像的无缝拼接,使图像直观、清晰和完整。

4 x' g' p6 ?! r% l, c0 E9 I* d

受声呐有效作用距离和作业航速的影响,双频识别声呐探测技术适用于对确定目标的精确探测,而不适用于大范围和较粗略的目标探测。

9 Z( _* {2 e8 l j A- g. t2 K

⒍合成孔径声呐探测技术

. q {, |1 V' g' L( \

合成孔径声呐探测技术利用小孔径的声呐换能器阵,通过运动形成虚拟大孔径,经延时补偿生成较均匀的高分辨率图像,适用于对水下地形地貌的全覆盖探测,且可提取确定目标的精确信息和完成三维成像。由于合成孔径声呐的分辨率与探测频率和距离无关,可比常规声呐高出1~2个数量级。常见的合成孔径声呐探测系统如表6所示。

7 y! K7 w. F' X# V0 Z

表6 合成孔径声呐探测系统的主要参数

) i& g! E4 C( l& u: F& N
0 [; W" {7 L! L; o$ L
$ x& `6 U: A2 j4 e% h$ m

系统

( j5 ?- s4 s! E6 U5 t; V/ O& L A

探测中心频率/KHz

. \ o( F, |4 P$ C: \

航向分辨率/cm

$ G1 I# Q8 W2 ?5 ]1 y N; S8 [

距离分辨率/cm

* t: @: F1 x- N- k

典型探测距离/m

8 R; l) s: A; c D. a, z% S. p, b

UCSB

' B/ s3 @* l3 S; D# c0 T7 z8 z! V( r

600

7 h. |$ g- J6 A# i, l+ ^5 ~5 \

5

# A1 N% V" @; L

19

% f2 I3 z: W+ S* h v! Y7 `$ [/ i

70

3 A; l1 m1 H2 T3 z2 e

HISAS1030

! O, _2 w9 a! w& U5 H3 v e- S

90

1 K: z5 P+ |5 [/ P/ ]* Q6 P8 V

5~50

, Z T! r* r5 z. _9 {) \

4

# u; T$ Q# p3 l# W3 J

200

. B s' N x" k3 g# P

DAPPA

. T0 Z! t/ f* w$ ]

50

& e5 s" T |6 o" r, \

10

* C! v$ H. i+ ~# ^ c5 _3 g( l

7.5

/ e9 B: v4 B. g- h

50~1000

8 s _' ^! B2 g

SAMI

& W2 \! u8 c; @6 a0 R% { l D: S

8

+ ?; Y# Z# _6 S% h( Y7 r: n

100

: N6 q$ z7 S' b; v: q2 f5 [

13

' ]- \. k7 P9 }. b5 ]& Z" a7 L

100~2500

% S- j% o$ S0 ~( |0 R
/ j+ V8 R# ~* ~/ F& T' v0 _
1 w# ~; t" U Y* e* z' {

中国科学院声学研究所于2000年成功研发国内第一代合成孔径声呐。李海森等研究多波束合成孔径声呐的机理,结合多波束探测技术和合成孔径成像技术,很好地解决了侧扫声呐探测精度较低和存在探测盲区等问题。

8 W2 ^6 J% w3 R9 u

为获取高精度和高分辨率水下地形地貌信息,合成孔径声呐探测技术仍须改进,如:对于水声信道稳定性和时空相关性以及物体散射性等造成的干扰,可通过对目标回波信号的相位补偿减小或消除;对于载体平台,可通过惯性导航设备的姿态补偿以减小误差;对于较低的作业航速,可采用多重接收基阵提高探测效率;对于数据后期处理,须完善多子阵成像算法、运动补偿模型和自聚焦算法。

% |3 N% b: N' h! i' m" V

⒎水下三维扫描声呐探测技术

7 }" q& C0 }# f! S! H$ m5 t6 u

水下三维扫描声呐探测系统类似于三维激光扫描探测系统,但其利用声呐设备发射声线,通过对回波信号的接收和处理,实现三维成像和定位。

. j8 I& J. B2 E. S' c- k

与常规声呐相比,水下三维扫描声呐可在水下实时探测,并生成较完整、精确、清晰和立体的高分辨率三维图像;而与三维激光扫描相比,水下三维扫描不受水体能见度的影响,即可在能见度极低的情况下生成高质量图像。常见的水下三维扫描声呐探测系统如表7所示。

' `/ X$ D8 y, D( Y. F0 S

表7 水下三维扫描声呐探测系统的主要参数

" b! z2 d, o6 g8 w) @0 s0 F

注:最大探测深度分别是探测系统在不同供电电压(120~240VAC和24VDC)下的参数值

& W% d W1 _, W0 L

通过水下三维扫描声呐探头的转动,可实现水平方向360°和竖直方向130°的立体扫描探测,获取高精度和高分辨率的水下地形地貌信息,但须关注声呐设备的固定和声速剖面的精度补偿。

: e6 a+ ^6 T& R, I! T

二、探测平台

) J: @( v: f$ p: i: u

常规的船载式探测不能最大限度地发挥探测设备的优势和性能,在中深层水域探测中的劣势尤其明显。为解决这一问题,一系列探测平台不断涌现。

1 A) w- U |8 q/ r

⒈无人船

. \. V8 C+ s1 K

载人船通常不便进入浅层水域和岛礁区水域,而无人船可通过搭载单波束探测仪、浅水多波束探测仪和侧扫声呐等设备,与母船协同作业,高效完成探测任务。无人船的探测效果取决于其搭载设备的性能,通过对数据的后期处理,获取高精度和高分辨率水下地形地貌信息。

; ]2 S$ F! q6 b! C: I" h

随着定位、导航、控制和自动化等技术的快速发展,无人船可通过遥控或自主的方式,连续走航式探测水下地形地貌。由于无人船具备自动GPS导航和自动规避障碍物等功能,可提前设定或随时接收探测任务,随后自动完成任务,极大地保障工作人员的安全和效率。

$ P2 }! V6 z/ z

⒉水下机器人

1 P3 G2 T+ E, B% z

水下机器人是水下地形地貌探测的重要平台,具有中深层水域探测优势,主要分为无缆自治潜水器(AUV)、载人潜水器(HOV)和缆控潜水器(ROV)。

6 f" u7 L: a# \. _

水下机器人可作为搭载各种设备的移动平台,获取高精度和高分辨率水下地形地貌信息,其性能主要体现在下潜深度、续航时间和航行速度等方面。目前水下机器人的最大下潜深度可达6000米级(AUV)、7000米级(HOV)和11000米级(ROV),提高续航能力是关键技术之一。

3 ~) z: g5 E) Q2 ~0 O

⒊海底观测网

/ [! v# K; c# J+ b5 l& l% B' V

通过在海底铺设以光缆和接驳盒连接的传感器,形成海底观测网,从而实时、在线和长期探测海底地形地貌和监测海洋环境。其中,采用光-电复合缆连接海底观测设备和陆地基站的海底观测网,适用于离岸较近的海域,具有数据传输较快的优势,但造价较高;无缆连接的海底观测网通过电池供电,并通过海面浮标发送卫星数据,适用于离岸较远的海域,且造价较低,但电力供应和数据传输受限较大。

, i/ M8 a8 s: O' W9 B

海底观测网可为其他探测技术提供基站,如为水下机器人提供中转平台和蓄能服务等,实现大范围、长期性和立体化的海底地形地貌探测。

8 f8 W @% z4 o( w; u5 Y

三、结语

9 Q/ @. G# X& `2 d2 f5 i8 s4 n

为满足军用或民用需求,开展高精度和高分辨率水下地形地貌探测研究具有重要的战略意义和广阔的应用前景。未来应重点发展多维度、多设备、多来源和多手段的探测技术,尽可能真实、完整、全局和精确地反映水下地形地貌,为进一步的研究和应用提供可靠的数据支持。

- E4 h! B1 N$ x) z. C& K

【作者简介】本文作者/宋帅 周勇 张坤鹏 范孝忠,分别来自32033部队和北京航天宏图信息技术股份有限公司。第一作者宋帅,32033部队,高级工程师,博士,研究方向为大气海洋和地理测绘;通讯作者周勇,北京航天宏图信息技术股份有限公司,工程师,硕士,研究方向为海洋水文数据处理。文章来自《海洋开发与管理》(2019年第6期),参考文献略,用于学习与交流,版权归作者及出版社共同拥有,图片系我平台加载,转载也请备注由“溪流之海洋人生”微信公众平台整理。

- V, ^$ ]$ h! l6 O

来源于:溪流之海洋人生

2 R* @2 t2 Q( [7 i( _

权威 | 专业 | 学术 | 前沿

- `' ^7 m, P4 p& Q: j6 q

微信投稿邮箱 | song_qi_fan@163.com

6 F' e0 z. s3 {& [# n

欢迎加入《测绘学报》作者QQ群: 297834524

" w; H0 J* M6 k$ i" F+ }& O

进群请备注:姓名+单位+稿件编号

, M; {+ A* I7 F; e0 L. x2 K+ w + D2 l4 A; P: n ( N6 A4 C- m( S/ [- O ( d/ [: [! S8 w! [. h8 j F* w, V% e0 h3 ]1 ?/ j& e" u; L
回复

举报 使用道具

相关帖子

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