Magnetic properties of sediments and their implications for sedimentary dynamic environment in the middle and lower reaches of the Qiantang River
-
摘要: 利用2010年1月与2010年8月采集的钱塘江中下游河床沉积物粒度和磁性测量数据以及2010年8月测量的流速数据，分析了沉积物粒度和磁性的时空分布特征，探讨了粒度和磁性参数对沉积动力环境的指示意义。结果表明：（1）冬季河床沉积物以粉砂和黏土为主，夏季以粉砂和砂为主，沉积物粒度呈现从中游到富春江水库逐渐变细、近口段到河口段逐渐变粗的规律，并且冬季粒度总体上细于夏季，指示夏季以及河口段较强的水动力环境。（2）磁性矿物含量从中游到下游呈现逐渐减少的趋势，夏季磁性矿物含量比冬季高。磁性矿物中亚铁磁性矿物占主导地位，夏季亚铁磁性矿物的含量高于冬季。（3）χfd%、χARM、χARM/SIRM和χARM/χ可以作为这个区域冬季＜16 μm的细颗粒沉积物的代用指标，χfd%和χARM/χ可以作为夏季＜32 μm细颗粒沉积物的代用指标。（4）沉积物磁性参数χfd%、χARM、χARM/χ和χARM/SIRM同样呈现出从中游到富春江水库逐渐变大、而从近口段到河口段逐渐变小的趋势，这与沉积物磁性矿物晶粒对沉积动力环境的响应是密切相关的，其中参数χfd%与χARM更能有效地反映沉积动力环境。Abstract: Based on the particle size and magnetic characteristics of the sediments collected in the middle and lower reaches of the Qiantang River in January 2010 and August 2010, as well as the flow velocity data measured in August 2010, the temporal and spatial distribution of particle size and magnetism of the sediments are analyzed, and their indication significance to the sedimentary dynamic environment dicussed. Results show that the river bottom sediments are mainly composed of silt and clay in winter, but dominated by silt and sand in summer. The grain size of the sediments gradually gets fining from middle reaches to the Fuchunjiang Reservoir and becomes coarse from near-mouth section to estuary; the particle size is smaller in winter in general, indicating a stronger hydrodynamic environment in the estuary in summer. The content of magnetic minerals gradually decreases from the middle reaches to the lower reaches. The content of magnetic minerals in summer is higher than that in winter. Ferrimagnetic minerals dominate magnetic minerals, and there are more ferrimagnetic minerals in summer. χfd%,χARM, χARM/SIRM and χARM/χ can be used as proxies for the sediments finer than 16 μm in the study area in winter, while χfd% and χARM/χ be used as indicators for sediments finer than 32 μm in summer. χfd%, χARM, χARM/χ and χARM/SIRM also show an increasing trend from middle reaches to the Fuchunjiang Reservoir, while show a decreasing trend from near-mouth section to estuary. This is closely related to the response of the sediment magnetic crystals to the sedimentary dynamic environment. χfd% and χARM can reflect the sedimentary dynamic environment more effectively.
-
1. 研究背景
S变换[1]是由美国的地球物理学家Stockwell等[2-3]在短时傅里叶变换和小波变换的基础上提出的。与短时傅里叶变换和小波变换相比，S变换既解决了短时傅里叶变换时窗固定的问题，同时具有小波变换的多分辨率分析。由于S变换所具有的独特优势，其在信号分析和处理领域受到广泛关注。由于高斯时窗和频率的函数关系较为单一，降低了S变换在实际应用中的灵活性，因此很多学者从不同角度对S变换进行了推广和改进得到广义S变换[4-8]，使得广义S变换更加适用于地震资料处理、油气检测[9-11]等方面的信号处理。
广义S变换作为时频域的滤波工具，能够压制噪音从而提高地震资料信噪比[12]。赵淑红等[13]首先在频域和时域确定干扰波范围，再利用广义S变换处理VSP数据；李雪英等[14]考虑到地震信号有效频带随时间逐渐变窄，在时频域采用交互式的方法确定高频噪音切除边界；王云专等采用时频谱叠加的方法计算共中心点道集中的滤波器进而压制随机干扰[15]。陈学华等引入λ和p两个参数控制高斯窗函数的变化趋势，在应用过程中通过试验确定最优参数[5]。本文对陈学华[5]的广义S变换方法进行了改进，使其兼顾对高低频端信号的分析能力，能够更加适用于随机噪音压制；其次由于标准S逆变换存在能量泄露，通过修正S逆变换[16-19]，以消除滤波噪音。在压制噪音时，首先计算信号瞬时信噪比，根据瞬时信噪比的高低结合不同处理策略实现噪音压制。通过处理合成数据和实际资料，验证改进广义S正反变换在随机噪音压制方面的有效性。
2. 方法原理
2.1 改进广义S变换
信号u(t)的标准S正变换的表达式为
$$ \int S\left(\tau ,f\right)={\int }_{-\infty }^{\infty }u\left(t\right)w(\tau -t,f){\mathrm{e}}^{-i2{{\text π} }ft}\mathrm{d}t $$ (1) 式中，
$ S\left(\tau ,f\right) $ 是信号$ u\left(t\right) $ 的时频谱，$ w(t,f) $ 是高斯窗函数，其具体表达式为$$ w\left(t,f\right)=\frac{1}{\sigma \sqrt{2\mathrm{{\text π} }}}{\mathrm{e}}^{\frac{-{t}^{2}}{2{\sigma }^{2}}} $$ (2) 式中，
$ \sigma =\dfrac{1}{\left|f\right|} $ 是控制高斯窗宽度的尺度因子。广义S变换就是通过调节尺度因子大小控制对应的高斯窗函数的有效宽度。针对标准S变换高斯窗形态随频率变化趋势单一进而导致S变换时频谱聚焦度低的问题，陈学华等[5]在标准S变换的基础上引入了控制参数λ和p，目的为丰富尺度因子随频率的变化趋势，针对不同的资料提供更多的参数组合。其改进的尺度因子为
$$ \sigma =\frac{1}{\lambda {\left|f\right|}^{p}} $$ (3) 不同的控制参数λ和p影响尺度因子变化的快慢（图1）。不同的λ值对尺度因子的主要影响集中在低频端，即随着λ值的增大，同一频率对应的尺度因子幅值增大。不同的p值对尺度因子的主要影响在高频端，即随着p值的增大，同一频率对应的尺度因子幅值减小。为了进一步说明两个参数对时变信号的聚焦作用，选取理论信号进行效果分析。该信号是由两个不同的宽频带调频信号组成，即使得高斯窗尺度因子随频率的变化范围比较大，可以有效检验方法的适用性。
![]()
图 1 理想时频谱与不同控制参数的广义S变换的时频谱a. 理想时频；b. p=0.7，λ=1对应的时频谱；c. p=0.9, λ=1对应的时频谱；d. p=1，λ=1标准S变换时频谱；e. p=1.2，λ=1对应的时频谱；f. p=1，λ=1.2对应的时频谱；g. p=1，λ=0.9对应的时频谱；h. p=1，λ=0.6对应的时频谱。Figure 1. Ideal time-frequency spectrum and time-frequency spectrum of S transform using different control factorsa. ideal time-frequency spectrum, b. p=0.7, λ=1 corresponding to time-frequency spectrum, c. p=0.9, λ=1 corresponding to time-frequency spectrum, d. p=1, λ=1 corresponding to time-frequency spectrum, e. p=1.2, λ=1 corresponding to time-frequency spectrum, f. p=1, λ=1.2 corresponding to time-frequency spectrum, g. p=1, λ=0.9 corresponding to time-frequency spectrum, h. p=1, λ=0.6 corresponding to time-frequency spectrum.$$\begin{split} y=&{\sin}\left\{3*\left[140{{\text π} }t+40{\sin}\left(3{{\text π} }t\right)\right]\right\}+\\&{\sin}\left\{3*\left[140{{\text π} }t+40{\cos}\left(3{{\text π} }t+120\right)\right]\right\} \end{split}$$ (4) 采用不同控制参数λ和p对上述信号进行广义S变换，分析对应的时频谱聚焦程度。当p值由小变大时，信号时频谱在低频端与高频端的聚焦程度均响应明显。其中低频端的聚焦程度随p值增大而变好，高频端随p值增大而变差，甚至当p＞1时，高频端出现严重失焦（图1a-d）。当p值由小变大时，信号时频谱在低频端的聚焦程度无响应，而高频端的聚焦程度也变化不大（图1c、f、g、h）。广义S变换中参数λ和p只能取固定值，其单一性的变化趋势不能兼顾高低频端的信号聚焦程度。综合以上实验结果可知，p值对时频域聚焦程度影响更大。改进广义S变换主要是更改p值变化趋势，而固定λ=1。
p值由小变大时，低频端与高频端聚焦程度响应不同。高频端信号需要小的p值，而低频端信号需要大的p值。为实现广义S变换在时频域的宽频信号扫描，将参数p改为频率的线性函数，即
$$ p=p\left(f\right)=a+bf $$ (5) 其中，初始值a＞0，变化率b＜0。
在S变换计算过程中，计算的最小频率为0 Hz，最大频率为待分析信号的Nyquist频率，在这个频带范围内规定p值线性递减，同时调整p值的初始值a和变化率b，使p值能够自动匹配信号，增加信号时频分析的聚焦程度。
为了检验该方法的效果，以理论模拟信号为例进行时频分析对比。该信号是由两个频率变化趋势相反的线性调频信号和一个正弦调频信号组成，见公式（6）。图2是标准S变换、广义S变换的优选方案（p=0.8，λ=1）和改进广义S变换的时频谱图（p=0→1，λ=1）。由图2可知，标准S变换在频率方向出现了条带状能量团，模糊了真实时频谱。广义S变换改善了时频谱图，增强了时频谱的聚焦度，而改进的广义S变换其高频的线性调频信号时频谱聚焦度更高，正弦型调频信号能量分布更均匀，分析效果要优于广义S变换。
![]()
图 2 合成信号时频谱图a. 标准S变换，b. 广义S变换（p=0.8，λ=1），c. 改进的广义S变换（p=0→1，λ=1）。Figure 2. Synthetic signal time-frequency spectruma. the standard S transform, b. conventional generalized S transform (p=0.8, λ=1), c. improved generalized S transform (p=0→1, λ=1).$$\begin{split} y=&{\cos}\left(68{{\text π} }t-20{{\text π} }{t}^{2}\right)+{\cos}\left[2{{\text π} }{\sin}\left(5{{\text π} }t\right)+120{{\text π} }t\right]+\\&{\cos}(168{{\text π} }t+28{{\text π} }{t}^{2}) \end{split}$$ (6) 为了进一步检验改进的广义S变换在地震信号分析中的适用性，选取一道添加高斯白噪声的实际地震信号进行分析。在高频噪声存在的情况下，利用上述三种方法得到时频谱的能量分布（图3）。由图3中立体图与平面图对比可知，标准S变换对低频部分的信号分析结果聚焦度比较好，但其高频噪音的时频谱的能量被过分放大，如果用这种方法在时频域采用阈值滤波的方法进行去噪，显然会把部分信号当成噪音同时也无法把高频噪音去掉。广义S变换的时频谱改善了高频噪音部分的聚焦度，但在低频端能量团几乎连在一起，时间分辨率变差。改进的广义S变换和标准S变换在低频端的能量分布相近，在高频端使高频噪音的能量分散聚焦，而且各个时频点的能量比较稳定，没有被过分放大。可以看出改进的广义S变换对地震资料进行分析或者时频域滤波，效果显然要优。
![]()
图 3 含随机噪声的地震信号时频谱图a. 标准S变换，b. 广义S变换（p=0.8, λ=1），c. 改进的广义S变换（p=0→1, λ=1）。Figure 3. Time-frequency spectrum of seismic signal with random noisea. the standard S transform, b. conventional generalized S transform (p=0.8, λ=1), c. improved generalized S transform (p=0→1, λ=1) .2.2 改进S域时频滤波
由于广义S变换具有逆变换，因此，可以利用广义S正反变换实现时频滤波。其常规滤波思路是利用二维滤波函数F(τ，f )与S变换时频谱相乘，对乘积结果沿时间轴积分，再反傅氏变换回时间域。滤波公式为：
$$ {u}_{\mathrm{f}\mathrm{i}\mathrm{l}\mathrm{t}1}\left(t\right)={\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}{S}(\tau ,f){F}(\tau ,f){\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}\tau \mathrm{d}f $$ (7) 式中，F(τ, f )是时频域滤波函数。通过计算发现当F(τ, f ) =1时，即不对时频谱作任何处理可以完全重构信号，但当F(τ, f ) ≠1时，实际滤波结果和理想滤波结果存在差异，无法保证时频滤波的准确性[10]，如图4所示，在时频域将合成信号为500～600 ms的时频谱滤掉，反变换回时间域就出现了能量泄露。
![]()
图 4 标准S逆变换时频滤波效果图a. 合成信号，b. 合成信号时频域滤波，c. 时间域滤波效果。Figure 4. Effect of time and frequency filtering of standard S inverse transforma. synthetic signal, b. time and frequency filtering of synthetic signal, c. filtering effect of time field.为了解决这个问题，需要修正变换公式得到精确滤波的改进S域时频滤波公式。
$$\begin{split} S\left(\tau ,f\right)=&{\int }_{-\infty }^{\infty }u\left(t\right)w\left(\tau -t,f\right){\mathrm{e}}^{-i2\mathrm{{\text π} }ft}\mathrm{d}t=\\&{\int }_{-\infty }^{\infty }u\left(t\right)\frac{\lambda {\left|f\right|}^{p}}{\sqrt{2\mathrm{{\text π} }}}{\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{(\tau -t)}^{2}}{2}}{\mathrm{e}}^{-i2\mathrm{{\text π} }ft}\mathrm{d}t\end{split}$$ (8) 式中，u(t)是待分析信号，w(τ-t, f)是窗函数。
根据高斯窗函数的表达式，引入一个新的变量x(τ, t)。
$$ x\left(\tau ,t\right)=u\left(t\right){\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{(\tau -t)}^{2}}{2}} $$ 对式（8）针对变量t做傅里叶变换：
$$ x\left(\tau ,f\right)={\int }_{-\infty }^{\infty }u\left(t\right){{\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{(\tau -t)}^{2}}{2}}\mathrm{e}}^{-i2\mathrm{{\text π} }ft}\mathrm{d}t $$ (9) 对比式（7）和式（9）可得
$$ x\left(\tau ,f\right)=\frac{\sqrt{2\mathrm{{\text π} }}}{\lambda {\left|f\right|}^{p}}S\left(\tau ,f\right) $$ (10) 对式（10）针对变量 f 作反傅里叶变换：
$$ x\left(\tau ,t\right)=\frac{\sqrt{2\mathrm{{\text π} }}}{\lambda }{\int }_{-\infty }^{\infty }\frac{S\left(\tau ,f\right)}{{\left|f\right|}^{p}}{\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}f $$ (11) 由式（8）可知，当τ=t时，u(t)和x(τ, t)之间具有对应关系：
$$ u\left(t\right)=x(t,t) $$ (12) 因此，令τ=t代入式（11）可得：
$$ u\left(t\right)=x\left(t,t\right)=\frac{\sqrt{2\mathrm{{\text π} }}}{\lambda }{\int }_{-\infty }^{\infty }\frac{S\left(t,f\right)}{{\left|f\right|}^{p}}{\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}f $$ (13) 对u(t)进一步展开：
$$ \begin{split}u\left(t\right)=&\frac{\sqrt{2\mathrm{{\text π} }}}{\lambda }{\int }_{-\infty }^{\infty }\frac{S\left(t,f\right)}{{\left|f\right|}^{p}}{\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}f=\\&{\int }_{-\infty }^{\infty }u\left(t\right){\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{\left(\tau -t\right)}^{2}}{2}}{\mathrm{e}}^{-i2\mathrm{{\text π} }ft}\mathrm{d}t=\\&{\int }_{-\infty }^{\infty }u\left(\tau \right){\int }_{-\infty }^{\infty }{\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{\left(t-\tau \right)}^{2}}{2}}{\mathrm{e}}^{-i2\mathrm{{\text π} }f(t-\tau )}\mathrm{d}f\mathrm{d}\tau =\\&{\int }_{-\infty }^{\infty }u\left(\tau \right)m\left(t-\tau \right)\mathrm{d}\tau =u\left(t\right)\otimes m\left(t\right)\end{split} $$ (14) 式中，
$ \otimes $ 代表褶积，函数m(t)的表达式：$$ m\left(t\right)={\int }_{-\infty }^{\infty }{\mathrm{e}}^{\frac{-{\lambda }^{2}{f}^{2p}{t}^{2}}{2}}{\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}f $$ (15) 由式（14）可知，计算得到的结果是原始信号和信号m(t)的褶积，因此在时频域滤波结果的基础上需要再做一步滤波修正处理，最终的滤波表达式为：
$$ {u}_{filt2}\left(t\right)=\frac{\sqrt{2\mathrm{{\text π} }}}{\lambda }{\int }_{-\infty }^{\infty }\frac{S(t,f)F(t,f)}{{\left|f\right|}^{p}}{\mathrm{e}}^{i2\mathrm{{\text π} }ft}\mathrm{d}f $$ (16) $$ {u}_{\mathrm{f}\mathrm{i}\mathrm{l}\mathrm{t}3}\left(t\right)={\rm{IFT}}\left(\frac{{U}_{\mathrm{f}\mathrm{i}\mathrm{l}\mathrm{t}2}\left(f\right)}{M\left(f\right)}\right) $$ (17) 式中，F(t, f)是时频域滤波函数，Ufilt2(f)是ufilt2(t)的频谱，M(f)是m(t)的频谱，IFT表示反傅里叶变换。式（16）和式（17）是最终的滤波公式。
为了检验对S逆变换的改进效果，分别用标准逆变换、改进逆变换将图4b中的时频谱逆变换，并计算与准确滤波结果的差值，如图5所示。对比结果显示，改进逆变换消除了滤波区信号的畸变，避免产生滤波噪音。
![]()
图 5 准确滤波结果和不同滤波方法计算的差值a. 准确滤波结果，b. 标准逆变换计算差值（红色）和改进逆变换（蓝色）计算差值。Figure 5. Accurate filtering result and difference calculated by different filtering methodsa. accurate filtering result, b. difference calculated by standard S inverse transform (red) and difference calculated by improved S inverse transform (blue).基于上述改进广义S变换和修正的S逆变换，可以将随机噪音压制概括为3步：① 对含噪信号进行改进广义S变换得到信号的S谱；② 在S谱的基础上利用谱比法[20]计算信号的瞬时信噪比；③ 根据瞬时信噪比的不同，采用指数衰减和阈值滤波两种策略压制随机噪音。具体处理步骤如下：
（1）计算信号的S变换时频谱，利用该时频谱计算瞬时信噪比和阈值。阈值的计算思路是：计算初至前两个频率段内噪音的时频谱的平均值，作为初始阈值，具体的频率段可以分成有效频带内和有效频带外两段，计算两个阈值。由于高频随机干扰的时频谱幅值的均值随时间变化不大，而相对浅层数据深层数据的高频成分衰减快，能量弱，高频噪音能量占优，因此在无法利用初至前记录计算阈值时也可利用深层数据较高频段的时频谱计算初始阈值。
$$ {E}_{\mathrm{m}\mathrm{e}\mathrm{a}\mathrm{n}}\left(t,f\right)=\frac{1}{T\times F}\sum _{t={t}_{1}}^{{t}_{2}}\sum _{f={f}_{1}}^{{f}_{2}}\left|S(t,f)\right| $$ (18) 式中，Emean(t, f)是阈值，S(t, f)是时频谱，T是参与计算时间点个数，F是参与计算频率点个数，t1到t2是参与计算的时间段，f1到f2是参与计算的频率段。
（2）实际地震数据需要考虑地层的吸收衰减，且高频衰减更快，因此深部有效频带变窄，有效信号能量弱[21]，此时如果将深部地震记录按照阈值处理，会损失很多有效信号，因此为了保留一部分有效信号以供后续处理和进一步使用，此时对瞬时信噪比低的数据段采用指数衰减的方法处理，即对信噪比小于1的时间点，对该时间点大于有效频带高频端的S谱的幅值进行一个指数衰减的加权，即将S(t, f)乘以一个以f为变量的指数衰减函数H(f)，如式（19）：
$$ H\left(f\right)={\mathrm{e}}^{-\frac{f}{\mu }} $$ (19) 式中，μ是一个衰减系数，f是频率，μ越大对高频衰减越小，反之对高频衰减越大。
（3）当瞬时信噪比大于1时，有效信号能量占优，且随机噪音的时频谱分布在整个频带上，对该时间点所有S谱采用阈值滤波法，当某频率点的幅值小于r倍的Emean(t, f)的时候，定义该点的S谱为噪音，将其充零，当其幅值大于r倍的Emean(t, f)的时候，将其时频谱的幅值减去r倍的Emean(t, f)，并将处理之后的时频谱反变换回时间域。
在实际应用中，应根据实际地震资料的情况，进行反复试验，选取合适的μ和r值，既能压制随机噪音，又能保护高频端的有效信号，为数据拓频留有空间。
3. 合成信号试验
为了检验方法的有效性，首先利用零相位雷克子波模拟48道经地层吸收衰减的地震记录，加入随机白噪声后，利用改进广义S变换对地震记录压制随机噪音（图6）。抽取单道地震模拟数据对比，能够清晰地看出改进广义S变换压制随机噪声的效果，尤其深部低信噪比信号也得到有效的恢复（图7）。分析时频域中滤波效果，可看到有效信号时频谱得以保留，随机噪音时频谱得到有效压制（图8）。
![]()
图 6 含噪地震记录去噪效果图a. 48道模拟地震记录， b. 含白噪声的模拟地震记录， c. 去除噪声的模拟地震记录。Figure 6. Denoising effect for seismic records with noisea. 48 channels simulated seismic record, b. simulated seismic record with white noise, c. simulated seismic record after denoising.![]()
图 7 单道信号对比图a. 加白噪声后与未加噪声单道信号对比， b. 去噪后与未加噪声单道信号对比。Figure 7. Comparison of single channel signalsa. comparison of single channel signals with noise and those without noise, b. comparison of single channel signals without noise and those after noise reduction.![]()
图 8 时频域滤波效果对比a. 不含白噪声模拟地震道的时频谱，b. 含白噪声模拟地震道的时频谱，c. 去除噪声模拟地震道的时频谱。Figure 8. Comparison of filtering effects in time-frequency domaina. time and frequency spectrum of simulated seismic channel without noise, b. time and frequency spectrum of simulated seismic channel with white noise, c. time and frequency spectrum of simulated seismic channel after noise reduction.4. 实际资料处理
为了检验本文改进的广义S变换方法在实际资料中的应用效果，选择高分辨率海洋二维地震叠前剖面进行去噪处理，并将其与传统广义S变换滤波去噪法进行对比。本文首先通过单炮初至前的噪声信号确定去噪阈值，再利用改进的广义S变换进行时频滤波。图9为地震叠前单炮剖面。传统广义S变换参数经过对比后确定为（p=0.6, λ=1），改进后的广义S变换应用参数（p=0.6→1, λ=1）。图10是两种方法压制噪声效果图。从整体上可以看出，两种广义S变换时频域去噪均能达到一定效果，随机噪声得到合理去除。但改进的广义S变换在远偏移距处的弱信号去噪效果更好，深部反射同相轴连续性更优，有效信号的损失较小，对弱信号起到保护作用。
![]()
图 10 两种方法压制噪声效果对比a. 改进广义S变换压噪后剖面，b. 传统广义S变换压噪后剖面，c. 改进广义S变换去除的噪声，d. 传统广义S变换去除的噪声。Figure 10. Comparison of noise suppression effects between two methodsa. denoised profile by improved generalized S transform, b. denoised profile by traditional generalized S transform, c. noise removal by improved generalized S transform, d. noise removal by generalized S transform.5. 结论
本文针对广义S变换高斯窗参数λ和p固定、逆变换存在信号泄露的问题，提出了改进的思路和计算过程，并将其应用于随机噪音压制。通过对合成信号和实际资料的应用，得到以下认识：
（1）改进广义S变换时频谱聚焦度要优于广义S变换，而且参数p与频率f呈线性变化，以适应不同的应用目的；改进后的逆变换可以实现精确滤波，消除了滤波噪音。
（2）利用单炮初至前或深部较高频段时频谱计算阈值，根据每个时刻时频谱的幅值与阈值的差异压制噪音，同时考虑到地震信号的有效频带和瞬时信噪比随时间变化的情况，对瞬时信噪比较低的采用指数衰减的策略压制噪音，可以在保护有效信号的基础上实现去噪目的。
（3）高斯窗参数
$ p $ 的范围选择，应当根据高低频端的聚焦程度而改变。通过试验确定的参数$ p $ ，既能够有效压制噪音又能避免弱有效信号的损失，为后续处理打好基础。 -
![]()
图 2 钱塘江中下游流域河床沉积物中值粒径沿程变化
Figure 2. Changes of median grain size of bed sediments along the middle and lower reaches of the Qiantang River
![]()
图 3 钱塘江中下游流域河床沉积物质量磁化率（χ）沿程变化
Figure 3. Changes of mass magnetic susceptibility （χ） of bed sediments along the middle and lower reaches of the Qiantang River
![]()
图 4 钱塘江中下游流域河床沉积物饱和等温剩磁（SIRM）沿程变化
Figure 4. Changes of saturation isothermal remanent magnetization (SIRM) of bed sediments along the middle and lower reaches of the Qiantang River
![]()
图 5 钱塘江中下游流域河床沉积物饱和硬剩磁（HIRM）沿程变化
Figure 5. Changes of hard remanent magnetization (HIRM) of bed sediments along the middle and lower reaches of the Qiantang River
![]()
图 6 磁性参数
$ \chi $ 和SIRM之间的相关关系Figure 6. Correlation between magnetic parameter
$ \chi $ and SIRM![]()
图 7 钱塘江中下游流域冬季河床沉积物S-100插值图
Figure 7. S-100 interpolation map of bed sediments in the middle and lower reaches of the Qiantang River in winter
![]()
图 8 钱塘江中下游流域夏季河床沉积物S-100插值图
Figure 8. S-100 interpolation map of bed sediments in the middle and lower reaches of the Qiantang River in summer
![]()
图 9 钱塘江中下游流域冬季
$ {\chi }_{\rm{fd \%}} $ 、$ {\chi }_{\rm{ARM}} $ 、$ {\chi }_{\rm{ARM}} $ /$ \chi $ 和$ {\chi }_{\rm{ARM}} $ /SIRM参数的沿程变化Figure 9. Changes of
$ {\chi }_{\rm{fd \%}} $ ,$ {\chi }_{\rm{ARM}} $ ,$ {\chi }_{\rm{ARM}} $ /$ \chi $ and$ {\chi }_{\rm{ARM}} $ /SIRM parameters along the middle and lower reaches of the Qiantang River in winter![]()
图 10 钱塘江中下游流域夏季
$ {\chi }_{\rm{fd \%}} $ 、$ {\chi }_{\rm{ARM}} $ 、$ {\chi }_{\rm{ARM}} $ /$ \chi $ 和$ {\chi }_{\rm{ARM}} $ /SIRM参数的沿程变化Figure 10. Changes of
$ {\chi }_{\rm{fd \%}} $ ,$ {\chi }_{\rm{ARM}} $ ,$ {\chi }_{\rm{ARM}} $ /$ \chi $ and$ {\chi }_{\rm{ARM}} $ /SIRM parameters along the middle and lower reaches of the Qiantang River in summer![]()
图 11 钱塘江中下游流域夏季流速大小分布图
Figure 11. Flow velocity distribution of the middle and lower reaches of the Qiantang River in summer
![]()
图 12 钱塘江中下游冬夏季河床沉积物各组分含量分布图
Figure 12. Content distribution of each component of bed sediments in the middle and lower reaches of the Qiantang River in winter and summer
表 1 钱塘江中下游流域河床沉积物粒度参数统计
Table 1 Statistics of grain size parameters of bed sediments in the middle and lower reaches of the Qiantang River
位置 季节 砂/% 粉砂/% 黏土/% 中值粒径/Φ 河口段 冬季 26.23 65.45 6.99 4.53 夏季 38.23 55.25 5.52 4.33 近口段 冬季 9.33 61.33 25.07 6.35 夏季 25.74 58.22 13.58 4.95 中游段 冬季 12.31 56.00 27.58 6.70 夏季 10.47 60.44 25.11 6.59 
下载: 导出CSV
表 2 钱塘江中下游流域冬季河床沉积物磁性参数与粒度参数的相关性
Table 2 Correlations between magnetic parameters and grain size characteristics of bed sediments in the middle and lower reaches of the Qiantang River in winter
粒级/μm $ \chi $ SIRM HIRM $ {\chi }_{\rm{fd \%}} $ $ {\chi }_{\rm{ARM}} $ $ {\chi }_{\rm{ARM}} $/SIRM $ {\chi }_{\rm{ARM}} $/$ \chi $ ＜2 0.211 0.198 0.183 0.494** 0.580** 0.413** 0.556** ＜4 0.182 0.180 0.189 0.551** 0.604** 0.471** 0.609** ＜8 0.159 0.169 0.209 0.604** 0.632** 0.521** 0.654** ＜16 0.164 0.192 0.231* 0.620** 0.653** 0.529** 0.662** ＜32 0.179 0.205 0.204 0.570** 0.618** 0.496** 0.605** ＜63 0.125 0.082 0.008 0.308** 0.303** 0.325** 0.283* 2～4 0.138 0.143 0.194 0.611** 0.615** 0.535** 0.662** 4～8 0.107 0.141 0.234* 0.671** 0.650** 0.588** 0.706** 8～16 0.170 0.243* 0.269* 0.576** 0.637** 0.483** 0.593** 16～32 0.029 0.018 −0.165 −0.369** −0.314** −0.283* −0.410** 32～63 −0.164 −0.261* −0.345** −0.623** −0.712** −0.470** −0.714** 砂 −0.125 −0.082 −0.008 −0.308** −0.303** −0.325** −0.283** 粉砂 0.030 −0.084 −0.192 −0.189 −0.247** −0.073** −0.279** 黏土 0.180 0.182 0.199 0.564** 0.611** 0.475** 0.617** 平均粒径 0.129 0.131 0.124 0.447** 0.482** 0.443** 0.487** 中值粒径 0.143 0.143 0.141 0.489** 0.478** 0.418** 0.454** 注：**表示在0.01水平（双侧）上显著相关；*表示在0.05水平（双侧）上显著相关。
下载: 导出CSV
表 3 钱塘江中下游流域夏季河床沉积物磁性参数与粒度参数的相关性
Table 3 Correlations between magnetic parameters and grain size characteristics of bed sediments in the middle and lower reaches of the Qiantang River in summer
粒级/μm $ \chi $ SIRM HIRM $ {\chi }_{\rm{fd \%}} $ $ {\chi }_{\rm{ARM}} $ $ {\chi }_{\rm{ARM}} $/SIRM $ {\chi }_{\rm{ARM}} $/$ \chi $ ＜2 0.222 0.318* 0.420** 0.386** 0.190 0.252* 0.406** ＜4 0.219 0.302* 0.409** 0.410** 0.190 0.267* 0.417** ＜8 0.214 0.280* 0.396** 0.429** 0.183 0.280* 0.422** ＜16 0.182 0.252* 0.384** 0.491** 0.179 0.313* 0.436** ＜32 0.120 0.194 0.362** 0.517** 0.168 0.337** 0.428** ＜63 −0.028 −0.009 0.230 0.393** 0.105 0.289* 0.290** 2～4 0.214 0.281* 0.393** 0.436** 0.189 0.282* 0.426** 4～8 0.193 0.230 0.349** 0.428** 0.160 0.281* 0.400** 8～16 0.074 0.143 0.021 0.560** 0.141 0.340** 0.398** 16～32 −0.216 −0.178 0.023 0.269* 0.005 0.206 0.102 32～63 −0.289* −0.404** −0.317* −0.342** −0.151 −0.165 −0.346** 砂 0.025 0.004 −0.232 −0.397** −0.106 −0.290* −0.293* 粉砂 −0.184 −0.213 0.047 0.271* 0.019 0.222 0.124 黏土 0.221 0.302* 0.408** 0.411** 0.191 0.267* 0.418** 平均粒径 0.111 0.153 0.349** 0.436** 0.158 0.297* 0.371** 中值粒径 0.109 0.151 0.340** 0.414** 0.141 0.279* 0.354** 注：**表示在0.01水平（双侧）上显著相关；*表示在0.05水平（双侧）上显著相关。
下载: 导出CSV
-
[1] Krom M D, Stanley J D, Cliff R A, et al. Nile River sediment fluctuations over the past 7000 yr and their key role in sapropel development [J]. Geology, 2002, 30(1): 71-74. doi: 10.1130/0091-7613(2002)030<0071:NRSFOT>2.0.CO;2
[2] Revel M, Ducassou E, Grousset F E, et al. 100, 000 years of African monsoon variability recorded in sediments of the Nile margin [J]. Quaternary Science Reviews, 2010, 29(11-12): 1342-1362. doi: 10.1016/j.quascirev.2010.02.006
[3] Montero-Serrano J C, Bout-Roumazeilles V, Sionneau T, et al. Changes in precipitation regimes over North America during the Holocene as recorded by mineralogy and geochemistry of Gulf of Mexico sediments [J]. Global and Planetary Change, 2010, 74(3-4): 132-143. doi: 10.1016/j.gloplacha.2010.09.004
[4] Yang S L, Xu K H, Milliman J D, et al. Decline of Yangtze River water and sediment discharge: Impact from natural and anthropogenic changes [J]. Scientific Reports, 2015, 5(1): 12581.
[5] 杨云平, 张明进, 李松喆, 等. 三峡大坝下游粗细颗粒泥沙输移规律及成因[J]. 湖泊科学, 2017, 29(4):942-954. [YANG Yunping, ZHANG Mingjin, LI Songzhe, et al. Transport patterns of the coarse and fine sediments and its causes in the downstream of the Three Gorges Dam [J]. Journal of Lake Sciences, 2017, 29(4): 942-954. doi: 10.18307/2017.0418 [6] 肖晓. 南海北部湾底质沉积物粒度和泥沙运移趋势研究[D]. 中国海洋大学硕士学位论文, 2015. XIAO Xiao. Surface sediment grain size and transport of the Beibu gulf in South China Sea[D]. Master Dissertation of Ocean University of China, 2015.
[7] Maher B A. Environmental magnetism and climate change [J]. Contemporary Physics, 2007, 48(5): 247-274. doi: 10.1080/00107510801889726
[8] Liu Q S, Roberts A P, Larrasoaña J C, et al. Environmental magnetism: Principles and applications [J]. Reviews of Geophysics, 2012, 50(4): RG4002.
[9] 时连强, 李九发, 张卫国, 等. 黄河三角洲飞雁滩HF孔沉积物的磁性特征及其环境意义[J]. 海洋学研究, 2007, 25(4):13-23. [SHI Lianqiang, LI Jiufa, ZHANG Weiguo, et al. Magnetic properties of Core HF from Feiyantan tidal flat, the Huanghe River Delta and its environmental significance [J]. Journal of Marine Sciences, 2007, 25(4): 13-23. doi: 10.3969/j.issn.1001-909X.2007.04.002 [10] 彭俊, 陈洪全, 马随随, 等. 黄河三角洲潮滩沉积物磁性特征与沉积环境分析[J]. 地理科学, 2014, 34(10):1262-1269. [PENG Jun, CHEN Hongquan, MA Suisui, et al. Magnetic properties of sediment and sedimentary environment in tidal flat of the Yellow River Delta [J]. Scientia Geographica Sinica, 2014, 34(10): 1262-1269. [11] 潘大东, 王张华, 陈艇, 等. 长江口表层沉积物矿物磁性分区特征及其沉积环境指示意义[J]. 海洋学报, 2015, 37(5):101-111. [PAN Dadong, WANG Zhanghua, CHEN Ting, et al. Mineral magnetic characteristics of surficial sediments and their implications for identifying sedimentary environments at the Changjiang River mouth [J]. Haiyang Xuebao, 2015, 37(5): 101-111. [12] 李金婵, 陈秀玲, 方红, 等. 福州河道表层沉积物磁学特征及其环境意义[J]. 地球环境学报, 2015, 6(1):17-25. [LI Jinchan, CHEN Xiuling, FANG Hong, et al. Magnetic properties of river sediments in Fuzhou and their environmental significance [J]. Journal of Earth Environment, 2015, 6(1): 17-25. doi: 10.7515/JEE201501003 [13] 杨小强, Grapes R, 周厚云, 等. 珠江三角洲沉积物的岩石磁学性质及其环境意义[J]. 中国科学D辑: 地球科学, 2008, 51(1):56-66. [YANG Xiaoqiang, Grapes R, ZHOU Houyun, et al. Magnetic properties of sediments from the Pearl River Delta, South China: Paleoenvironmental implications [J]. Science in China Series D: Earth Sciences, 2008, 51(1): 56-66. doi: 10.1007/s11430-007-0151-4 [14] 钱塘江志编纂委员会. 钱塘江志[M]. 北京: 方志出版社, 1998. Qiantang River Records Compilation Committee. Qiantang River Journal[M]. Beijing: Local Records Press, 1998.
[15] 刘朝. 钱塘江流域河流表层沉积物特征及物源分析[D]. 华东师范大学硕士学位论文, 2016. LIU Chao. Study on the characteristics of surface sediment in the Qiantang River and analysis of provenance[D]. Master Dissertation of East China Normal University, 2016.
[16] 中国海湾志编纂委员会. 中国海湾志(第十四分册)[M]. 北京: 海洋出版社, 1998. China Gulf Records Compilation Committee. Chinese Gulf Journal (Part 14)[M]. Beijing: China Ocean Press, 1998.
[17] Folk R L, Ward W C. Brazos River Bar: a study on the significance of grain size parameters [J]. Journal of Sedimentary Research, 1957, 27(1): 3-26. doi: 10.1306/74D70646-2B21-11D7-8648000102C1865D
[18] 赖智荣. 钱塘江河流表层沉积物矿物特征及其物源指示意义[D]. 华东师范大学硕士学位论文, 2019. LAI Zhirong. Study on the characteristics of minerals of surface sediment in the Qiantang River and analysis of provenance[D]. Master Dissertation of East China Normal University, 2019.
[19] 张卫国, 俞立中. 长江口潮滩沉积物的磁学性质及其与粒度的关系[J]. 中国科学D辑: 地球科学, 2003, 46(9):954-966. [ZHANG Weiguo, YU Lizhong. Magnetic properties of tidal flat sediments of the Yangtze Estuary and its relationship with particle size [J]. Science in China Series D: Earth Sciences, 2003, 46(9): 954-966. doi: 10.1007/BF02991341 [20] 褚慧敏, 周立旻, 黄静, 等. 岷江上游干流边滩沉积物岩石磁学特征及其影响因素[J]. 海洋地质与第四纪地质, 2016, 36(4):57-66. [CHU Huimin, ZHOU Limin, HUANG Jing, et al. Rock-magnetic properties of the point bat deposits in the upper mainstream of Minjiang River and their origin [J]. Marine Geology & Quaternary Geology, 2016, 36(4): 57-66. [21] 董瑞斌, 张卫国, 卢升高, 等. 土壤和沉积物的磁参数及其在环境科学中的应用[J]. 科技通报, 2000, 16(6):479-483. [DONG Ruibin, ZHANG Weiguo, LU Shenggao, et al. The magnetic indexes for soil and sediment and their applications on environmental studies [J]. Bulletin of Science and Technology, 2000, 16(6): 479-483. doi: 10.3969/j.issn.1001-7119.2000.06.014 [22] 陈晖, 刘坤松, 郭晓娟, 等. 珠江磨刀门河口表层沉积物磁性特征及其动力沉积环境意义[J]. 海洋学报, 2017, 39(3):44-54. [CHEN Hui, LIU Kunsong, GUO Xiaojuan, et al. Magnetic properties of surficial sediment and its implication for sedimentation dynamic environment in the Modaomen Outlet of the Pearl River Estuary [J]. Haiyang Xuebao, 2017, 39(3): 44-54. [23] 林炳尧, 黄世昌, 毛献忠, 等. 钱塘江河口潮波变化过程[J]. 水动力学研究与进展(A辑), 2002, 17(6):665-675. [LIN Bingyao, HUANG Shichang, MAO Xianzhong, et al. Deformation process of tidal waves in Qiantang Estuary [J]. Journal of Hydrodynamics, 2002, 17(6): 665-675. [24] 中国海湾志编纂委员会. 中国海湾志(第五分册)[M]. 北京: 海洋出版社, 1992. China Gulf Records Compilation Committee. Chinese Gulf Journal (Part 5)[M]. Beijing: China Ocean Press, 1992.
[25] Pejrup M. Flocculated suspended sediment in a micro-tidal environment [J]. Sedimentary Geology, 1988, 57(3-4): 249-256. doi: 10.1016/0037-0738(88)90032-2
[26] 周开胜, 孟翊, 刘苍字, 等. 长江口北支潮流沉积物磁性特征与沉积环境分析[J]. 海洋通报, 2008, 27(5):47-55. [ZHOU Kaisheng, MENG Yi, LIU Cangzi, et al. Magnetic properties and sedimentary environment of the Xinglong Sand in the North Branch, the Yangtze Estuary [J]. Marine Science Bulletin, 2008, 27(5): 47-55. doi: 10.3969/j.issn.1001-6392.2008.05.008 -
期刊类型引用(3)
1. 刘景良，苏杰龙，戴逸宸，李宇祖，黄永，郑文婷. 基于滑动窗宽优化的LMSSGST识别非平稳信号瞬时频率. 振动与冲击. 2024(19): 183-193 . 
百度学术
2. 李雪英，王鑫. 二阶同步提取变换的沉积旋回界面定位与追踪识别. 黑龙江科技大学学报. 2023(04): 614-621 . 百度学术
3. 吴耀乐，邓敏. 基于残差卷积神经网络的随机噪声压制网络的海洋地震资料去噪技术研究. 中国石油和化工标准与质量. 2022(22): 175-177 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 1535
- HTML全文浏览量: 268
- PDF下载量: 59
- 被引次数: 7
