Spatiotemporal distribution pattern and transport mechanism of suspended sediments in Longkou offshore under the action of tidal current
-
摘要: 根据2015年1月与2017年5月两期大潮期27小时海流连续同步观测及悬浮泥沙资料，利用Morlet 小波分析与单宽悬沙通量机制分解法，结合潮流场特征，分析了龙口湾海域悬浮泥沙的时空分布特征，研究了潮流作用下的悬沙输运机制。结果表明，平面上，悬沙浓度在龙口湾内较低，湾口、湾外海域的的悬沙浓度相对较高；垂向上，研究区整体表现为悬沙浓度由底层向表层递减。潮周期尺度上，悬浮泥沙浓度峰值较流速峰值存在0.5～2 h的时间滞后；悬浮泥沙浓度在时间上的变化以12～16 h尺度为主要周期。季节尺度上，1月悬沙浓度远远大于5月。研究区内单宽输沙通量为1.21～239.77 g·s−1·m−1，整体上悬沙通量呈现出湾口最高、湾外开阔海域其次、湾内最低的分布格局。平流输运项（T1+T2）在悬沙输移中占绝对优势，潮致余流在悬沙输运起主要作用。季节上，1月份单宽输沙通量远远高于5月份，整体呈现出夏季落淤储沙而冬季再悬浮输沙的季节性输运特征。湾外海域悬浮泥沙大致向SW、S向输运，屺坶岛以北海域沿W、SW向屺坶岛头输运，进入湾内后，悬浮泥沙的输运沿人工岛北侧和西侧分为两支。研究区内，在海流、地形等的综合作用下，整体表现为由渤海向莱州湾的净输沙趋势。Abstract: Based on the data collected from tidal current observation and suspended sediment sampling at the coastal Longkou Bay, which had lasted for a total of 27 hours in January, and May 2017, the characteristics of tidal current field and the transport mechanism of suspended sediments were studied with the method of Morlet Wavelet Transform and flux mechanism analysis. The results show that the concentration of suspended sediments is lower within the bay, but relatively higher at the mouth and outside areas of the bay. Vertically, the concentration of suspended matters gradually decreases from the bottom to top. The peak of suspended sediment concentration has a half to two hour time lag behind the peak of current velocity, and the variation in suspended sediment concentration shows a cycle of 12～16 h variation. The net sediment transport rate ranges from 1.21 to 239.77 g·s−1·m−1. In general, the suspended sediment flux outside the bay is higher than that inside. In terms of seasonal variation, the suspended matters are far higher in January than that in May, indicating the seasonal transport characteristics: storage in summer and transport in winter. Suspended sediments move towards SW and S in offshore waters, while in the water to the north of the Qimu Island follow a W and SW path to the head of the island. When it comes into the bay, the transport of suspended sediment is separated into two branches along the north and west sides of the artificial island respectively. Under the joint effect of ocean current, topography and so on, the net sediment transport is mainly from the Bohai Sea to Laizhou Bay in the study area.
-
随着海洋地震勘探的发展，勘探目标逐渐趋于复杂化，高精度成像方法已经成为海洋油气勘探的瓶颈技术。各向异性现象普遍存在于地下介质中，基于各向同性的传统偏移方法难以满足实际生产中高精度成像的需求，会严重影响最终偏移结果及后续的偏移速度分析精度[1]。由于海洋拖缆地震勘探观测系统的特殊性，为了确保偏移精度，对偏移方法需提出更高要求。因此，有必要发展一种适应于海洋观测系统的各向异性高精度偏移成像方法，进一步提高地震资料解释的精度，解决海洋复杂油气藏精细勘探的问题。
近年来在三维实际资料处理中，以射线理论为基础的高斯束偏移方法取得了较好的应用效果，已成为各大油田公司关注的焦点，同时也是众多地球物理学家研究的热点之一。一方面，高斯束方法提供了地震波在地下传播的物理解释，更易在复杂的地球物理模型中描述波动现象，且由于射线追踪角度无限制，适合陡倾角成像问题；另一方面，与波动方程偏移相比，高斯束方法较为高效，非常适合三维实际资料成像处理。目前，工业界常用的Kirchhoff偏移也具有很高的成像效率，但其用于复杂探区成像处理时，射线焦散现象[2-4]及多路径问题[5-8]难以避免，从而影响成像精度。而高斯束偏移采用复值初始束参数动力学射线追踪, 能避免焦散点附近振幅奇异问题，且不同方向的局部平面波分解使其能够独立地从不同方向传播到地下成像点，不会产生多路径问题。总之，高斯束偏移不仅具有接近波动方程偏移的精度，而且具有Kirchhoff偏移灵活高效的特点。
高斯束方法起源于电磁学，在20世纪80年代初期，Červený和Popov将其引入到地球物理学领域，并取得了飞速发展[9-13]。Hill[14]基于高斯束表征的格林函数推导了叠后深度偏移公式，奠定了高斯束偏移的理论框架。其成像思想为：基于局部倾斜叠加近似，将分解后的平面波分量沿着不同方向以半速度反传到地下成像点处进行成像。后来，考虑到地下介质的复杂性，各向异性介质中的高斯束偏移方法的相关研究不断涌现，Alkhalifah[15]最早通过引入各向异性射线追踪方程发展了各向异性高斯束叠后偏移方法。Zhu等[16]通过引入相速度和群速度简化的各向异性射线追踪体系，提出了共炮域各向异性高斯束叠前深度偏移方法，进一步降低了基于弹性参数算法的计算量。段鹏飞等[17]在局部角度域提出了一种相对经济的动力学射线追踪近似方案。Protasov[18]提出了各向异性介质中针对多分量数据的高斯束偏移方法。张凯等[19]利用各向异性介质运动学和动力学射线追踪体系，提出了角度域各向异性高斯束逆时偏移方法。李振春等[20]基于转换波射线追踪算法实现了各向异性介质角度域转换波高斯束偏移方法。随着勘探精度要求逐步提高，常规的高斯束偏移难以满足高精度成像的要求，因此许多新的方法及优化算法相继出现。一方面，基于波动方程偏移中逆时偏移波场延拓思想的高斯束逆时偏移[19,21-23]极大地提高了高斯束偏移在复杂构造中的成像精度，有利于地下复杂构造中层位及断层识别等解释工作；另一方面，基于振幅保真理论[24-26]及最小二乘理论[27-29]的高斯束偏移在保幅成像及提高成像分辨率方面取得重大突破，极大地提高了地震解释过程中的AVA分析及AVO反演精度。此外，黄建平等[30]通过修改高斯束局部倾斜叠加提出了双复杂条件下非倾斜叠加精确束偏移方法。张瑞等[31]结合高斯束传播过程中反射张角信息的计算，发展了高斯束高精度叠前深度偏移方法，并在复杂海洋模型中取得了较好的应用效果。Yang和Zhu[32]在传统高斯束偏移基础上基于菲涅尔带思想提出了一种优化的高斯束偏移方法，提高了稀疏数据的偏移成像精度。
为了弥补高斯束偏移在处理低信噪比资料上的不足，出现了一种改进的高斯束偏移方法—控制束偏移方法。Vinje等[33]提出了各向同性介质的控制束偏移方法，并将其应用到裂缝基底成像、盐丘底边界成像、低覆盖次数及宽方位角数据成像等方面。Zhou等[34]将各向同性介质控制束偏移推广到各向异性介质中，发展了一种各向异性高保真度控制束偏移方法。Casasanta等[35]实现了基于共偏移距的PS转换波高斯束偏移。黄建平等[36]认为可以在τ-p域内设定阈值，针对不同信噪比的实际数据进行适当的切除来实现控制束偏移。
本文在传统高斯束偏移的基础上，考虑到地下介质的各向异性和实际资料中含有较强噪声干扰，根据有效信号和干扰信号在τ-p域中的相干性差异，基于各向异性运动学和动力学射线追踪，发展了一种VTI介质共偏移距域数据驱动控制束偏移理论方法。在高斯束偏移成像过程中，计算τ-p道集的相似系数，通过设定相似系数阈值对干扰信号进行控制，从而降低偏移剖面中的随机噪声，提高同相轴的连续性。通过各向异性洼陷模型、修改的SEG/Hess VTI模型对新方法进行测试，并与传统高斯束偏移结果进行对比验证本文方法的正确性及有效性；最后通过对中国东部某海上探区的实际资料进行处理，探讨了本文方法的优势及其对海上数据的适应性。
1. 方法原理
1.1 各向异性射线追踪
各向异性高斯束偏移方法的核心在于利用各向异性运动学和动力学射线追踪求解高斯束的复值走时和振幅。Červený[37]由各向异性介质无源弹性波波动方程经过一系列数学推导，最后得到各向异性运动学射线追踪方程组
$$\left\{ {\begin{aligned} &{\frac{{{\rm{d}}{x_i}}}{{{\rm{d}}\tau }} = {a_{ijkl}}{p_l}{g_j}{g_k}}\\ &{\frac{{{\rm{d}}{p_i}}}{{{\rm{d}}\tau }} = - \frac{1}{2}\frac{{\partial {a_{njkl}}}}{{\partial {x_i}}}{p_n}{p_l}{g_j}{g_k}} \end{aligned}} \right.$$ (1) 其中，aijkl=cijkl/ρ表示密度归一化的弹性参数；τ表示沿中心射线路径的旅行时；pi=∂τ/∂xi为慢度矢量；gi为极化矢量。
在二维情况下，式（1）中特征向量可以表示为
$$\left\{ {\begin{aligned} & {{g_1}{g_1} = \frac{{{\Gamma _{33}} - G}}{{{\Gamma _{11}} + {\Gamma _{33}} - 2G}}} \\ &{{g_3}{g_3} = \frac{{{\Gamma _{11}} - G}}{{{\Gamma _{11}} + {\Gamma _{33}} - 2G}}} \\ & {{g_1}{g_3} = \frac{{ - {\Gamma _{13}}}}{{{\Gamma _{11}} + {\Gamma _{33}} - 2G}}} \end{aligned}} \right.$$ (2) 式中，Γik=aijklpjpl为Christoffel矩阵，G为程函方程，当其取不同特征值时代表不同的程函方程。
根据式（1）即可求得中心射线路径从而构建射线中心坐标系，在各向异性介质中射线中心坐标系不再正交，因此需要在求解过程中引入一个沿中心射线路径的权值变量以校正这种非正交性所带来的误差。基于各向同性介质动力学射线追踪方程组，Hanyga[38]推导出了各向异性介质中的动力学射线追踪方程组
$$\left\{ {\begin{aligned} &{\frac{{{\rm{d}}Q}}{{{\rm{d}}\tau }} = MP + VQ}\\ &{\frac{{{\rm{d}}P}}{{{\rm{d}}\tau }} = - VP - HQ} \end{aligned}} \right.$$ (3) 其中，P和Q是动力学射线追踪参量，M、V、H表示程函方程对n和pn的导数，满足下列关系式
$$\left\{ {\begin{aligned} &{H = \frac{1}{2}\frac{{{\partial ^2}{G_m}}}{{\partial p_n^2}} - \frac{1}{4}{{\left( {\frac{{\partial {G_m}}}{{\partial {p_n}}}} \right)}^2}} \\ & {M = \frac{1}{2}\frac{{{\partial ^2}{G_m}}}{{\partial {n^2}}} - \frac{1}{4}{{\left( {\frac{{\partial {G_m}}}{{\partial n}}} \right)}^2}} \\ & {V = \frac{1}{2}\frac{{{\partial ^2}{G_m}}}{{\partial {p_n}\partial n}} - \frac{1}{4}\frac{{\partial {G_m}}}{{\partial {p_n}}}\frac{{\partial {G_m}}}{{\partial n}}} \end{aligned}} \right.$$ (4) 上式中的Gm=aijklpiplgjgk为Christoffel方程中行列式det（Γjk−Gmδjk）=0的特征值，代表各向异性介质中存在的三种不同地震波的程函方程，当m取值1、2和3时分别表示qP波、qSV波和qSH波。式（1）和式（3）中的各向异性介质运动学和动力学射线追踪方程均以弹性参数表示，目前人们习惯使用物理意义更易理解的Thomsen参数来表征弱各向异性介质，同时可以通过Thomsen参数与弹性参数转化关系来消除这种不便。
1.2 数据驱动控制束偏移
叠前高斯束偏移是通过高斯束构建炮点和检波点在地下成像点处的格林函数，然后通过它们来构建上行波场和下行波场，最后利用上、下行波场的互相关得到偏移结果。本文从频率域的互相关成像条件出发来推导出数据驱动控制束偏移公式，根据Zhang等[39]的文章可知，二维频率域互相关成像条件可以表示为
$$I({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}}) = - i\int {{\rm{d}}\omega {\rm{sgn}}(\omega ){P_{{\rm{up}}}}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega )P_{{\rm{down}}}^*} ({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega )$$ (5) 式中，Pup和Pdown分别表示上、下行波场，*表示复共轭。Gray和Bleistein[24]给出了其表达式如下：
$$\begin{aligned} {P_{{\rm{up}}}}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega ) =& 2\int {{\rm{d}}{x_{\rm{r}}}{P_{{\rm{up}}}}({\bf{x}_{\bf{g}}},{\bf{x}_{\bf{s}}},\omega )\frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}}\times\\&{\sqrt { - i\omega } {A_{{\rm{ART,g}}}}{\rm{exp}}[ - i\omega {\tau _{\rm{g}}}]} \\=& - 2\int {{\rm{d}}{x_{\rm{g}}}{P_{{\rm{up}}}}({\bf{x}_{\bf{g}}},{\bf{x}_{\bf{s}}},\omega )\frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}}\times\\&i\omega {G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{g}}},\omega ) \end{aligned}$$ (6) $$\begin{aligned} P_{{\rm{down}}}^*({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega ) = & - 2\frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}\sqrt { - i\omega } {A_{{\rm{ART,s}}}}{\rm{exp}}[ - i\omega {\tau _{\rm{s}}}]\\ = &2\frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}i\omega {G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega ) \end{aligned}$$ (7) 式中，Vs和Vg分别表示炮点xs和检波点xg处的速度，θs和θg分别表示炮点和检波点处的出射角，τs、τg和AART,s、AART,g分别表示渐进射线理论下的实值走时和振幅，G（x0,x,ω）是波动方程中的格林函数，将式（6）和式（7）带入到式（5）中可得如下共炮域偏移成像公式：
$$\begin{aligned} I({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}}) = & - 4\int {{\rm{d}}\omega i\omega \left| \omega \right|\frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}{G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega )}\times \\ & \int {{\rm{d}}{x_{\rm{g}}}{P_{{\rm{up}}}}({\bf{x}_{\bf{g}}},{\bf{x}_{\bf{s}}},\omega )\frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}{G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{g}}},\omega )} \end{aligned}$$ (8) 在高斯束方法中，格林函数是通过一系列由震源点x出射的且具有不同出射方向的高斯束的叠加积分表示的：
$$G({\bf{x}_{\bf{0}}},{\bf{x}},\omega ) = \frac{i}{{4\pi }}\int {\frac{{{\rm{d}}{p_x}}}{{{p_z}}}} A({\bf{x}_{\bf{0}}},{\bf{x}}){\rm{exp}}[i\omega T({\bf{x}_{\bf{0}}},\bf{x})]$$ (9) 其中，px和pz分别为射线参数的水平和垂直分量，T（x0,x）和A（x0,x）分别是复值走时和振幅，它们可以通过式（1）和式（3）表示的各向异性射线追踪方程得到。对于海上接收的共偏移距地震记录，可以通过如下变换：
$${\bf{m}} = \frac{1}{2}({\bf{x}_{\bf{g}}} + {\bf{x}_{\bf{s}}}),{\bf{h}} = \frac{1}{2}({\bf{x}_{\bf{g}}} - {\bf{x}_{\bf{s}}})$$ (10) 将式（8）中的积分变量由炮点-检波点域变换到中心点-半偏移距域，
$$\begin{aligned}I({\bf{x}_{\bf{0}}},{\bf{h}}) =& - 4\int {{\rm{d}}\omega i\omega \left| \omega \right|\frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}{G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega )}\times\\&\int {{\rm{d}}{m_x}P({\bf{h}},{\bf{m}},\omega )\frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}{G^*}({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{g}}},\omega )} \end{aligned}$$ (11) 其中，P（h,m,ω）表示中心点坐标为m，半偏移距为h的共偏移距道集，这种类型的道集可以理解为半偏移距h固定，沿中心点坐标m变化的地震记录。接着，需要沿中心点坐标来确定中心点位置Lm，在确定Lm后，结合式（9）表示的格林函数，根据如下关系：
$${\bf{x}_{\bf{s}}} = {\bf{m}} - {\bf{h}},{\bf{x}_{\bf{g}}} = {\bf{m}} + {\bf{h}}$$ (12) 将中心点Lm有效范围内的炮点、检波点格林函数分别用Lm−h、Lm+h处出射的高斯束来近似表示，
$$\begin{aligned} G({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{s}}},\omega ) =& \frac{i}{{4\pi }}\int {\frac{{{\rm{d}}{p_{sx}}}}{{{p_{sz}}}}} A({\bf{x}_{\bf{0}}},{\bf{L}_{\bf{m}}} - {\bf{h}})\times\\&{\rm{exp}}[i\omega T({\bf{x}_{\bf{0}}},{\bf{L}_{\bf{m}}} - {\bf{h}})]\\ G({\bf{x}_{\bf{0}}},{\bf{x}_{\bf{g}}},\omega ) =& \frac{i}{{4\pi }}\int {\frac{{{\rm{d}}{p_{gx}}}}{{{p_{gz}}}}} A({\bf{x}_{\bf{0}}},{\bf{L}_{\bf{m}}} + {\bf{h}})\times\\&{\rm{exp}}[i\omega T({\bf{x}_{\bf{0}}},{\bf{L}_{\bf{m}}} + {\bf{h}})] \end{aligned}$$ (13) 其中，psx、psz和pgx、pgz分别为炮点和检波点射线参数的水平和垂直分量。由于高斯束具有一定的宽度，所以检波点处会接收到不同波束的能量贡献。为了区分这些波束，可以结合高斯束在初始位置波前面为平面的特点，采用局部倾斜叠加方法来将这些波束的能量反传到地下成像点。根据Hill[14]文章中给出的一维高斯窗函数：
$$\frac{{\Delta {x_L}}}{{\sqrt {2\pi } {w_0}}}\sqrt {\left| {\frac{\omega }{{{\omega _r}}}} \right|} \sum\limits_{{x_L}} {{\rm{exp}}\left[ { - \left| {\frac{\omega }{{{\omega _r}}}} \right|\frac{{{{\left( {{x_g} - {x_L}} \right)}^2}}}{{2w_0^2}}} \right]} \approx 1$$ (14) 式中，xL为高斯窗中心位置，w0是高斯窗宽度，ωr为参考频率，将式（13）和式（14）代入到公式（11）中可得
$$\begin{aligned} &I({\bf{x}_{\bf{0}}},{\bf{h}}) = \frac{{\Delta L}}{{4{\pi ^2}\sqrt {2\pi } {w_0}}}\sum\limits_{{x_L}} {\int {{\rm{d}}\omega i\omega \frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}\int {\frac{{{\rm{d}}{p_{sx}}}}{{{p_{sz}}}}} A_s^*} }\times\\ & {\rm{exp}}\left[ { - i\omega T_{\rm{s}}^*} \right]\times \frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}\int {\frac{{{\rm{d}}{p_{gx}}}}{{{p_{gz}}}}} A_{\rm{g}}^*{\rm{exp}}\left[ { - i\omega T_{\rm{g}}^*} \right]{D_{\rm{S}}}({x_L},{p_m},\omega ) \end{aligned}$$ (15) 其中，pm=ps+pg表示中心点射线参数，DS（xL,pm,ω）表示共偏移距道集加窗局部倾斜叠加，
$$\begin{aligned}{D_{\rm{S}}}({x_L},{p_m},\omega ) =& {\left| {\frac{\omega }{{{\omega _r}}}} \right|^{3/2}}\int {{\rm{d}}{m_x}} P({\bf{h}};{\bf{m}};\omega )\times\\&{\rm{exp}}\left[ { - i\omega {p_m}({m_x} - {x_L}) - \left| {\frac{\omega }{{{\omega _r}}}} \right|\frac{{{{({m_x} - {x_L})}^2}}}{{2w_0^2}}} \right]\end{aligned}$$ (16) 式（15）表示的是先对频率、炮点出射的高斯束以及检波点处出射的高斯束进行三重积分，然后再对不同中心点所对应的三重积分求和的过程。考虑到高斯束中的振幅和走时均为复数，可以令
$$\begin{array}{l} {T_{{\rm{re}}}} = {\rm{Re}}(T_{\rm{s}}^* + T_{\rm{g}}^*),{A_{{\rm{re}}}} = {\rm{Re}}(A_{\rm{s}}^*A_{\rm{g}}^*)\\ {T_{{\rm{im}}}} = {\rm{Im}}(T_{\rm{s}}^* + T_{\rm{g}}^*),{A_{{\rm{im}}}} = {\rm{Im}}(A_{\rm{s}}^*A_{\rm{g}}^*) \end{array}$$ (17) 直接求解公式（15）需要面对一系列频率求和的巨大计算量，根据高成等[40]，可以通过改变积分次序，并且应用傅里叶反变换和希尔伯特变换，推导出式（15）对应的时间域表达式
$$\begin{aligned} & I({\bf{x}_{\bf{0}}},{\bf{h}}) =\frac{{\Delta L{\omega _r}}}{{4\pi \sqrt {2\pi } {w_0}}}\sum\limits_{{x_L}} {\int {\frac{{{\rm{d}}{p_{sx}}}}{{{p_{sz}}}}} \int {\frac{{{\rm{d}}{p_{gx}}}}{{{p_{gz}}}}} } \frac{{\cos {\theta _{\rm{s}}}}}{{{V_{\rm{s}}}}}\frac{{\cos {\theta _{\rm{g}}}}}{{{V_{\rm{g}}}}}\times\\ &[{A_{{\rm{re}}}}{{\bar D}_{\rm{S}}}({x_L},{p_m},{T_{{\rm{re}}}},{T_{{\rm{im}}}}) - {A_{{\rm{im}}}}{{\bar D}_{{\rm{S}},{\rm{H}}}}({x_L},{p_m},{T_{{\rm{re}}}},{T_{{\rm{im}}}})] \end{aligned}$$ (18) 根据Hale[41]的研究可知，成像过程中可首先评估虚值走时Tim的范围，然后对Tre和Tim进行离散采样，并计算
${\bar D_{\rm{S}}}$ 和${\bar D_{{\rm{S,H}}}}$ 的离散样板表，最后根据地下成像点的旅行时信息与离散样板表进行匹配成像。该过程可以理解为一个数据驱动的过程，即从道集中提取每一样点采集“属性特征”和旅行时，然后找到这些特定数据点所指向的成像点。对含有噪声干扰的海上共偏移距道集而言，噪声会在一定程度上影响样点信息的正确提取，从而影响最终的偏移结果。通过对含噪声的共偏移距道集进行加窗局部倾斜叠加可以发现，τ-p坐标图版上的峰值显示在局部共偏移距道集中存在相干同相轴，并且这个同相轴来自地下的某一特定反射段，而一些相对较弱的能量团可能是一些背景和随机噪声等干扰信号。为了降低随机噪声对偏移结果产生的影响，本文采用相干体技术来提高偏移剖面的信噪比。相干体技术不仅能用于检测地震波同相轴的不连续性，也可以用来识别断层、特殊岩性体以及河道等，在地震资料解释中扮演着重要角色[42]。目前，相干体算法主要可以划分为三类：基于相关的算法、基于相似系数的算法以及基于本征结构分析的算法[43-45]。本文采用基于相似系数的算法，推导τ-p道集的相似系数计算公式。根据Neidell等[46]可知，相似系数可以定义为归一化的输出/输入能量比值，其中输出道能量表示为输入道的一个简单合成或者求和。在信号分析中相似系数的表达式
$${S_{\rm{C}}} = \frac{1}{N}\sum\limits_{i = 1}^N {{{\left(\sum\limits_{j = 1}^J {{f_{i,j}}} \right)}^2}} \bigg/\sum\limits_{i = 1}^N {\sum\limits_{j = 1}^J {f_{i,j}^2} } $$ (19) 式中，N表示在一个短时窗内的样点数，J表示总道数，下标i和j分别表示样点序号和道号，fij为样点（i, j）处的振幅值。在此基础上，根据式（19）可以推导出τ-p道集的相似系数计算公式（20），同时根据有效信号和干扰信号在τ-p域中的相干性差异，可以应用相似系数阈值滤波公式（21）对离散样板表进行滤波处理，去除τ-p道集中的部分干扰能量，从而提高偏移剖面中同相轴的连续性。
$${S_{\rm{C}}}(\tau ,{x_L},{p_m}) = \frac{{\displaystyle\sum\nolimits_{i = 1}^N {\left[ {\displaystyle\sum\nolimits_{j = 1}^J {u(\tau + i\Delta t - {p_m}\Delta x,{x_j})} } \right]} }}{{N\displaystyle\sum\nolimits_{i = 1}^N {\displaystyle\sum\nolimits_{j = 1}^J {{u^2}(\tau + i\Delta t - {p_m}\Delta x,{x_j})} } }}$$ (20) $$W({\tau _i},{x_L},{p_j}) = \left\{ {\begin{array}{*{20}{c}} {1,{S_{\rm{C}}}{\text{≥}} {\rm{tp}}}\\ {0,\;{\rm{others}}\;} \end{array}} \right.$$ (21) 其中，u表示局部共偏移距道集，xj为高斯窗内第j道的笛卡尔坐标，xL为参考道的位置；τ表示走时，Δt和Δx分别表示时间采样间隔和空间采样间隔；tp为相似系数阈值参数，范围为0～1。
如图1所示，VTI介质共偏移距域数据驱动控制束偏移过程可以概括为如下4个步骤：（1）计算τ-p道集对应的相似系数并选取大于相似系数阈值参数对应的τ-p道集中的能量；（2）利用各向异性射线追踪分别计算来自炮点和检波点的中心射线以及旁轴射线；（3）组合两组旁轴射线来计算两个“胖射线”相交区域上所有点的走时和振幅；（4）根据旅行时及射线参数关系选取τ-p道集中的能量进行匹配得到偏移后的子波束，所有偏移后的子波束叠加得到最终的成像结果。单偏移距数据驱动控制束具体的实现流程如图2所示。
![]()
图 1 共偏移距域数据驱动控制束偏移原理图a. 偏移后的子波束，b. τ-p道集，c. 相似系数。Figure 1. Diagram of the data-driven CBM in common-offset domaina. migrated sub-beam，b. τ-p gather，c. semblance coefficient.![]()
图 2 单偏移距数据驱动控制束偏移程序流程图Figure 2. Flow chart of the single-offset data-driven controlled beam migration program2. 应用实例
从三个方面来进一步测试本文方法的正确性、有效性及实用性。首先，在编程实现算法的基础上，对简单的各向异性洼陷模型进行偏移处理，并与常规高斯束偏移结果对比，验证本文方法的正确性及其在处理存在各向异性地层中的优势；其次，对修改的SEG/Hess VIT复杂模型进行测试，验证本文方法的有效性及处理低信噪比数据的优势。最后，通过对中国东部某海上探区实际资料进行成像处理，验证本文方法对海上数据的适应性。
2.1 各向异性洼陷模型
为了验证本文方法的正确性，建立了如图3a所示的各向异性洼陷模型。正演模拟采用交错网格各向异性高阶有限差分正演方法：震源为主频20 Hz的Ricker子波；接收方式采用共偏移距接收方式，偏移距范围从−2 000 m到2 000 m，偏移距间隔为20 m，共201个共偏移距道集，每个共偏移距道集有401道，道间距为20 m，总记录时间为3 s，采样率为4 ms。图3b和3c分别显示了一个近偏移距（200 m）和远偏移距（1 500 m）道集。基于相同的成像参数，分别采用各向同性和各向异性共偏移距域高斯束偏移方法对其进行偏移成像，结果如图4所示。从图4a和4b对比可以发现，对于近偏移距数据，采用各向同性方法时，由于忽略了各向异性的影响，弯曲底界面成像位置不准确，曲率发生了明显的变化（黑色箭头所示）。对比图4c和4d的远偏移距结果可以看出，各向异性参数对共偏移距道集中的大偏移距信息成像影响较大，各向同性结果中弯曲底界面处存在明显的能量发散问题。图4e和4f为所有偏移距叠加后成像结果，采用各向异性偏移方法可以很好地聚焦成像（白色箭头所示），界面同相轴清晰且位置准确，界面弯曲曲率与实际模型更加吻合。
![]()
图 3 各向异性洼陷模型a. 各向异性参数场，b. 200 m偏移距道集，c. 1 500 m偏移距道集。Figure 3. Anisotropic sag modela. anisotropic parameters，b. common-offset gather （200 m），c. common-offset gather （1 500 m）.![]()
图 4 各向异性洼陷模型偏移结果a. 各向同性（偏移距为200 m），b. 各向异性（偏移距为200 m），c. 各向同性（偏移距为1 500 m），d. 各向异性（偏移距为1 500 m），e. 各向同性所有偏移距叠加，f. 各向异性所有偏移距叠加。Figure 4. Results of GBM on anisotropic sag modela. isotropic （offset=200 m），b. anisotropic （offset=200 m），c. isotropic （offset=1 500 m），d. anisotropic （offset=1 500 m）， e. all offset stack isotropic result，f. all offset stack anisotropic result.2.2 修改的SEG/Hess VTI模型
为了进一步验证本文方法的有效性，接着对不包含盐丘构造的部分SEG/Hess VTI模型进行测试，其各向异性参数如图5所示，可以看出，该模型具有典型的正断层、岩性低速体及各向异性透镜体，具有一定的构造代表性。共偏移距道集通过交错网格各向异性高阶有限差分正演模拟得到，震源为主频20 Hz的Ricker子波，为了得到完整的共偏移距道集，采取模型扩边处理，总共201个共偏移距道集，偏移距范围从−2 000 m到2 000 m，偏移距间隔为20 m，每个共偏移距道集有451道，道间距为20 m，总记录时间为3 s，采样率为4 ms。
基于相同的成像参数，分别采用各向同性和各向异性共偏移距域高斯束偏移对该模型进行处理，结果如图6a和6b所示。通过对比可以看出，各向异性相比各向同性而言，在断层附近成像位置更加准确且断点绕射能量收敛（图6红色箭头），地层的成像能量更加聚焦且同相轴更加清晰（图6蓝色箭头）。为了更加清楚地对比，将图6中绿色矩形框内成像结果进行放大显示（图6c和图6d）。从放大图中可以看出，由于各向异性高斯束偏移考虑地下各向异性的影响，可以大幅改善断点的成像效果（图6c和图6d蓝色箭头），并且图6c中透镜体未能很好成像，而各向异性的成像效果（图6d红色箭头）远好于各向同性偏移结果。
![]()
图 6 修改的SEG/Hess VIT模型共偏移距域高斯束偏移结果及局部放大显示对比a. 各向同性，b. 各向异性，c. 各向同性局部放大，d. 各向异性局部放大。Figure 6. Results of common-offset GBM and magnified views on modified SEG/Hess VTI modela. isotropic，b. anisotropic，c. isotropic magnified view，d. anisotropic magnified view.为了测试共偏移距域控制束偏移的成像效果，选取了偏移距为200 m的共偏移距道集（图7a）作为测试对象，并使用Seismic Unix软件包中的suaddnoise命令在原始道集中加入随机噪声得到低信噪比共偏移距道集（图7b）,加噪后的共偏移距道集中的同相轴变得模糊并混叠在背景噪声中。图8为图7中局部共偏移距道集（蓝色矩形框）相干性分析结果，通过对低信噪比局部共偏移距道集（图8a）进行倾斜叠加得到图8b所示的τ-p道集，与图8d所示的原始道集生成的标准τ-p道集对比可以发现，随机噪音在τ-p道集中表现为一些比较弱的能量，同时可以注意到即使是在标准τ-p道集能量也不能很好聚焦（图8d红色箭头），这主要是因为局部τ-p变换本身存在噪声泄露问题，然而这些能量却可以在控制束中得到压制（图8e、8f红色箭头）。随后通过计算图8b对应的相似系数，结果如图8c所示，可以发现不同时刻、不同中心点射线参数对应大小不同的相似系数，而不同中心点射线参数又分别对应着不同出射方向的一组炮检点高斯束，对于其中某一时刻，与最大相似系数对应的中心点射线参数控制着高斯束在炮检点主要的出射方向，即优势出射方向。而一些较小相似系数可能对应一些随机噪音的出射方向。因此，如果想要在偏移结果中压制由随机噪音引起的干扰，需要计算得到高斯束的优势出射方向，即找到最大相似系数对应的τ-p道集。为了解决上面的问题，我们通过设定相似系数阈值对τ-p道集进行滤波，图8e和图8f是分别采用阈值参数tp=0.4和0.6进行滤波得到的结果，对比可以发现，当取合适的阈值参数时随机噪声在τ-p道集中的能量可以得到一定压制（蓝色椭圆）；但是，过大的阈值参数也可能压制了部分有效能量（红色椭圆），使得有效能量得不到成像。
![]()
图 7 交错网格各向异性高阶有限差分正演模拟的共偏移距（200 m）道集a. 无噪共偏移距道集，b. 加噪共偏移距道集。Figure 7. Common-offset gathers （200 m） obtained by staggered grids anisotropic finite difference forwarda. no-noise common-offset gathers，b. common-offset gathers with noise.![]()
图 8 图7中局部共偏移距道集相干性分析a. 图7b蓝框中局部共偏移距道集，b. τ-p道集，c. 相似系数，d. 图7a蓝框中局部共偏移距道集生成的标准τ-p道集，e. 相似系数阈值tp=0.4滤波后的τ-p道集，f. 相似系数阈值tp=0.6滤波后的τ-p道集。Figure 8. Semblance analysis of local common-offset gathers in Fig.7a. local common-offset gathers marked by blue box in Fig. 7b，b. τ-p gather，c. semblance coefficient，d. standard τ-p gather，e. filtering by tp=0.4，f. filtering by tp=0.6.图9为偏移距为200 m共偏移距域各向异性高斯束及控制束偏移结果，整体上看，常规高斯束偏移结果（图9b），可以发现存在明显的偏移噪声，控制束偏移结果（图9c和9d）中噪声得到明显压制，成像剖面的信噪比显著提高，这对于复杂地区构造解释是非常有用和不可或缺的。同时可以注意到，tp=0.6偏移剖面（图9d）与tp=0.4偏移剖面（图9c）相比而言，阈值较大时信噪比更高，特别是浅层同相轴更加清晰（蓝色矩形框），更接近于图9a原始道集偏移结果，但同时也会损失部分有效地层能量，如模型深部蓝色箭头所示。
![]()
2.3 海上数据处理结果
最后测试新方法对海洋共偏移距实际资料的适应性。图10a、10b和10c分别为中国东部某海上探区的二维各向异性模型的VP0、ε和δ参数场，该地区包含一个大断层及洼陷构造。此各向异性模型参数场是通过传统射线层析方法构建的，海上接收的共偏移距道集已经过带通滤波等预处理，偏移距范围从150 m到4 650 m，偏移距间隔为100 m，共46个共偏移距道集，记录总时间为3.6 s，采样间隔为2 ms，其中偏移距为150 m的共偏移距道集如图10d所示。由图10d可知，该共偏移距道集中含有较为明显的背景及随机噪声，资料信噪比较低，适合进行控制束偏移处理。
![]()
图 10 中国东部某海上探区实际数据a. VP0，b. ε，c. δ，d. 150 m偏移距道集。Figure 10. The field data of marine survey in Eastern Chinaa. VP0，b. ε，c. δ，d. common-offset gather （150 m）.基于相同的成像参数，分别采用各向同性和各向异性共偏移距域高斯束偏移对该实际资料进行处理，结果如图11a和11b所示。为了更加清楚地对比，将图11中白色矩形框内成像结果进行放大显示（图11c和图11d）。从放大图中可以看出，由于考虑地下各向异性参数（更接近于真实地下介质情况）对偏移结果的影响，可以大幅改善高角度构造的成像质量，其结果更能反映地下实际地质构造，陡倾断面更加清晰、连续。为了测试不同相似系数阈值参数对偏移结果的影响，分别采用阈值参数tp=0、0.4、0.6和0.8进行测试得到如图12a-d所示的偏移结果。同样通过对其蓝色矩形框内的偏移结果局部放大对比（图12e-h）可以看出，当阈值参数较小时，成像噪音较严重，当取合适的阈值参数时可大幅提高成像结果的信噪比。但是，过大的阈值参数也可能压制了部分有效能量，使得有效能量得不到成像。所以针对不同的数据资料，需要通过局部共偏移距道集相干性分析来选取合适的阈值参数，使噪音得到压制，有效能量得到保持，最终得到一个较为理想的高精度偏移剖面。
![]()
图 11 海上数据共偏移距域高斯束偏移结果及局部放大对比a. 各向同性，b. 各向异性，c. 各向同性局部放大，d. 各向异性局部放大。Figure 11. Results of GBM and magnified views of the marine surveya. isotropic，b. anisotropic，c. isotropic magnified view，d. anisotropic magnified view.![]()
图 12 海上数据共偏移距域控制束偏移结果及局部放大对比a. tp=0，b. tp=0.4，c. tp=0.6，d. tp=0.8，e. tp=0局部放大，f. tp=0.4局部放大，g. tp=0.6局部放大，h. tp=0.8局部放大。Figure 12. Results of CBM and magnified views of the marine surveya. tp=0，b. tp=0.4，c. tp=0.6，d. tp=0.8，e. tp=0 magnified view，f. tp=0.4 magnified view，g. tp=0.6 magnified view，h. tp=0.8 magnified view.3. 结论
本文根据有效信号和干扰信号在τ-p域中的相干性差异，基于各向异性运动学和动力学射线追踪，发展了一种VTI介质共偏移距域数据驱动控制束偏移理论方法。通过各向异性洼陷模型、修改的SEG/Hess VTI模型及海上实际资料成像处理，得到如下几点认识：
（1）各向异性参数对共偏移距道集中的大偏移距信息成像影响较大。
（2）当地层各向异性不能忽略时，新方法能够更加准确地恢复地下的复杂构造。
（3）新方法能够在一定程度上提高低信噪比数据的偏移成像效果。同样也应该注意到，在实际资料处理中，相似系数阈值参数的选取也十分关键。阈值较小时偏移噪音较强，但过大阈值也可能压制了部分有效能量或者产生偏移假象，选取合适的阈值参数才能够得到较为理想的控制束偏移剖面。因此，下一步拟针对相似系数阈值参数的选取开展深入研究。
-
![]()
图 2 各站位垂向平均潮流矢量图
Figure 2. The vector graph of depth-averaged tidal current at each station
![]()
图 3 各站位悬浮泥沙平面分布图
Figure 3. The average suspended sediment concentration of each layers
![]()
图 4 各站位流速和悬浮泥沙浓度随时间变化分布图
Figure 4. The time-based variation in current velocity and suspended sediment concentration
![]()
![]()
图 7 悬沙输运项与垂向平均余流叠置图
（注：站位A1、A2、A3、A6、B4的悬沙输运项分别缩小至原来的1/10、1/5、1/2、1/5，站位A4、A5、A6扩大2.5倍。）
Figure 7. The vector graph of suspended sediment transport item and vertical average residual flow
表 1 各站位表中底层涨-落潮平均流速（cm/s）
Table 1 The average velocity of high and low tide at each station （cm/s）
站位 A1 A2 A3 A4 A5 A6 B1 B2 B3 B4 B5 B6 涨潮 表层 44.6 25.8 24.5 14.7 9.1 8.1 31.7 44.6 30.3 43.6 25.9 23.4 中层 42.6 23.8 21.9 7.1 8.8 7.9 30.3 28.4 25.5 38.7 20.4 21.4 底层 36.3 16.7 16.7 7.6 7.0 6.7 21.5 19.1 17.5 28.8 13.0 13.3 落潮 表层 32.8 14.0 19.4 8.7 7.2 6.1 27.7 39.6 32.8 41.9 25.7 32.9 中层 31.4 22.4 20.5 11.8 9.5 7.1 24.6 28.4 23.5 39.5 20.0 19.6 底层 28.2 16.6 16.3 9.0 8.4 5.7 19.7 18.7 17.8 26.0 11.4 8.9 
下载: 导出CSV
表 2 各层位平均悬浮泥沙浓度（mg/L）
Table 2 The average suspended sediment concentration of each layers（mg/L）
层位 A1 A2 A3 A4 A5 A6 表层 53.3 65.8 49.8 23.6 22.9 16.4 中层 62.2 76.2 65.8 24.6 33.1 16.5 底层 79.7 90.0 76.9 41.2 31.9 19.9 层位 B1 B2 B3 B4 B5 B6 表层 13.6 15.4 19.7 15.6 10.7 15.9 中层 16.0 18.6 20.2 18.1 14.0 22.8 底层 19.3 20.3 22.8 20.1 17.4 19.2
下载: 导出CSV
表 3 各站位最大流速、最大悬浮泥沙浓度
Table 3 The maximum suspended sediment concentration and velocity of each station
站位 A1 A2 A3 A4 A5 A6 B1 B2 B3 B4 B5 B6 流速/（cm/s） 91.1 126.1 46.5 23.0 19.0 24.0 55.33 65.33 55.67 121.0 47.0 69.67 悬沙浓度/（mg/L） 148.0 160.2 141.8 88.8 79.4 77.8 29.0 31.8 40.0 37.8 28.0 46.0
下载: 导出CSV
表 4 各站位实测平均悬浮泥沙浓度与平均流速
Table 4 The average suspended sediment concentration and velocity of each station
站位 A1 A2 A3 A4 A5 A6 B1 B2 B3 B4 B5 B6 流速/（cm/s） 36.3 21.1 19.51 7.77 8.40 6.30 26.4 24.9 24.57 36.70 20.32 20.15 悬沙/（mg/L） 63.6 74.5 62.5 29.5 29.3 19.4 16.3 18.1 20.92 17.94 14.1 19.27
下载: 导出CSV
表 5 1月各站悬沙输运项及单宽净输沙率
Table 5 Flux decomposition and sediment net transport rate in January
g·s−1·m−1 站位 T1 T2 T1+T2（%） T3+T4 T5（%） T6+T7+T8 T A1 230.77 9.83 240.22（100.2） 0.65 3.77（1.6） 0.15 239.77 A2 99.05 4.67 103.09（104.6） 1.27 6.12（6.2） 0.18 98.60 A3 19.49 10.91 24.40（103.6） 1.76 1.78（7.6） 0.69 23.55 A4 1.65 0.54 1.16（95.9） 0.17 0.29（7.6） 0.20 1.21 A5 3.81 0.10 3.86（97.6） 0.03 0.17（4.3） 0.04 3.96 A6 2.56 1.53 1.59（81.5） 0.59 0.04（2.2） 0.03 1.95 注：T1-T8为具有方向的矢量值，因此，部分悬沙输运项在净输沙率（T）中占的比例大于100%。
下载: 导出CSV
表 6 5月各站悬沙输运项及单宽净输沙率
Table 6 Flux decomposition and sediment net transport rate in May
g·s−1·m−1 站位 T1 T2 T1+T2/% T3+T4 T5/% T6+T7+T8 T B1 20.84 0.47 21.19（102.6） 0.08 1.63（7.9） 0.02 20.66 B2 18.81 0.96 18.53（112.5） 0.07 2.06（12.5） 0.07 16.47 B3 12.63 0.63 13.18（100.0） 0.15 0.24（1.8） 0.02 13.17 B4 64.72 1.98 64.86（105.0） 1.79 1.97（2.4） 0.14 61.78 B5 5.60 0.68 5.73（110.0） 0.11 0.85（16.4） 0.03 5.21 B6 5.68 0.72 5.33（127.6） 0.05 1.95（46.7） 0.05 4.17 注：T1-T8为具有方向的矢量值，因此，部分悬沙输运项在净输沙率（T）中占的比例大于100%。
下载: 导出CSV
-
[1] 陈勇, 韩震, 杨丽君, 等. 长江口水体表层悬浮泥沙时空分布对环境演变的响应[J]. 海洋学报, 2012, 34(1):145-152. [CHEN Yong, HAN Zhen, YANG Lijun, et al. The spatial and temporal distribution of suspended sediment concentration and its response to environmental change in the Changjiang Estuary [J]. Acta Oceanologica Sinica, 2012, 34(1): 145-152. [2] 秦蕴珊, 赵一阳, 赵松龄, 等. 渤海地质[M]. 北京: 科学出版社, 1985: 1-232. QIN Yunshan, ZHAO Yiyang, ZHAO Songling, et al. Geology of the Bohai Sea[M]. Beijing: Science Press, 1985: 1-232.
[3] Webster T, Lemckert C. Sediment resuspension within a microtidal estuary/embayment and the implication to channel management [J]. Journal of Coastal Research, 2002, 36(sp1): 753-759.
[4] 乔璐璐. 冬季大风事件下渤黄海环流及泥沙输运过程研究[D]. 中国海洋大学博士学位论文, 2008. QIAO Lulu. Circulation and sediments transport due winter storms in the Bohai Sea and Yellow Sea[D]. Doctor Dissertation of Ocean University of China, 2008.
[5] 李九发. 长江河口南汇潮滩泥沙输移规律探讨[J]. 海洋学报, 1990, 12(1):75-82. [李九发. Study on sediment transport law of Nanhui tidal flat in Yangtze estuary [J]. Acta Oceanologica Sinica, 1990, 12(1): 75-82. [6] Hu R J, Ma F, Wu J Z, et al. Sediment transport in the nearshore area of Phoenix Island [J]. Journal of Ocean University of China, 2016, 15(5): 767-782. doi: 10.1007/s11802-016-2967-z
[7] XIe D F, Gao S, Wang Z B, et al. Numerical modeling of tidal currents, sediment transport and morphological evolution in Hangzhou Bay, China [J]. International Journal of Sediment Research, 2013, 28(3): 316-328. doi: 10.1016/S1001-6279(13)60042-6
[8] Liu J H, Yang S L, Zhu Q, et al. Controls on suspended sediment concentration profiles in the shallow and turbid Yangtze Estuary [J]. Continental Shelf Research, 2014, 90: 96-108. doi: 10.1016/j.csr.2014.01.021
[9] Stephens J D, Allison M A, Di Leonardo D R, et al. Sand dynamics in the Mekong River channel and export to the coastal ocean [J]. Continental Shelf Research, 2017, 147: 38-50. doi: 10.1016/j.csr.2017.08.004
[10] Yang Y P, Zhang M J, Li Y T, et al. The variations of suspended sediment concentration in Yangtze River Estuary [J]. Journal of Hydrodynamics, 2015, 27(6): 845-856. doi: 10.1016/S1001-6058(15)60547-9
[11] 秦蕴珊, 李凡. 渤海海水中悬浮体的研究[J]. 海洋学报, 1982, 4(2):191-200. [QIN Yunshan, LI Fan. Study on the suspended matter of the sea water of the Bohai gulf [J]. Acta Oceanologica Sinica, 1982, 4(2): 191-200. [12] 刘潇, 冯秀丽, 刘杰, 等. 山东半岛靖海湾及其附近海域悬沙粒度特征及再悬浮作用研究[J]. 泥沙研究, 2013(4):68-73. [LIU Xiao, FENG Xiuli, LIU Jie, et al. Characteristics of grain size distribution and resuspension process of suspended sediment in Jinghai Bay and its adjacent waters, Shandong Peninsula [J]. Journal of Sediment Research, 2013(4): 68-73. doi: 10.3969/j.issn.0468-155X.2013.04.011 [13] Yang Z S, Ji Y J, Bi N S, et al. Sediment transport off the Huanghe (Yellow River) delta and in the adjacent Bohai Sea in winter and seasonal comparison [J]. Estuarine Coastal & Shelf Science, 2011, 93(3): 173-181.
[14] 刘潇, 冯秀丽, 刘杰, 等. 山东半岛靖海湾及其附近海域悬浮泥沙分布与变化[J]. 海洋地质与第四纪地质, 2012, 32(6):9-15. [LIU Xiao, FENG Xiuli, LIU Jie, et al. Spatial and temporal distribution pattern of suspended sediment in the Jinghai Bay and adjacent waters, Shandong peninsula [J]. Marine Geology & Quaternary Geology, 2012, 32(6): 9-15. [15] 陈斌, 刘健, 高飞. 莱州湾悬沙输运机制研究[J]. 水科学进展, 2015, 26(6):857-866. [CHEN Bin, LIU Jian, GAO Fei. Suspended sediment transport mechanism in Laizhou Bay [J]. Advances in Water Science, 2015, 26(6): 857-866. [16] 庞重光, 李坤, 于炜. 渤海表层悬沙的时空分布特征及其动力成因[J]. 海洋科学进展, 2014, 32(4):450-458. [PANG Chongguang, LI Kun, YU Wei. Distribution characteristics, seasonal variability and dynamical mechanism of suspended sediment in the surface layer of the Bohai Sea [J]. Advances in Marine Science, 2014, 32(4): 450-458. doi: 10.3969/j.issn.1671-6647.2014.04.002 [17] Bian C W, Jiang W S, Greatbatch R J, et al. The suspended sediment concentration distribution in the Bohai Sea, Yellow Sea and East China Sea [J]. Journal of Ocean University of China, 2013, 12(3): 345-354. doi: 10.1007/s11802-013-1916-3
[18] 杜家笔, 裴艳东, 高建华, 等. 弱动力浅海中的悬沙输运机制: 以天津港附近海域为例[J]. 海洋学报, 2012, 34(1):136-144. [DU Jiabi, PEI Yandong, GAO Jianhua, et al. The suspended sediment transport associated with low flow patterns in shallow waters: a case study from the Tianjin subtidal area [J]. Acta Oceanologica Sinica, 2012, 34(1): 136-144. [19] 韩彬, 宋转玲, 曹磊, 等. 龙口湾近岸海域水质状况调查与评价[J]. 海洋科学进展, 2010, 28(2):186-192. [HAN Bin, SONG Zhuanling, CAO Lei, et al. Survey and assessment of coastal seawater quality in Longkou Bay [J]. Advances in Marine Science, 2010, 28(2): 186-192. doi: 10.3969/j.issn.1671-6647.2010.02.008 [20] Miao J B, Liu X Q. A numerical study of the wintertime circulation in the Northern Huanghai Sea and the Bohai Sea part I: Basic characteristics of the circulation [J]. Chinese Journal of Oceanology and Limnology, 1988, 6(3): 216-226. doi: 10.1007/BF02846499
[21] 蓝先洪, 蜜蓓蓓, 李日辉, 等. 渤海东部和黄海北部沉积物中重金属分布特征[J]. 中国环境科学, 2014, 34(10):2660-2668. [LAN Xianhong, MI Beibei, LI Rihui, et al. Distribution of heavy metals in surface sediments in the eastern Bohai Sea and the northern Yellow Sea [J]. China Environmental Science, 2014, 34(10): 2660-2668. [22] 赵保仁, 方国洪, 曹德明. 渤海、黄海和东海的潮余流特征及其与近岸环流输送的关系[J]. 海洋科学集刊, 1995(1):1-11. [ZHAO Baoren, FANG Guohong, CAO Deming. Characteristics of tidal residual current sand the irrelations with coastal current transports in the Bohai sea, Yellow sea and East China Sea [J]. Studia Marina Sinica, 1995(1): 1-11. [23] Yang S Y, Jung H S, Lim D I, et al. A review on the provenance discrimination of sediments in the Yellow Sea [J]. Earth-Science Reviews, 2003, 63(1-2): 93-120. doi: 10.1016/S0012-8252(03)00033-3
[24] Guan B X. Patterns and structures of the currents in Bohai, Huanghai and East China Seas[M]//Oceanology of China Seas. Dordrecht: Springer, 1994: 17-26.
[25] 韦钦胜, 于志刚, 冉祥滨, 等. 黄海西部沿岸流系特征分析及其对物质输运的影响[J]. 地球科学进展, 2011, 26(2):145-156. [WEI Qinsheng, YU Zhigang, RAN Xiangbin, et al. Characteristics of the western coastal current of the Yellow Sea and its impacts on material transportation [J]. Advances in Earth Science, 2011, 26(2): 145-156. [26] 刘成, 胡日军, 朱龙海, 等. 庙岛群岛海域沉积动力环境分区及沉积物输运趋势[J]. 海洋地质前沿, 2018, 34(8):24-33. [LIU Cheng, HU Rijun, ZHU Longhai, et al. Dynamic environment division and sediment transport trend in the area off Miaodao Islands [J]. Marine Geology Frontiers, 2018, 34(8): 24-33. [27] 马晓红, 韩宗珠, 艾丽娜, 等. 中国渤黄海的沉积物源及输运路径研究[J]. 中国海洋大学学报: 自然科学版, 2018, 48(6):96-101. [MA Xiaohong, HAN Zongzhu, AI Li’na, et al. Research on provenance and transport pattern in the Bohai Sea and Yellow Sea [J]. Periodical of Ocean University of China, 2018, 48(6): 96-101. [28] 冯兴如, 杨德周, 尹宝树. FVCOM在龙口海域潮汐潮流模拟中的应用研究[J]. 海洋科学, 2010, 34(6):94-99. [FENG Xingru, YANG Dezhou, YIN Baoshu. Application of FVCOM in tidal modeling of the seas adjacent to Longkou City [J]. Marine Sciences, 2010, 34(6): 94-99. [29] 王从敏, 彭马川, 张启龙. 龙口湾潮流场的数值模拟[J]. 海岸工程, 1996, 15(4):1-8. [WANG Congmin, PENG Machuan, ZHANG Qilong. Numerical simulation for the tidal current field in the Longkou bay [J]. Coastal Engineering, 1996, 15(4): 1-8. [30] 安永宁, 吴建政, 朱龙海, 等. 龙口湾冲淤特性对人工岛群建设的响应[J]. 海洋地质动态, 2010, 26(10):24-30. [AN Yongning, WU Jianzheng, ZHU Longhai, et al. Response of erosion-deposition pattern to artificial islands construction in Longkou bay [J]. Marine Geology Letters, 2010, 26(10): 24-30. [31] 吕双燕, 金秉福, 贺世杰, 等. 莱州湾-龙口湾表层沉积物有机质特征及来源分析[J]. 环境化学, 2017, 36(3):650-658. [LYU Shuangyan, JIN Bingfu, HE Shijie, et al. Characteristics and sources of organic matter in surface sediments of Laizhou Bay and Longkou Bay [J]. Environmental Chemistry, 2017, 36(3): 650-658. doi: 10.7524/j.issn.0254-6108.2017.03.2016070103 [32] 任鹏, 孙志高, 王传远, 等. 人工岛建设对龙口湾表层沉积物粒度及黏土矿物组成特征的影响[J]. 海洋科学进展, 2016, 34(4):578-587. [REN Peng, SUN Zhigao, WANG Chuanyuan, et al. Impacts of construction of artificial islands on the flow-sediment regulation scheme on grain and clay compositions in the Longkou Bay [J]. Advances in Marine Science, 2016, 34(4): 578-587. doi: 10.3969/j.issn.1671-6647.2016.04.014 [33] 中国海湾志编纂委员会. 中国海湾志(第三分册)[M]. 北京: 海洋出版社, 1991. China Bay Record Committee. The Bay Chorography in China: Bays in the South Shandong Peninsula and Jiangsu Province[M]. Beijing: China Ocean Press, 1993: 73-94.
[34] 杨晓东. 乐清湾悬沙输移特性研究[D]. 浙江大学硕士学位论文, 2010. YANG Xiaodong. Study on transport characteristics of suspended sediment in Yueqing bay[D]. Master Dissertation of Zhejiang University, 2010.
[35] 吴德安. 江苏辐射沙洲水道潮流及悬沙动力研究[D]. 南京师范大学博士学位论文, 2004. WU Dean. Hydrodynamic studies to the tidal current and suspended sediments in the channels of Jiangsu radial sand ridges[D]. Doctor Dissertation of Nanjing Normal University, 2004.
[36] 王文圣, 丁晶, 向红莲. 小波分析在水文学中的应用研究及展望[J]. 水科学进展, 2002, 13(4):515-520. [WANG Wensheng, DING Jing, XIANG Honglian. Application and prospect of wavelet analysis in hydrology [J]. Advances in Water Science, 2002, 13(4): 515-520. doi: 10.3321/j.issn:1001-6791.2002.04.021 [37] Venugopal V, Foufoula-Georgiou E. Energy decomposition of rainfall in the time-frequency-scale domain using wavelet packets [J]. Journal of Hydrology, 1996, 187(1-2): 3-27. doi: 10.1016/S0022-1694(96)03084-3
[38] Kumar P, Foufoula-Georgiou E. A multicomponent decomposition of spatial rainfall fields: 2. Self-similarity in fluctuations [J]. Water Resources Research, 1993, 29(8): 2533-2544. doi: 10.1029/93WR00549
[39] 胡日军, 吴建政, 朱龙海, 等. 东海舟山群岛海域表层沉积物运移特性[J]. 中国海洋大学学报: 自然科学版, 2009, 39(3):495-500, 442. [HU Rijun, WU Jianzheng, ZHU Longhai, et al. Characteristic of surface sediment transport in Zhoushan Archipelago Sea Area, East China Sea [J]. Periodical of Ocean University of China, 2009, 39(3): 495-500, 442. [40] Ingram R G. Characteristics of the Great Whale River Plume [J]. Journal of Geophysical Research: Oceans, 1981, 86(C3): 2017-2023. doi: 10.1029/JC086iC03p02017
[41] Uncles R J, Elliott R C A, Weston S A. Dispersion of salt and suspended sediment in a partly mixed estuary [J]. Estuaries, 1985, 8(3): 256-269. doi: 10.2307/1351486
[42] 范恩梅, 陈沈良, 张国安. 连云港海域水文泥沙运动特征[J]. 世界科技研究与发展, 2009, 31(4):703-707. [FAN Enmei, CHEN Shenliang, ZHANG Guoan. The hydrological and sediment characteristics in Lianyungang coastal waters [J]. World Sci-Tech R & D, 2009, 31(4): 703-707. doi: 10.3969/j.issn.1006-6055.2009.04.040 [43] 孔亚珍, 朱传芳, 丁平兴. 水流作用下的悬沙垂向扩散系数[J]. 泥沙研究, 2009(6):43-47. [KONG Yazhen, ZHU Chuanfang, DING Pingxing. Vertical diffusion coefficient of suspended sediment in flows [J]. Journal of Sediment Research, 2009(6): 43-47. doi: 10.3321/j.issn:0468-155X.2009.06.007 [44] 李远. 近期长江口北槽水沙特性及悬沙浓度垂向分布规律[D]. 华东师范大学硕士学位论文, 2018. LI Yuan. Recent characteristics of flow and suspended sediment and vertical sediment concentration distribution pattern in the North Passage of Yangtze Estuary[D]. Master Dissertation of East China Normal University, 2018.
[45] 鲍献文, 万修全, 吴德星, 等. 2000年夏末和翌年初冬渤海水文特征[J]. 海洋学报, 2004, 26(1):14-24. [BAO Xianwen, WAN Xiuquan, WU Dexing, et al. Hydrologic characteristics of the end summertime in 2000 and early wintertime in 2001 in the Bohai Sea [J]. Acta Oceanologica Sinica, 2004, 26(1): 14-24. [46] 陈斌, 黄海军, 严立文, 等. 小清河口附近海域泥沙运动特征及风场对泥沙输运的影响[J]. 海洋学报, 2009, 31(2):104-112. [CHEN Bin, HUANG Haijun, YAN Liwen, et al. Sediment movement and the effect of wind to sediment transport near Xiaoqing River Estuary [J]. Acta Oceanologica Sinica, 2009, 31(2): 104-112. [47] 杨晓东, 姚炎明, 蒋国俊, 等. 乐清湾悬沙输移机制分析[J]. 海洋通报, 2011, 30(1):53-59. [YANG Xiaodong, YAO Yanming, JIANG Guojun, et al. Study on the transport mechanism of suspended sediment in Yueqing Bay [J]. Marine Science Bulletin, 2011, 30(1): 53-59. doi: 10.3969/j.issn.1001-6392.2011.01.009 [48] 杨作升, 郭志刚, 王兆祥, 等. 黄东海陆架悬浮体向其东部深海区输送的宏观格局[J]. 海洋学报, 1992, 14(2):81-90. [YANG Zuosheng, GUO Zhigang, WANG Zhaoxiang, et al. The macro pattern of the suspended body of the Yellow Sea and the East China Sea to the eastern deep sea area [J]. Acta Oceanologica Sinica, 1992, 14(2): 81-90. [49] 孙效功, 方明, 黄伟. 黄、东海陆架区悬浮体输运的时空变化规律[J]. 海洋与湖沼, 2000, 31(6):581-587. [SUN Xiaogong, FANG Ming, HUANG Wei. Spatial and temporal variations in suspended particulate matter transport on the Yellow and East China Sea shelf [J]. Oceanologia Et Limnologia Sinica, 2000, 31(6): 581-587. doi: 10.3321/j.issn:0029-814X.2000.06.001 [50] Bian C W, Jiang W S, Quan Q, et al. Distributions of suspended sediment concentration in the Yellow Sea and the East China Sea based on field surveys during the four seasons of 2011 [J]. Journal of Marine Systems, 2013, 121-122: 24-35. doi: 10.1016/j.jmarsys.2013.03.013
[51] 冷星, 朱龙海, 胡日军. 山东半岛东部海域泥质区冬季悬浮泥沙时空变化及输运机制[J]. 中国海洋大学学报: 自然科学版, 2019, 49(4):106-116. [LENG Xing, ZHU Longhai, HU Rijun. The spatiotemporal change and transport mechanism of suspended sediment in the mud area of the eastern sea of Shandong Peninsulain winter [J]. Periodical of Ocean University of China, 2019, 49(4): 106-116. -
期刊类型引用(1)
1. 杨华臣，张建中. 拖缆与OBN资料联合成像域最小二乘逆时偏移成像. 地球物理学报. 2022(10): 4099-4110 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 2609
- HTML全文浏览量: 558
- PDF下载量: 33
- 被引次数: 3
