岸礁地形对港内低频波浪共振影响的数值研究

港湾共振是指当近岸海域的低频波浪传播到港湾内,如果其周期和港湾本征周期相一致,那么就会激起港湾内剧烈的水体共振的现象．港湾共振可以诱发港内泊船的剧烈运动,显著影响系泊安全和货物装卸[1]．实际港口的本征周期大多在几分钟至1 h范围内．在自然界中,能够诱发显著的港湾共振的外力有多种．这些外力包括低频波浪、大气压力扰动、海啸波以及近岸滑坡或结构物坠落产生的冲击波等[2]．然而,对于大部分港口,最常见的外力往往是从外海入射的波群．典型的风浪和涌浪总是以波群的形式向前传播,尽管它们的周期仅为几秒至十几秒,但由于波浪非线性,波群会产生二阶锁相长波,而这些长波往往会诱发港内长周期的水体共振现象[3]．

文献[4]中通过对澳大利亚的Two Rocks港进行实地观测,证实了Two Rocks港内的水体共振现象是由伴随着入射波群传播的锁相长波导致的．他们还通过Boussinesq模型模拟发现:当外海波群向岸传播时,相对于外海低频波浪的能量,在近岸珊瑚礁系统上的低频波浪能量大约增大8～10倍,进而显著增大了Two Rocks港内水体发生长周期共振的风险[5]．文献[6-7]中通过实验和数值模拟的方法也证实了近岸珊瑚礁地形能够显著增大低频波浪的能量．

刚度与贮藏时间的关系，如图4所示。圣女果刚度随贮藏时间延长，呈逐渐下降趋势，表明圣女果抵抗变形的能力在逐渐下降。贮藏前6天下降较快，6天之后下降速度开始逐渐减缓。

文献[8]中提出了基于最小二乘法的分离港内低频波浪的方法,并基于此波浪分析技术,研究了在港湾发生第1共振模态条件下短波波长的变化对于港内锁相长波和自由长波波幅以及它们的相对成分的影响．随后,文献[9]中将此研究拓展到最低的4个共振模态,并进一步研究了入射短波波幅的变化对于它们的影响．然而,文献[8-9]中均将港外地形假设为平底,并未考虑近岸地形对于港湾长周期共振的影响．与此不同，文中不再假设港外地形为平底,而是将近岸珊瑚礁地形考虑在内．

文中将系统研究近岸珊瑚礁地形外坡坡度和形状变化对于港湾共振条件下港内自由长波、锁相长波以及它们相对成分的影响．受到文献[10]的启发，文中将使用斜坡地形、反正切型地形和双曲余弦型地形等不同形状地形作为近岸珊瑚礁外坡地形．文中的研究将对我国南海岛礁上的港口建设提供理论指导,以期能够降低这些港口发生港湾共振的风险．

1 数值模型和分析技术

1.1 数值模型

Boussinesq方程是近些年来被广泛地应用于模拟波浪在近岸的传播过程的一种数值模型．它可以精确地模拟波浪的色散性、非线性、浅水变形、破碎和波成流等一系列现象[11]．文中采用FUNWAVE模型[12]来进行数值模型实验．该模型以文献[13]中提出的色散性精确到二阶的完全非线性Boussineq方程为基础开发．

采用笛卡尔坐标系,x和y轴位于静水面, z轴以垂直向上为正,则控制方程表达式为:

 

(1)

式中:

(·uαt)+zα[·(huαt)]-

(·u)+

[·(hu)]

(2)

M= (h+η)u+(h+η)·

(3)

V2= [(zα-η)(uα·)(·(huα))+

[(·(huα)+η·uα)2]

采用SPSS 16.0软件进行统计学分析。计数资料分析采用χ2检验或Fisher 精确概率法检验；生存分析及生存曲线绘制采用Kaplan Meier法和Log-rank 检验；患者预后多因素相关分析采用COX比例风险回归模型。以P＜0.05为差异有统计学意义。

(4)

（68）圆叶鞭苔 Bazzania conophylla（Sande Lac.）Schiffn. 熊源新等（2006）；杨志平（2006）

文献[16]中采用FUNWAVE模型模拟了文献[17]中的物理实验,通过将模拟得到的港湾共振条件下1～3倍频波浪成分沿港内中轴线的分布与物理实验结果进行比较,发现它们之间吻合很好,表明该数值模型能够精确地模拟港湾共振现象．

(1)岩棉复合型保温模板的组成材料为岩棉保温板、冷拔低碳钢丝、聚合物保温砂浆，该模板具有免拆模、使用安全、耐火等级高、施工便捷等优点。

1.2 分离方法

文中采用文献[8]中提出的港内低频波浪的分离程序来将港内的低频波浪成分进一步分解为低频锁相驻波和低频自由驻波成分．

图1呈现了港口平面和测点布置情况．港口长度L=50.0 m,宽度W=5.0 m．港池内的水深常数hI=1.0 m．沿着港口中轴线等间距地布置N=21个测点．相邻两个测点的间距D=2.5 m．测点G01布置在港口底墙处．海岸线和港内边墙均设置为全反射直墙．笛卡尔坐标系统(o,x, y, z)的原点取在静水面,z轴以垂直向上为正．

  

图1 实验布置和测点位置示意图Fig.1 Definition sketch of the experimental setupand the location of the wave gauges

在港外,由频率分别为f1和f2的自由短波组成的双色波群垂直于口门入射．根据文献[3]中对双色波垂直口门入射情况下狭长矩形港内波浪的公式描述,港内的低频波浪成分可以表示为:

ηL= ζbcos[(2πΔf)t-αb]cos(Δkx)+

ζfcos[(2πΔf)t-αf]cos(κx)

初识老胡，是在一次图书采访研讨会期间。会上，他沉默寡语，一如他平日做派，敏事慎言，以“行”代“言”，稍微点头便是算对发言者的首肯了。会下，他同我交流地方志采访事宜，我便聊起上海图书馆方志收藏数量、特点及开发现状。未及说完，他悠悠一句：何日有暇，请来敝馆看看我们收藏的戏剧方志，那可是上图不及的啊！虽令我语塞，但尔后几番打听，得知老胡任职的上海戏剧学院图书馆在其独到眼力和勤勉采集之下，专题类方志文献量大质高，特色独具，并且关于戏剧类专业资料独霸沪上，傲视群雄，便不由得对老胡有格外起敬之意。

(5)

此外，新一轮农地改革政策正在加紧制定。土地承包期再延长30年政策制定工作已启动，不久后将出台。还将深化农村集体产权制度改革，全面开展农村集体资产清产核资、集体成员身份确认，推进集体经营性资产股份合作制改革。资源变资产、资金变股金、农民变股东“三变”改革红利将加快释放。

Δf=|f1-f2|

式中:

Δk=|k1-k2|

(2πf1)2=gk1tanh(k1hI)

(2πf2)2=gk2tanh(k2hI)

表1列出了文中数值实验所用的珊瑚礁外坡的几何形状参数．表中表示珊瑚礁外坡段的平均水深,可用下式计算:

式(5)等号右端第1项为低频锁相驻波,第2项为低频自由驻波．ζb、ζf分别为低频锁相驻波和低频自由驻波的波幅．Δk、κ分别为低频锁相驻波和低频自由驻波的波数．αb、αf分别为低频锁相驻波和低频自由驻波的初始相位．k1和k2分别为港内短波f1和f2成分的波数．

对式(5)进行一系列的三角变换,可以得到港内低频波浪波幅沿着中轴线变化的解析表达式:

A(x)2=aφ0(x)+bφ1(x)+cφ2(x)

(6)

式中:

 

(7)

通过对各测点测得的波面时间序列进行离散序列的傅里叶级数展开可计算出差频成分在各测点处的波幅Ai,其中i=1, 2,…, N．显然φ0(x)、φ1(x)和φ2(x)是线性无关的,可以构成1组基底,则参数a、b和c可以通过运用最小二乘法拟合式(6)和实测值来得到,然后通过式(7)即可以求出ζb、ζf和|αb-αf|．

文献[8-9]中使用已知的解析信号对该低频波浪分离程序的精确性进行检验,证明了它能够在不同的共振模态条件下精确地将港内的低频波浪成分进一步分解为锁相长波和自由长波．

2 数值实验布置

使用针对狭长型矩形港口发生共振的线性解析解[18],可以计算出文中所用港口的放大因子曲线(图2)．

  

图2 港口的放大因子曲线Fig.2 Amplification factor curve for the rectangle harbor

可见,图1所示港口的最低5个模态的共振频率分别为0.014、0.043、0.073、0.103和0.133 Hz．文中仅研究双色波群垂直于港口口门入射并诱发港内产生第1共振模态的情况．文中所有数值实验中,入射双色波群中的短波频率f1和f2均分别设置为0.250、0.264 Hz,则伴随波群传播的锁相长波的频率Δf=0.014 Hz,其与港口的第1模态的共振频率相一致,因此在双色波群的作用下港内将发生第1模态的水体共振[3]．对于入射双色短波波幅,则考虑3种情况:① a1=a2=0.02 m；② a1=a2=0.05 m；③ a1=0.05 m、a2=0.02 m．a1、a2分别为入射f1 和f2 短波成分的波幅．前两种情况下,入射双色波群为完全调制的,而对于第3种情况,波群为弱调制的．

图3为文中所用的数值波浪水池和港口内外地形布置示意图．数值波浪水池尺寸为1 257.8 m×2 084.2 m．港外的珊瑚礁礁坪的长度为50.0 m,礁坪处的水深与港内水深相同,均为hI=1.0 m．为了研究珊瑚礁外坡坡度对于港内低频波浪的影响,外坡设置为常坡度,在100.0 m≤x≤(100.0+b)m范围内,水深可解析表达为:

h(x)=hI+γ(x-100) (常坡度地形)

(8)

外坡宽度b在5.0～100.0 m之间变化．外海的水深常数h0=5.0 m,则外坡坡度γ在0.04～0.80的范围内变化．根据线性色散关系,频率Δf=0.014 Hz的低频自由长波在水深h0=5.0 m时的波长为499.9 m,为了充分吸收向外海辐射的低频长波,在数值波浪水池的右边界和上下边界分别布置了宽度分别为1 002.8 m和1 000.6 m的海绵层,均略大于低频自由长波波长的2倍．因为图1所示港口沿中轴线对称,为了节省计算网格和时间,仅模拟港口和外海的一半区域．港内等间距布置了N=21个测点．x方向网格尺寸Δx从港内至造波区均为0.5 m,而在海绵层内则从0.5 m逐渐增大到18.7 m;y方向网格尺寸Δy从港内的0.5 m逐渐增大至19.9 m．计算域中包含285 200个网格数和286 431个节点数．所有的实验中时间步长均取为0.025 s．数值模型模拟到2 000.0 s,总计80 000个时间步．

由图1可知，卤汁1在复卤过程中总酸含量总体上比卤汁2要高，卤汁1复卤过程中总酸最低含量为0.91 g/kg，最高含量为3.88 g/kg；卤汁2复卤过程中总酸最低含量为1.13 g/kg，最高含量为3.79 g/kg。卤汁中总酸变化的原因可能是，在豆干卤制过程中，工厂每天要进行补料和补水，所以总酸会在一定的范围内波动。卤汁1在卤制过程中，总酸含量在2.71 g/kg上下波动；卤汁2在卤制过程中，总酸含量在2.31 g/kg上下波动。而变化不大的原因是工厂每天会对卤汁总酸进行监控。所以总酸只在一定的范围内波动[9]。

uα为zα=-0.531h时x和y方向的水质点速度,η为自由水面．对以上的控制方程采用文献[14]中提出的有限差分格式进行求解,对方程中空间的一阶和二阶导数分别采用五点和三点中心差分格式,边界附近的点采用相应的偏心差分格式．控制方程的时间离散则使用预估—校正法进行处理．预估步骤采用三阶显式Adams-Bashforth格式,校正步骤采用四阶隐式Adams-Moulton格式．为方便边界的处理,采用域内源造波与海绵层结合的方法．域内源造波采用文献[15]中提出的方法,即在计算域内设置造波区,产生向造波区一侧传播的波浪．海绵层即在固壁边界前设置1～2倍波长的人工阻尼区使波浪衰减．

  

图3 数值波浪水池示意图Fig.3 Sketch of the numerical wave basin

  

图4 当b=50.0 m时珊瑚礁外坡不同形状地形比较Fig.4 Comparison of various bottom shapes over the reef slope under the condition of b=50.0 m

为了研究珊瑚礁外坡形状对港内低频波浪的影响,在珊瑚礁外坡宽度b=50.0 m和100.0 m时,文中还考虑了反正切型地形和双曲余弦型地形,即在100.0 m≤x≤(100.0+b)m范围内,水深可解析表达为:

图8为当珊瑚礁外坡宽度b=50.0 m时,图4所示的不同珊瑚礁外坡形状对港内锁相长波和低频自由长波波幅的影响．可以发现,对于文中考虑的3种不同的入射短波波幅,港内低频自由长波波幅均随着珊瑚礁外坡平均水深的增大而显著增大,而港内锁相长波波幅则随着的增大而逐渐减小．以a1=a2=0.05 m的情况为例(图8(b)),当从1.94 m增大到4.42 m时,港内自由长波的相对波幅ζf/a从0.83逐渐增大到0.91,而港内锁相长波的相对波幅ζb/a则从0.39减小到0.32．图9呈现了不同珊瑚礁外坡形状对振幅比ζb/ζf的影响．可以发现,振幅比ζb/ζf随着珊瑚礁外坡平均水深的增大而逐渐减小,并且随着入射短波波幅的增大而略有减小，同图7．

(9)

h(x)=hIcoshκ[μ(x-100)] (双曲余弦型地形)

现有文献中，Sarkar曾提出企业环境行为是企业为了平衡环境和经济效益而推行的一系列战略措施,其目的可能是响应来自外界的压力,或是主动减轻企业经营生产对社会环境造成的危害。[2]在此之后,学者们又将企业环境行为的内涵逐渐细化至具体方向及措施。学者王宇露指出企业为响应环境问题而采取的绿色行为包括处置废弃物、降低能源消耗、采购生态友好型产品等。[3]结合国内外学者的观点，企业环境行为就是企业在其生产经营中对环境造成影响的行为，同时为挽救这些负面影响所采取的一系列措施以及为平衡企业经济利益与社会环境效益而施行的决策。

(10)

几何参数α、β、ν和μ用于决定珊瑚礁外坡的形状,它们满足如下关系:

在饲养中牦牛易患出血性败血症，严重时极易导致牦牛死亡。而牦牛饲养作为高原地区的主要经济来源，一旦出现疫情，就会对当地的经济产生重大影响。因此，在牦牛饲养中，一定要注意病害的防治工作，从思想上进行重视，保证畜牧业的健康发展。

 

(11)

 

(12)

(2πΔf)2=gκtanh(κhI)

 

(13)

为了直观显示不同珊瑚礁外坡地形的差别,图4为b=50.0 m时不同珊瑚礁外坡形状的比较．对于反正切型地形,因为α取值不同,使得珊瑚礁外坡具有不同的平均水深．α=2.7 m和3.3 m对应外坡平均水深4.42 m和3.70 m．对于双曲余弦型地形,不同的κ值使得珊瑚礁外坡具有不同的平均水深．ν=0.003和50对应的外坡平均水深分别为2.48 m和1.94 m．常坡度外坡地形的平均水深为3.0 m,介于α=3.3 m的反正切型地形和ν=0.003的双曲余弦型地形之间．当b=100.0 m时,各种珊瑚礁外坡地形和图4中相应地形相似,仅珊瑚礁的外坡宽度变为原来的2倍．

除了服饰，教堂里还展出了拉祜族、苗族、彝族等5种少数民族文字的《圣经》版本。这些民族或支系中的一部分群众，因为历史原因较早地接受了基督教信仰。近年来，云南省基督教“两会”与联合圣经公会合作，将《圣经》翻译为多种少数民族文字版本，对保护少数民族语言和民族文化遗产做出了贡献。这是一项巨大的文化工程，北辰教堂也在积极参与，给予支持。这样的翻译工作，对于少数民族文字、文化的传承有积极的意义。

 

表1 珊瑚礁外坡形状参数Table 1 Geometric parameters of the reef slope inall the numerical experiments

  

外坡形状类型形状参数常坡度地形b=5.0～100.0m,γ=0.04～0.80,h=3.0m反正切型地形b=50.0mα=2.7m,h=4.42mα=3.3m,h=3.70mb=100.0mα=2.7m,h=4.42mα=3.3m,h=3.70m双曲余弦型地形b=50.0mν=0.003,h=2.48mν=50,h=1.94mb=100.0mν=0.003,h=2.48mν=50,h=1.94m

3 结果和分析

3.1 珊瑚礁外坡坡度对于港内低频波浪的影响

当波群传播到港内时,水体将从初始的平静状态逐渐达到稳定状态．文中对于港内初始反应阶段的波况不予考虑,仅对达到稳定状态的港内波况进行研究[8-9]．对各测点达到平稳的波面时间序列采用离散序列傅里叶级数展开的方法提取出差频成分的波幅,然后使用文中的波浪分离程序将差频成分进一步分解为锁相驻波和低频自由驻波成分．图5为当b=50.0 m时在珊瑚礁外坡形状分别为α=2.7 m的反正切型、常坡度型和ν=50的双曲余弦型条件下，由不同波幅双色波群入射产生的港内各测点实测的差频成分波幅和采用最小二乘法拟合得到波幅包络线．图中AL表示差频成分波幅,a=(a1a2)1/2为入射双色短波的平均波幅．可以看出,对于不同波幅的入射双色波和不同形状的珊瑚礁外坡,港内实测的数据点与拟合得到的曲线均吻合得很好,这表明使用文中的波浪分离方法计算得到港内低频自由长波和锁相长波的结果是精确可靠的．文中以a1=a2=0.02 m时的情况为例．当珊瑚礁外坡形状为α=2.7 m的反正切型时(图5(a)),分离得到的锁相驻波波幅ζb=4.61×10-3m,低频自由驻波波幅ζf=9.99×10-3m,相位差|αb-αf|=97.16°;当珊瑚礁外坡形状为常底坡型时(图5(b)),ζb=4.64×10-3m, ζf=9.12×10-3m,|αb-αf|=86.70°; 当珊瑚礁外坡形状为ν=50的双曲余弦型时(图5(c)), ζb=4.75×10-3m, ζf=9.02×10-3m,|αb-αf|=74.09°．文中用振幅比ζb/ζf来定量表示锁相驻波和低频自由驻波的相对成分．在以上3组实验中,振幅比分别为46.12%、50.75%和52.72%．

  

图5 当b=50.0 m时港内各测点实测的差频成分波幅(点)和使用最小二乘法拟合得到的波幅包络线Fig.5 Subharmonic component amplitudes (dot) in all gauges and their envelope line obtained by usingthe separation procedure under the condition of b=50.0 m

图6为当珊瑚礁外坡为常坡度地形时港内低频自由驻波波幅ζf和锁相驻波波幅ζb随着珊瑚礁外坡坡度系数(即坡度的倒数)的变化．由图可知,对于文中所考虑的3种入射波波幅,珊瑚礁外坡坡度的变化对于港内低频自由驻波成分的影响较大．具体来说,港内低频自由驻波的波幅随着珊瑚礁外坡坡度的增大而增大,而港内锁相驻波成分则对珊瑚礁外坡坡度的变化并不敏感．在文中所研究的外坡坡度变化范围内,港内锁相驻波波幅仅随着外坡坡度的增大略有减小．同时,由于低频自由长波在港内被共振放大,其波幅始终显著大于港内锁相长波的波幅,这与文献[8-9]中的发现相一致．图7为在不同入射波波幅条件下振幅比ζb/ζf随着外海珊瑚礁外坡坡度系数的变化．可以发现,无论入射双色短波波幅的大小以及双色波群是否为完全调制的,振幅比ζb/ζf均随着外坡坡度系数的减小(即珊瑚礁外坡坡度的增大)而减小．另外,还可以发现随着入射短波波幅的增大,振幅比ζb/ζf略有减小．对于入射波波幅a1=a2=0.02 m的情况,ζb/ζf在44.39%～55.29%;当a1=0.05 m、a2=0.02 m时,ζb/ζf则在38.25%～53.17%;当a1=a2=0.05 m时,ζb/ζf则在33.37%～48.52%．

  

图6 港内低频自由驻波和锁相驻波波幅随着珊瑚礁外坡坡度系数的变化Fig.6 Variations of amplitudes of bound and free longstanding waves with the slope coefficient of thereef slope

3.2 珊瑚礁外坡形状对于港内低频波浪的影响

h(x)=hI+αatan[β(x-100)] (反正切型地形)

  

图7 振幅比ζb/ζf随着珊瑚礁外坡坡度系数的变化Fig.7 Variations of the amplitude ratios, ζb/ζf, with the slope coefficient of the reef slope

  

图8 当b=50.0 m时,港内低频自由驻波和锁相驻波波幅随着珊瑚礁外坡平均水深的变化Fig.8 Variations of amplitudes of bound and free longstanding waves with the mean water depth overthe reef slope, under the condition of b=50.0 m

  

图9 当b=50.0 m时,振幅比ζb/ζf随着珊瑚礁外坡平均水深的变化Fig.9 Variations of the amplitude ratios, ζb/ζf,withthe mean water depth over the reef under the condition of b=50.0 m

图10为当珊瑚礁外坡宽度b=100.0 m时,珊瑚礁外坡形状对港内锁相长波和低频自由长波波幅的影响．此时,珊瑚礁外坡形状对港内自由长波的影响则与b=50.0 m时有显著的不同．港内自由长波波幅不再随着的增大而单调增加,而是出现了先减小后增大的变化规律．港内锁相长波随的变化规律则与b=50.0 m时相同．图11呈现了相应的振幅比ζb/ζf随着的变化．与图9所示不同,当b=100.0 m时,振幅比ζb/ζf随着的增大出现了先增大后减小的变化特征．并且,当 m时,越大的入射短波波幅反而得到越大的振幅比,该现象并未在图7、9中出现过．

  

图10 当b=100.0 m时,港内低频自由驻波和锁相驻波波幅随着珊瑚礁外坡平均水深的变化Fig.10 Variations of amplitudes of bound and free longstanding waves with the mean water depth overthe reef under the condition of b= 100.0 m

  

图11 当b=100.0 m时,振幅比ζb/ζf随着珊瑚礁外坡平均水深的变化Fig.11 Variations of the amplitude ratios, ζb/ζf, with the mean water depth over the reef under the condition of b=100.0 m

4 结论

使用FUNWAVE模型模拟了由双色波群诱发的狭长矩形港内的二阶长波共振现象．文中考虑了近岸珊瑚礁地形对于港湾共振的影响．采用一个基于最小二乘法的分离港内低频波浪的方法,对港口发生第一共振模态时，港内锁相长波和自由长波成分及其它们的相对成分随着珊瑚礁外坡坡度和形状变化进行了系统研究,并得到以下结论:

数控维修专业中的相关知识和内容都有极强的实践性与操作性，在科学技术发展之下，将数字化技术使用在教学中，以此建设数字化课程。经过对课程内容的数字化、实训环境的数字化以及网络课程的建设，将数控维修专业数字化课程建设落到实处，以此提升学生的综合能力。

(1) 对于文中所研究的特定港口和珊瑚礁外坡坡度的变化范围,港内自由长波的波幅随着珊瑚礁外坡坡度的增大而增大,而港内锁相驻波成分则对珊瑚礁外坡坡度的变化并不敏感．前者由于在港内被共振放大,其波幅始终要显著大于后者的波幅．振幅比ζb/ζf随着珊瑚礁外坡坡度的增大而减小．

(2) 珊瑚礁外坡形状对港内自由长波和锁相长波的影响有所不同．对于港内自由长波,珊瑚礁外坡形状对它的影响还与珊瑚礁外坡宽度有关．当b=50.0 m时,港内自由长波波幅随着珊瑚礁外坡平均水深的增大而增大;而当b=100.0 m时,港内自由长波波幅随着珊瑚礁外坡平均水深的增大出现了先减小后增大的变化趋势．对于港内锁相长波,当b=50.0、100.0 m时,它均随着珊瑚礁外坡平均水深的增大而逐渐减小．

以上结论仅针对文中所使用的特定港口、珊瑚礁外坡坡度和形状变化范围．

参考文献(References)

[ 1 ] LPEZ M, IGLESIAS G. Long wave effects on a vessel at berth[J]. Applied Ocean Research, 2014, 47: 63-72. DOI:10.1016/j.apor.2014.03.008.

[ 2 ] RABINOVICH A B.Seiches and harbor oscillations[M].Handbook of Coastal and Ocean Engineering, 2009: 193-236.DOI:10.1142/9789812819307-0009.

[ 3 ] BOWERS E C. Harbour resonance due to set-down beneath wave groups[J]. Journal of Fluid Mechanics,1977,79：71-92.DOI:10.1017/s0022112077000044.

[ 4 ] THOTAGAMUWAGE D T, PATTIARATCHI C B. Observations of infragravity period oscillations in a small marina[J]. Ocean Engineering, 2014, 88: 435-445. DOI:10.1016/j.oceaneng.2014.07.003.

[ 5 ] THOTAGAMUWAGE D T, PATTIARATCHI C B. Influence of offshore topography on infragravity period oscillations in Two Rocks Marina, Western Australia[J]. Coastal Engineering, 2014, 91: 220-230. DOI:10.1016/j.coastaleng.2014.05.011.

[ 6 ] NWOGU O, DEMIRBILEK Z. Infragravity wave motions and runup over shallow fringing reefs[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2010, 136(6): 295-305. DOI:10.1061/(asce)ww.1943-5460.0000050.

[ 7 ] YAO Y, HUANG Z, MONISMITH S G, et al. 1DH Boussinesq modeling of wave transformation over fringing reefs[J]. Ocean Engineering, 2012, 47: 30-42. DOI:10.1016/j.oceaneng.2012.03.010.

[ 8 ] DONG G, GAO J, MA X, et al. Numerical study of low-frequency waves during harbor resonance[J]. Ocean Engineering, 2013, 68: 38-46. DOI:10.1016/j.oceaneng.2013.04.020.

[ 9 ] GAO J, JI C, GAIDAI O, et al. Numerical study of infragravity waves amplification during harbor resonance[J]. Ocean Engineering, 2016, 116: 90-100. DOI:10.1016/j.oceaneng.2016.02.032.

中国汽车维修行业协会原会长康文仲，交通运输部科学院原党委书记、中国汽车保修设备行业协会原副会长王晓曼，中国汽车维修行业协会原会长孙守仁，中国汽车保修设备行业协会原副会长罗大柱 ，中国汽车保修设备行业协会原副会长张熙伦，中国汽车保修设备行业协会原副会长兼秘书长田国华，中国汽车维修行业协会常务副秘书长王逢铃， 中国汽车保修设备行业协会原常务副秘书长刘蕴，中国汽车保修设备行业协会秘书长刘建农，中国汽车维修行业协会连锁工委秘书长严雪月，广东省道协机动车维修检测分会会长罗少泽，北京汽车维修行业协会副秘书长侯金凤及数十位国内知名汽保企业负责人作为嘉宾出席此次茶话。

[10] GAO J, JI C, LIU Y,et al. Numerical study on transient harbor oscillations induced by solitrary waves[J]. Ocean Engineering,2016,126:467-480.

[11] KIRBY J T. Boussinesq models and applications to nearshore wave propagation, surf zone processes and wave-induced currents[J]. Elsevier Oceanography Series,2003,67(3):1-41.

[12] KIRBY J T, LONG W, SHI F. Funwave 2.0 fully nonlinear boussinesq wave model on curvilinear coordinates [R]. [S.l.]:University of Delaware, 2003.

[13] WEI G, KIRBY J T, GRILLI S T, et al. A fully nonlinear Boussinesq model for surface waves. I. highly nonlinear, unsteady waves[J]. Journal of Fluid Mechanics, 1995,294(294):71-92.

[14] WEI G, KIRBY J T. Time-dependent numerical code for extended boussinesq equations[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1995, 121(5): 251-261. DOI:10.1061/(asce)0733-950x(1995)121:5(251).

[15] CHAWLA A, KIRBY J T. A source function method for generation of waves on currents in Boussinesq models[J]. Applied Ocean Research, 2000, 22(2): 75-83. DOI:10.1016/s0141-1187(00)00005-5.

[16] DONG G, WANG G, MA X, et al. Numerical study of transient nonlinear harbor resonance[J]. Science China Technological Sciences, 2010, 53(2): 558-565. DOI:10.1007/s11431-009-0409-5.

然后进行超声引导下盆腔肿块穿刺活检，未进行超声引导下盆腔肿块穿刺活检为排除病例。术前对盆腔肿块行常规超声评估，观察盆腔肿块与周围脏器及血管的关系，盆腔肿块内部血供分布。患者取平卧位或膀胱结截位，可选经腹部或经阴道穿刺。根据穿刺目标大小选取合适的切割槽，在超声实时引导下进针穿刺活检，见图2。通常采用单点多角度或多方向活检，穿刺活检次数2-3次，送检。术后观察病人情况，记录并发症发生情况，做好交接、病例资料随访等工作。

[17] ROGERS S R, MEI C C. Nonlinear resonant excitation of a long and narrow bay[J]. Journal of Fluid Mechanics, 1978, 88(1): 161. DOI:10.1017/s0022112078002037.

[18] MEI C C. The applied dynamics of ocean surface waves[M]. New York:Wiley,1983.