 点击上方“溪流之海洋人生”即可订阅哦乔正明1，郭晓明1，商少平2，刘轲2，赵小峰3
: A, k" S% ^8 T5 \ 1.国家海洋技术中心漳州基地筹建办公室；2.厦门大学海洋与地球学院；
7 `! s0 ?: \/ t$ N! r 3.国防科技大学气象海洋学院
' ]. U; `9 N; V9 @! f/ U 【摘要】为满足当前ADCP计量检定的迫切需求，根据ADCP在水池拖车试验检定时换能器发射声波的主瓣、旁瓣与检定水池池壁、池底形成的几何关系，推导出了抑制旁瓣、不抑制旁瓣时检定水池的尺寸需求。以市场占有率最高的ADCP产品为研究对象，计算得出了不同工作频率ADCP对检定水池宽度和深度的理论值。本文导出的公式，进一步完善了ADCP计量检定的理论基础，为ADCP在室内可控环境下的检定规范制定及检定平台设施建设提供参考与借鉴。 * U" t9 H7 \8 M( }/ e/ Q/ f Z( F" b1 {
【关键词】ADCP检定；水池拖车试验；池体计算；换能器；旁瓣影响  
3 Y6 v5 e& S3 v4 f 一、引言 
0 F3 ~/ ^/ T6 q8 K 声学多普勒流速剖面仪(Acoustic Doppler Current Profiler，ADCP)，是目前国际上测量多层剖面海流及河道流量的最有效设备[1]，广泛应用于海洋调查和水文监测领域。尽管近年来ADCP发展迅速且应用广泛，但国内外对其检定还处于起步阶段，目前行业内尚无被广泛认可的ADCP检定校准或相关规范[2]。
0 y0 P& o1 a$ e; S z3 C 目前，ADCP流速流向检测检定方法主要有水池拖车、同步比测和自身航行试验[3]。经前人试验验证，在室内可控环境下的水池拖车试验是目前检定ADCP最可靠最有效的方法之一[2,4-6]。2010年颁布实施的国家标准《GB/T 24558-2009 声学多普勒流速剖面仪》[3]（下文简称国标）规定：工作频率大于300kHz的ADCP应进行水池拖车试验，且水池长度应大于100m，宽度大于7m，深度大于4.5m。国际上，对中高频ADCP（500kHz以上）水池拖车试验所用水池宽度都在15m以上[7]。综合国内外的各项标准及应用可见，均未给出不同频率ADCP检定所需水池物理尺寸上的最佳值。
) r# H" z, ~" l5 M. ] 本文参照国外相关研究，以NOAA等相关机构提出的检定方案[8]为基础，以RDI公司3000m型Workhorse Long Ranger 75kHz（下文简称WHLR-75K）ADCP产品为研究对象，针对ADCP水池拖车试验中所需水池的深度和宽度需求提出了新的计算方法，进一步完善了ADCP计量检定的理论基础，为今后ADCP在室内可控环境下的水池拖车试验检定规范制定及检定平台设施建设提供参考与借鉴。 . ~' B, O( A; J6 ], ^
二、ADCP简介
, |5 e) O! q; ]3 _5 e ADCP由换能器阵、温度传感器、罗经、倾斜传感器、壳体及应用软件等组成，一般分为自容式、直读式、船用式（走航式）[3]。为测量三维流速，ADCP一般配置2个以上不同方向且彼此独立的声学换能器，最典型的4种换能器阵型结构：3波速Convex结构、4波束Cyclops结构、4波束Janus结构和5波束Janus结构[9]。每个换能器都与ADCP的轴线成一定夹角，且都具有发射和接收功能，工作时换能器发射声脉冲波信号，通过随水体流动的泥沙颗粒、浮游生物等散射体反射声波信号，由换能器接受反射的声波信号，经测定其多普勒频移而计算出不同层深的水体流速矢量[9]。ADCP具有不扰动流场、测验用时短、测速范围大等特点。 3 C- e9 F* f! `- Y- y0 ?6 h, H1 K3 T
WHLR-75K ADCP为4波束Janus结构，相邻换能器法线间夹角90°，每个换能器法线与ADCP纵轴夹角20°，是一款用于长期、大量程海流观测的产品，其采集的年度或季度的详细海流数据，主要用于环境监测浮标或潜标、海上石油钻井平台监测、极地或大洋锚系观测等海洋科学研究或海洋油气区海流观测。 ; e' H6 q" D1 i& g3 \; H$ P6 J6 ]; F3 B
三、水池拖车试验 4 c1 ~ k1 ]5 w) C* {
⒈试验原理
: ]8 ~: B/ E0 f0 w+ h# L5 I 水池拖车试验主要由大型水池、运动拖车、拖车车速测量设备、ADCP安装架、回转机构、测量数据同步采集处理装置和被测ADCP仪器等组成[3,10]。试验时将ADCP固定在拖车上，并向水体中播撒石灰粉粉末等颗粒物作为水中散射体，采用水静止而ADCP运动的方式形成ADCP与水的相对运动，以ADCP自身的运动速度值检验其流速测量值，以ADCP在其自身坐标系上的运动方向值检验其流向测量值[3,10]。
2 I& m. ]" L( N; {' G# k ⒉试验水池建设要求 , U! d! Z. [1 ~/ u, C0 K( i
由声学特性可知，ADCP的工作频率越低，在进行计量检定时要求水池长度、宽度和深度就越大。检定水池需满足以下条件：一是水池要足够长，否则在检定过程中，ADCP运动的速度还没有达到稳定或还未获取足够的采样数据就已到达水池端头；二是水池要足够深，不同型号ADCP的盲区大小各不相同，水池深度至少要大于盲区，且在声波信号到达池底前，应获取至少2个完整且质量较好的深度单元流速数据；三是水池要足够宽，须充分考虑池壁对旁瓣的影响[2]。 . Q+ p' ^' K5 z3 t5 f m( ?: N
四、试验水池的物理尺寸计算 1 i0 r: T/ S/ D% x& f1 q
⒈水池长度分析
% w5 L; i( X# B( y5 Z 水池长度分析分为两个部分：一是根据具体试验需求计算水池长度；二是通过任务量分析选取合适的水池长度。
8 a. r; \' [, n* `# t9 C ⑴以具体试验需求计算水池长度 % T2 F" M- g) j, n
参考NOAA及RDI公司的检测检定方案，水池长度的测算公式为：
6 ? b8 D! P5 V" g" d6 y* X9 u" {* W6 _" h 水池长度＝加速段长＋稳定段长＋采样段长＋制动段长＋安全段长
% C+ u1 U# O% `0 x 其中： 6 ?; e7 E3 ^& T& `
加速段长=拖车加速到采样段速度所走长度；
0 n9 }5 x( D6 ]1 M 稳定段长=拖车采样段速度×稳定时间； ) X1 G# V: y% t# Q
采样段长=拖车采样段速度×采样时间； , M% U0 w5 b' b4 p
制动段长=拖车减速至零所走长度；
/ o2 T Y- ? n$ o7 i; B- ? 安全段长=一般设置拖车滑道和液气缓冲装置，可按20m选取； ; ~ e: y* x6 k8 G2 S, i
采样时间=采样剖面数目×每个剖面脉冲数目/采样频率； ' ~9 C( M5 A2 _ x! c
各参数取值：
6 d. w* e( \/ K" m2 l 拖车采样段速度：国标规定检定流速的最大值为4m/s[3]。 p- l5 c, ]) R" A
拖车加速和制动的加速度：0.65m/s2（根据目前拖车试验时所能实现的理想加速度）；
" ~1 `) ?' |, c& R) i 到达采样速度后的稳定时间：2s（根据实际检定情况具体确定）； 0 {1 s, ?7 B! h. T" L
每个剖面脉冲数：为保证检定精度一般设为10个[3]；
0 G$ @( P* t$ l 采样剖面数：每个拖曳速度至少应有30个剖面流速测量值，国标允许多次重复拖曳来完成30个流速测量值[3]； 6 U. _8 @+ m8 L( \
采样频率：1Hz（以RDI公司WHLR-75K ADCP为例）。 , x/ g" c1 L7 W$ g# \0 a# |! Z
因ADCP的工作频率越低，最大采样频率越低，对水池长度要求也就越大。按上述计算方法，若要满足30个剖面流速测采样的检定要求，在拖曳速度分别为1m/s、2m/s、3m/s、4m/s时，对应水池长度应分别不小于324m、631m、940m、1253m（见表1）。 / c% D6 c: t3 ?
表1 不同拖车速度对应水池长度计算表 0 c# T/ \( H* b5 U
拖车速度（m/s）
. w* x4 X) d- Y7 x" g 加速段长（m） : R( M1 x4 x0 W% {' V; ^+ a1 E$ Z
稳定段长（m） + f! B% A9 w/ W2 b: ?9 P0 Z
采样段长（m）
: o$ s9 q& ?+ i% p0 K/ @- ? 制动段长（m）
; y! r/ ~! u u& g% W 安全段长（m） ! [+ s4 W, [. e
水池长（m） 9 c( a, ^/ f9 ~4 g5 l" I$ l4 S$ W5 [
1
3 t V& d7 }7 b0 [# d4 B* t 0.77
1 O! D; x+ u2 j. q( B 2 5 A! @' m/ j" A' J3 }# W
300
* [9 u! }3 Y7 d+ {6 A4 @; n% _" H- ^ 0.77 & a1 W5 B; Z& |
20
9 I4 w0 b# K% L" B 323.54 . P& ? z& ^" ^$ T2 [+ @
2 5 b7 I0 I0 F2 z& R. i+ e% x
3.08 - ~- m% ?* @3 h0 w9 L n
4
R/ T) D4 b; H: r, S+ e# f3 h4 h% V) ~ 600 + a8 p; F8 k: L) `; X2 o9 z2 W! u3 p7 c
3.08 6 }5 o# J. r( T! v c
20
1 [6 j" U+ j! S s' K$ r9 l# x 630.15
" C, |7 J/ u+ s, U& _# F% K 3
+ H( K& ~* J6 O5 ^, W6 J1 v3 y9 s 6.92
6 V1 [: `2 `- G, t+ w2 i 6
. |9 r4 Y7 _+ h5 I 900
6 Z3 p- j( P: ?% X1 U 6.92
+ l0 n$ c& D8 u, [2 ?: K+ A 20
3 ]+ W$ G! I5 L0 y( o& X 939.85
6 W" q, P- `4 D1 q$ u! \, t 4
" B7 r8 z% p( o! {# \ 12.31 # b- J5 T) N; s" j
8 + V7 Y$ m% M, ~, l V5 F8 y
1200 1 b) v/ ?( n5 f. i; x4 F
12.31
9 O! _5 ?5 J! V. \, j; u 20
; I& H* ~! a9 p0 z4 s: O 1252.62
$ @+ V4 t( ]' s! j( f3 ?* M; n& } ⑵以任务量对水池长度进行分析 ! V2 N) b/ J6 f7 \7 O ^5 Z& m
因国标允许多次重复采样完成检定，水池越长重复采样的次数就越少，检定时间就会越短，因此水池长度还要根据具体年度检定任务量进行合理推算[11]，本文不再赘述。
: S5 t1 _8 A; G {& M ⒉存在旁瓣影响时所需水池深度计算
2 \' X' p- O% ], [ 如图1，depth为水池深度，width为水池宽度，h为仪器入水深度，d为仪器直径，D为仪器至池底距离，θ为波束角，φ为换能器主旁瓣所在平面与池底面的交线和池底边夹角，OF⊥GF，OH⊥GH，D＝OA＝OG。 & U/ m8 z5 N/ k; z# _) v
为确保ADCP计量检定效果，应获取至少2个完整且质量较好的深度单元流速数据。水池深度的具体计算公式为：
- Z' o( |) p3 j, g( V( S, H depth＝h＋h1＋n×c＋h2 ⑴ * l8 j9 i! {' N1 r9 @/ t1 t5 Y
式中，h为仪器入水深度，取值0.5m；h1为第一层深度单元中心至换能器距离，可通过该型ADCP自带仪器软件获得，取值12.71m；c为最小深度单元长度，取值4m；取2个完整的流速数据，此处n＝2，h2为旁瓣影响区长度，即h2＝HA＝OA(1－cosθ)，其中θ＝20°。  3 X% A' k/ k% `$ d
图1 水池宽度和深度计算用示意图（考虑旁瓣） 8 Q) r6 q: x* y/ e3 \) y6 Q
由图2可得出： 6 f/ |5 X/ y' h5 \ s# T3 _
h2＝OH/cosθ－OH ⑵
$ h4 p! ^2 ^" y& | 将OH＝h1＋n×c带入式(2)，得： 8 T/ h' {: r; j* A4 j6 W
h2＝(h1＋n×c)/cosθ－(h1＋n×c) ⑶
# X0 T% I2 [& z2 m. O% r6 S' ?' \" K 将式⑶带入式⑴，并带入相关参数可得到水池最小深度depth为22.54m。 
' ]1 j4 L7 D0 H+ \7 j( C 图2 检定安装位置示意图 2 d" |/ E; ^4 y3 Z
⒊存在第一旁瓣影响时所需的水池宽度计算
2 \7 ]8 u& A8 t/ j& V ADCP计量检定对水池宽度主要有两个影响因素：池壁旁瓣影响、仪器换能器与池壁的夹角。由图1可以看出，为完全避免声束主瓣到水平方向之间旁瓣的影响，水平方向上换能器到池壁的距离OC必须大于等于OA的最小值（即水池最小深度-仪器入水深度）。检定时，ADCP在水池中安装位置见图2，图1和图2中φ相同。由图2可知，水池最小宽度计算公式为： 5 o3 u2 q" c4 v% }8 A" h/ e* D( K
width＝2×w＋d ⑷
) g! M* M/ L D& J- v 式中，w等于OG×sinφ与OG×sin(90－φ) 中的较小值，d取RDI型ADCP最大直径0.55m，可证明当φ＝45°，OC＝OA＝depth－h时，水池宽度最小。代入式⑷，可得水池最小宽度width为31.71m。 : Z3 i; e% Z3 O( O, J/ G9 O
⒋讨论 4 W$ E6 z5 A& X2 J# R* x
在检定ADCP时，若能在4个换能器端安装导流吸声罩来抑制旁瓣对池壁的影响，理论上可减小对池体宽度的需求。下面分析旁瓣被抑制后水池宽度需求，本文暂不讨论加装导流吸声罩抑制旁瓣的实现方法。  7 i6 r4 `. r; M: G3 e4 u/ s }
图3 水池宽度和水池深度计算用图（抑制水平旁瓣）
/ v+ R% K( F0 S, G9 y 由图3可知，若旁瓣得到抑制后，OC′为主波束的边缘线，为了完全避免池壁的影响，OC′必须大于等于OA的最小值（即水池最小深度-仪器入水深度），OC′与OA之间的角度为（θ＋β），β为一半的波束宽度，取值2°，水池最小宽度计算公式为：
* K0 y' I6 {& Q5 {$ ?5 U- h; b% c/ C6 a width＝2×OC′×sin（θ＋β）×sinφ＋d ⑸
: ~7 F4 Q0 @) t- I- v2 n3 E0 w" }! J 当旁瓣被抑制后，只考虑主瓣对水池宽度的需求时，可证明当OC′＝OA＝depth－h时，且φ＝45°时，水池宽度最小，代入⑸式可得水池最小宽度width为12.23m。
% E6 f, C5 ^% E/ B; T+ \* p 综上分析，为了满足WHLR-75K ADCP的计量检定需求，水池深度应不小于22.54m，不抑制旁瓣与抑制旁瓣时水池宽度分别应不小于31.71m、12.23m。 * [+ s9 E, w2 e8 ?' K
⒌不同工作频率ADCP水池宽度和深度计算 # Y& c& p: ?& e1 Q8 A
根据2008年实施的《HY/T 102-2007声学多普勒流速剖面仪检测方法》[12]及国标[3]提供的ADCP工作频率分类，以美国RDI公司生产的Workhorse系列ADCP产品参数为例，可计算得出不同工作频率ADCP对试验水池的宽度和深度需求，详见表2。 0 \7 m4 p7 W3 O1 R$ \, f
表2 骏马系列ADCP产品参数及水池宽度和深度需求计算
8 W+ q& \; p. v# g* L; z+ ^ 型号/频率（kHz） $ g& `) S8 x6 j( k
Long Ranger 75
" C% O8 h$ P; j. c5 R% c Quartermaster 150
8 G; f9 t5 t6 y- e) q Sentinel 300
8 I# E: v( ~! Q( u Sentinel 600 ( X( t' y' a% ^7 M9 G* r; L
Sentinel 1200 2 i- {$ O4 u, H' n! ~5 n
起始测量深度（盲区）/m 1 l* | P$ v R# f1 ?# l& H3 J; J
8.0
5 z( z1 X' y9 L3 Z 4.0 3 i! i" Z- e$ B' c# t0 v2 T
2.0 . Z; J i/ e1 C4 Q' E
1.0
6 r1 ]- S" n+ @ 0.5 3 Z; L3 u9 Y* W" x [0 m/ e5 l
波束角θ/°
1 T6 R$ [; f5 n7 n! L 20 $ m6 \- {3 f& P4 v" {6 S3 h
20 % J3 u8 G0 _9 S% q& g+ d$ T) V1 T
20
; M8 S! N( N; b- O0 |7 o2 l. _ 20
, M, f% |8 n/ `2 Q- F Q 20
. O4 p8 |7 q3 b# s3 H/ y 波束宽2β/° & c; j/ P) r9 J5 Q+ w# g( |5 U
4 # i) ~0 v8 l8 U% R
4
1 R' d) K+ u% U9 a- m# q 4
5 o5 Z: l. \( \6 d- M; j 4 % D, {9 l8 u6 T X
4 6 v; | T9 z4 ~
第一层深度单元至换能器距离h1/m
* m8 Y6 G0 r" u 12.71
4 O3 K1 E5 R* W9 F; c 6.35
# {3 B) y" x) A( j4 Q 3.21
. ~0 t( s1 ?, e7 G! U3 Y 1.6
6 _; T3 K7 ^: X/ [ t* U2 S 0.8 0 z2 u8 h! ^/ p! s0 t2 _
最小深度单元c/m 1 W D t+ e8 R# u1 M; y, y4 R
4
2 t% u& p) F0 n; {6 C 2 - d2 F: `* Q+ @
1 ) {/ F. e2 h: n5 d& b/ b
0.5 $ ^8 e3 t3 @% m8 F* p
0.25
5 j3 a! K/ W9 T( ^) F% r 仪器直径d/m
0 h. t8 L7 j% ]# g" X. R6 v3 z) i 0.55
. S* A) _$ J: R: m 0.49
% h! N- f$ A- A% Y 0.23
8 E1 i) p6 @# A- f5 A2 k 0.23
' y# W* @% G: C, f 0.23 . P y6 M1 K/ M4 Y0 n
水池深度depth/m & A; M; A5 ]+ L" W# }2 C
22.54 ( N' D* o( @" O3 P
11.51 ; m. a+ R* w* i$ _, v. {
6.04
, p$ k! d8 A& k9 }& k 3.27
4 O) b9 S7 Q4 ? G+ z; W6 ~ 1.88 " P7 K3 b+ x0 O H/ K2 k- g6 ] e
不抑制旁瓣水池宽度with/m
& O8 P! n/ W9 e9 J 31.71
, j, A/ J, M- A; `* S6 [ 16.06 5 z) r+ D: m0 |5 f; o& U
8.06 : k9 L2 M+ W! p: [: o9 o
4.15
7 J* o$ ~' U! L4 `/ g0 q 2.18 1 G- z1 P& f" p; P0 [2 I
抑制旁瓣水池宽度with/m 3 ?1 F! F& l4 U7 ]' O
12.23 # ?0 A* N5 n+ h
5.83
/ p% I( A) E* [8 s5 v' }# b 2.93 & F- Q8 @9 Q2 f) x5 _
1.70 - V2 ]$ E6 I# p$ v6 a% G) W% V
0.95 ! K# F4 ]* l- h8 y/ e7 l; _
五、结论 . H5 ~. D# W4 W. K! I
为满足当前ADCP计量检定的迫切需求，本文依据水声学及ADCP工作原理，参考国外相关经验，结合ADCP检定现状，对ADCP水池拖车试验中水池池体需求提出了新的计算方法，推导出了具体的计算公式，得出了不同工作频率ADCP检定所需的水池尺度需求值。其中，在计算水池宽度时，分别考虑了存在旁瓣和抑制旁瓣的情况，并进行了计算比较，得出了在检定时若能抑制ADCP换能器旁瓣，可大幅降低对水池宽度的需求。本文的研究量化了水池拖车试验对池体的具体需求，进一步完善了ADCP计量检定的理论基础，为ADCP在室内可控条件下检定规范的制定及检定平台设施建设提供一定的参考与借鉴作用。 - Z0 |* O1 F, {" ~# J3 M: N
参考文献： 5 v! r" e& M3 D/ s4 ~) A
[1]李雨,袁德忠,周波,等.ADCP在水文测验中的应用及其发展前景[J].人民长江,2013(z2):35-38.
: g# s U; H3 P [2]周庆伟,张松,汪小勇,等.声学多普勒剖面流速仪检测方法探讨[J].海洋技术学报,2016,35(4):31-35. * g' m" W+ p! F
[3]中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会.GB/T 24558-2009声学多普勒流速剖面仪[S].北京:中国标准出版社,2009. # Y/ B' x" i0 F" P
[4]赵越,英小勇.声学多普勒流速剖面评价方法介绍与探讨[J].水利信息化,2010(4):58-61．
" ~" Z2 `* X E5 r- x5 {' N [5]姚永熙.水文用声学多普勒剖面流速仪的计量检定[J].水利技术监督,2008,NO.2:45-49． - U, Y$ O F; G; ^* \0 o
[6]李占桥,田志光,吴宝勤,等.两种型号RDI ADCP测速精度分析[J].海洋测绘,2012,32(6):57-59．
& R4 {; V6 k2 O$ i; d2 o [7]voulgarisG,Trowbridge JH.Evaluationofthe acousticDoppler velocimeter (ADV) for turbulence measurements[J].Journal of Atmospheric&Oceanic Technoligy,1998,15(15):272-289. - M5 y4 |/ o, s1 b3 ]
[8]Shih H H, Payton C, Sprenke J, et al. Towing Basin Speed Calibration of Acoustic Doppler Current Profiling Instruments[R]. Joint Conference on Water Resource Engineering and Water Resources Planning and Management. 2000:1-10. ) o4 k; [/ b4 L
[9]刘德铸.声学多普勒流速测量关键技术研究[D].哈尔滨：哈尔滨工程大学，2010.
' v- e) w: Q- l* b$ n% d( K* Y [10]胡波.声学多普勒流速剖面仪检测与检验关键技术研究[D].天津:南开大学,2013. 8 E5 f% T! c4 z7 ^( u/ r" W6 S
[11]智永明,邵军,唐跃平,等.ADCP流速仪检定拖车驱动方案探讨[J].中国测试,2016,(42)11:79-93. # G" C: P( L6 @2 i4 D" [. M
[12]国家海洋局.HY/T 102-2007声学多普勒流速剖面仪检测方法[S].北京:中国标准出版社,2008. + h6 s; o( r( `- f
【作者简介】第一作者乔正明，1984年出生，男，山东商河人，工程师，本科，主要从事海洋装备计量检定研究。本文为基金项目，国家自然科学基金面上项目(41775027),国家自然科学基金青年科学基金项目(41405025)；文章来自《海洋测绘》（2018年第5期），版权归《海洋测绘》所有，转载请备注论文作者，说明文章来源，并请备注由“溪流之海洋人生”微信公众平台整理。   
! ]! E2 [3 z' L! x- i 相关阅读推荐 8 B4 i/ }# @. E3 y7 x0 x3 X& E* X, H
征文作品▏阳明：ADCP测流为什么不做声速剖面改正 6 U) w, w: P5 u1 T0 E. s3 J4 _
论文专区▏海流测量技术发展及应用 + Y+ N+ J( h+ ]/ R, R
海洋军事▏海洋环境对潜艇活动的影响
$ x4 {0 V( L( @. F; ^+ }/ m3 U, D z 院士论坛▏吴立新等：物理海洋观测研究的进展与挑战
. Z8 V/ b6 Z& T4 l1 ? 海洋技术▏欧洲国家的海洋观测系统及其对我国的启示      - N; S1 j; F3 ?/ v: j
公众号 . [/ m2 s8 z4 I; t( R4 J
溪流之海洋人生 4 O/ X9 |0 D* Z- ~' ?
微信号▏xiliu92899  4 {% O2 B3 c( L0 ^: n
用专业精神创造价值
: E' V% R3 Z8 {+ Z. T8 g 用人文关怀引发共鸣 您的关注就是我们前行的动力 $ a& V4 z8 _/ b' a, y+ k% a3 B
投稿邮箱▏452218808@qq.com  
: @$ J4 D( Y* s; E% W6 j7 P2 i; j; c
# }' E( c" [4 x! z0 J
, s9 R, b' D6 U2 i
. a; d) S2 [7 m# ?& _ | 
