CN107219551A - 拓宽地震数据频带的方法及装置 - Google Patents
拓宽地震数据频带的方法及装置 Download PDFInfo
- Publication number
- CN107219551A CN107219551A CN201610167841.0A CN201610167841A CN107219551A CN 107219551 A CN107219551 A CN 107219551A CN 201610167841 A CN201610167841 A CN 201610167841A CN 107219551 A CN107219551 A CN 107219551A
- Authority
- CN
- China
- Prior art keywords
- frequency
- seismic data
- mrow
- channel
- frequency domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims description 18
- 238000009499 grossing Methods 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 abstract description 21
- 238000004088 simulation Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种拓宽地震数据频带的方法及装置，所述方法包括以下步骤：步骤1：对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；步骤2：计算频率域单道地震数据的振幅；步骤3：对振幅进行拟合与光滑处理，获得拟合曲线，计算拟合曲线最大值步骤4：基于拟合曲线和拟合曲线最大值计算单道拓频滤波算子；步骤5：在频率域，基于单道拓频滤波算子计算经过拓频的频率域单道地震数据；步骤6：对经过拓频的频率域单道地震数据进行反傅里叶变换。本发明可以实现高频端和低频端的同时拓宽，通过拓频控制因子有效控制拓频，避免产生假象和高频噪声。此外，本发明以原始数据的频谱作为约束，可以最大程度的保持原有信号的主频特征，减少信号畸变。
Description
技术领域
本公开涉及地震勘探资料处理技术领域，特别涉及一种拓宽地震数据频带的方法及装置。
背景技术
在地震勘探中，地震数据的频谱分析是现代信号处理与分析中的重要手段之一。频谱分析是指将信号源发出的信号强度按频率顺序展开，使其成为频率的函数，并考察变化规律。将时域信号变换至频率域，其目的是把复杂的时间历程波形经过傅里叶变换分解为若干单一的谐波分量来研究，从而获得信号的频率结构以及各谐波和相位信息。
分辨率是分辨一个给定的地震剖面细节特征的能力。对于3D地震资料来说，地震资料分辨率包括纵向(时间上的)分辨率和横向(空间上的)分辨率。纵向分辨率主要取决于地震信号的频率，为了更好地分辨薄层，需要较宽的频带。由于大多数实际地震资料频带范围较窄，因此需要符合地质规律且有效的拓频方法对地震资料进行频带拓宽处理。常规处理方法中反褶积处理方法(如用于压缩子波的脉冲反褶积等)一般多用于拓展高频地震资料(有时也可以拓展低频地震资料)。总的来说，这些方法存在两个弊端：一是引入相当多的高频噪音；二是拓展的频宽一般均小于一个倍频程。同时拓宽地震信号的高频端和低频端的拓频技术一直是地震资料采集和处理所追求的目标。
发明内容
本公开的目的是提供一种拓宽地震数据频带的方法及装置，其针对地震资料频带范围较窄、分辨率较低的问题，基于单道地震数据，对地震数据进行拓频，能够同时拓宽地震信号的高频端和低频端，拓频后能有效提高地震数据分辨率，同时还能保持原始信号的基本特征。
本公开的一方面提供一种拓宽地震数据频带的方法，包括以下步骤：
步骤1：对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；
步骤2：计算所述频率域单道地震数据的振幅；
步骤3：对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；
步骤4：基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；
步骤5：在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；
步骤6：对经过拓频的频率域单道地震数据进行反傅里叶变换。
优选地，基于以下公式计算单道拓频滤波算子：
其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Ampmax'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。
优选地，ε的取值范围在0.01至1之间。
优选地，基于以下公式计算经过拓频的频率域单道地震数据：
F1(w)＝Fr1(w)+iFi1(w)
其中，
Fr1(w)＝Fr(w)×f(w)
Fi1(w)＝Fi(w)×f(w)
其中，F1(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr1(w)表示经过拓频的频率域单道地震数据的实部，Fi1(w)表示经过拓频的频率域单道地震数据的虚部。
本公开的另一方面提供一种拓宽地震数据频带的装置，包括：
傅里叶变换模块，用于对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；
频率域计算模块，用于计算所述频率域单道地震数据的振幅；
拟合模块，用于对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；
滤波算子确定模块，用于基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；
拓频计算模块，用于在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；
傅里叶反变换模块，用于对经过拓频的频率域单道地震数据进行反傅里叶变换。
优选地，所述滤波算子确定模块基于以下公式计算单道拓频滤波算子：
其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Ampmax'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。
优选地，ε的取值范围在0.01至1之间。
优选地，所述拓频计算模块基于以下公式计算经过拓频的频率域单道地震数据：
F1(w)＝Fr1(w)+iFi1(w)
其中，
Fr1(w)＝Fr(w)×f(w)
Fi1(w)＝Fi(w)×f(w)
其中，F1(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr1(w)表示经过拓频的频率域单道地震数据的实部，Fi1(w)表示经过拓频的频率域单道地震数据的虚部。
本公开的有益效果在于：(1)可以实现高频端和低频端的同时拓宽；(2)可以通过拓频控制因子有效控制拓频，避免产生假象和高频噪声；(3)由于以原始数据的频谱作为约束，可以最大程度的保持原有信号的主频特征，减少信号畸变。
附图说明
通过结合附图对本公开示例性实施例进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据示例性实施例的拓宽地震数据频带的方法的流程图；
图2(a)至图2(d)显示根据示例性实施例的拓宽地震数据频带的方法的应用示例，其中图2(a)显示原始雷克子波，图2(b)显示应用根据示例性实施例的拓宽地震数据频带的方法进行拓频处理后的子波，图2(c)显示原始雷克子波频谱，图2(d)显示拓频处理后的子波频谱；
图3(a)和图3(b)分别显示应用示例中的地震资料原始剖面和经过拓频处理后的剖面；
图4(a)和图4(b)分别示出图3(a)和图3(b)的波形显示；以及
图5(a)、图5(b)和图5(c)分别显示应用示例中的地震资料原始剖面频谱、经过拓频处理后的剖面频谱和处理前后剖面频谱叠合对比结果。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的优选实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
在信号处理领域，通过傅里叶变换和傅里叶反变换可以实现信号在时间域和频率域之间的转换。在频率域，傅里叶变换的典型用途是将信号分解成振幅谱——显示与频率对应的幅值大小，横轴一般为频率，纵轴为振幅。当原始地震信号低频化时，在频率域表现为子波主频小、频带较窄，因而具有较低的分辨率，将不利于信号的分析和进一步的地震解释工作。
针对上述问题，本公开提出一种拓宽地震数据频带的方法，其基于单道地震数据，将地震数据变换到频率域，并以该道地震数据的频谱特征作为约束，构建拓频滤波算子，然后将拓频滤波算子作用于原始数据频谱来实现原始单道地震数据的拓频处理。
图1示出了根据示例性实施例的拓宽地震数据频带的方法的流程图，其包括以下步骤：
步骤1：对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据
以f(t)表示时间域单道地震数据，以F(w)表示频率域单道地震数据，根据傅里叶变换，即有：
频率域单道地震数据F(w)也可表述为公式(2)：
F(w)＝Fr(w)+iFi(w) (2)
其中，Fr(w)为频率域单道地震数据F(w)的实部，Fi(w)为频率域单道地震数据F(w)的虚部。
步骤2：计算频率域单道地震数据的振幅
频率域单道地震数据F(w)的振幅表述为公式(3)，其相位可表述为公式(4)：
其中Amp(w)表示频率域单道地震数据F(w)的振幅，tan(θ)表示频率域单道地震数据F(w)的相位。
步骤3：对振幅进行拟合与平滑处理，获得拟合曲线，计算拟合曲线最大值
对振幅Amp(w)进行拟合处理，例如可采用基于误差最小的多项式拟合法进行拟合处理，然后进行平滑处理，得到平滑的拟合曲线Amp'(w)。通过进行拟合与平滑处理，可以得到振幅谱的振幅包络，便于求取最大值，消除局部极值的干扰。
通过公式(5)计算拟合曲线Amp'(w)的最大值Ampmax'(w)：
Ampmax'(w)＝arg max(Amp'(w)) (5)
步骤4：基于拟合曲线和拟合曲线最大值计算单道拓频滤波算子，作为单道滤波器
根据公式(6)计算单道拓频滤波算子：
其中，参数ε为拓频控制因子，取值范围在0.01至1之间，其值越小频谱的频带展宽越大，但频带展宽过大会引入假象和高频噪声，ε的具体取值可根据对实际数据进行参数测试确定。
步骤5：在频率域，基于单道拓频滤波算子计算经过拓频的频率域单道地震数据
在频率域，根据公式(7)和公式(8)将拓频滤波算子分别作用于频率域单道地震数据F(w)的实部和虚部：
Fr1(w)＝Fr(w)×f(w) (7)
Fi1(w)＝Fi(w)×f(w) (8)
其中，Fr1(w)表示经过拓频的频率域单道地震数据的实部，Fi1(w)表示经过拓频的频率域单道地震数据的虚部。
那么，经过拓频的频率域单道地震数据F1(w)可以表述为公式(9)：
F1(w)＝Fr1(w)+iFi1(w) (9)
步骤6：对经过拓频的频率域单道地震数据进行反傅里叶变换，得到经过拓频的时间域单道地震数据
经过拓频的时间域单道地震数据表述为公式(10)：
根据上述方法，可以对所有处理道进行步骤1到步骤6的处理，从而对所有道的地震数据进行频率拓宽。
本公开还提供一种拓宽地震数据频带的装置，根据示例性实施例的拓宽地震数据频带的装置包括：
傅里叶变换模块，用于对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；
频率域计算模块，用于计算频率域单道地震数据的振幅；
拟合模块，用于对振幅进行拟合与平滑处理，获得拟合曲线，计算拟合曲线最大值；
滤波算子确定模块，用于基于拟合曲线和拟合曲线最大值计算单道拓频滤波算子；
拓频计算模块，用于在频率域，基于单道拓频滤波算子计算经过拓频的频率域单道地震数据；
傅里叶反变换模块，用于对经过拓频的频率域单道地震数据进行反傅里叶变换。
应用示例
图2显示了根据示例性实施例的拓宽地震数据频带的方法的应用示例，其中图2(a)显示原始雷克子波，图2(b)显示应用根据示例性实施例的拓宽地震数据频带的方法进行拓频处理后的子波，图2(c)显示原始雷克子波频谱，图2(d)显示拓频处理后的子波频谱。
在该示例中，选取主频为30Hz的雷克子波，采样点数为101，采样率为4ms。在拓频处理中，选择拓频控制因子ε＝0.1。
对比图2(a)与图2(b)可以看出，经过拓频处理后，提高了子波的分辨率，同时子波旁瓣和原始子波旁瓣保持了很好的一致性，没有产生畸变和噪声干扰。对比图2(c)和图2(d)可以看出，处理后子波频谱保持了原始子波频谱特征，同时频带得到有效拓宽。
图3(a)和图3(b)分别显示应用示例中的地震资料原始剖面和经过拓频处理后的剖面，图4(a)和图4(b)分别示出图3(a)和图3(b)的波形显示。
对比图3(a)和图3(b)显示的剖面可以看出，经过拓频处理后剖面同相轴连续性增强，同相轴波形变“瘦”，表现为分辨率提高，处理后波形一致性问题得到一定程度的改善。图4中的波形显示更直观地体现了处理前后的效果对比。
图5(a)、图5(b)和图5(c)分别显示应用示例中的地震资料原始剖面频谱、经过拓频处理后的剖面频谱和处理前后剖面频谱叠合对比结果。
从图5(c)的频谱叠合对比结果可以看出，应用根据示例性实施例的拓宽地震数据频带的方法进行拓频处理后，保持了原始数据的主频以及频谱特征，数据频带得到很好的展宽，同时避免了产生高频噪声对数据质量的影响。
本公开的理论模型及实际地震数据试验结果表明，根据示例性实施例的拓宽地震数据频带的方法有效提高了地震数据的分辨率，保持了原始数据的基本频谱特征，同时避免了噪声干扰的产生和信噪比的降低，是一种符合地质规律且有效的地震数据拓频方法。
上述技术方案只是本发明的一种实施例，对于本领域内的技术人员而言，在本发明公开的原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施例的描述，因此前面的描述只是优选的，而并不具有限制性的意义。
Claims (8)
1.一种拓宽地震数据频带的方法，包括以下步骤：
步骤1：对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；
步骤2：计算所述频率域单道地震数据的振幅；
步骤3：对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；
步骤4：基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；
步骤5：在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；
步骤6：对经过拓频的频率域单道地震数据进行反傅里叶变换。
2.根据权利要求1所述的拓宽地震数据频带的方法，其中基于以下公式计算单道拓频滤波算子：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>Amp</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>Amp</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<mo>&times;</mo>
<msup>
<msub>
<mi>Amp</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Ampmax'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。
3.根据权利要求2所述的拓宽地震数据频带的方法，其中ε的取值范围在0.01至1之间。
4.根据权利要求1所述的拓宽地震数据频带的方法，其中基于以下公式计算经过拓频的频率域单道地震数据：
F1(w)＝Fr1(w)+iFi1(w)
其中，
Fr1(w)＝Fr(w)×f(w)
Fi1(w)＝Fi(w)×f(w)
其中，F1(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr1(w)表示经过拓频的频率域单道地震数据的实部，Fi1(w)表示经过拓频的频率域单道地震数据的虚部。
5.一种拓宽地震数据频带的装置，包括：
傅里叶变换模块，用于对时间域单道地震数据进行傅里叶变换，得到频率域单道地震数据；
频率域计算模块，用于计算所述频率域单道地震数据的振幅；
拟合模块，用于对所述振幅进行拟合与平滑处理，获得拟合曲线，计算所述拟合曲线最大值；
滤波算子确定模块，用于基于所述拟合曲线和拟合曲线最大值计算单道拓频滤波算子；
拓频计算模块，用于在频率域，基于所述单道拓频滤波算子计算经过拓频的频率域单道地震数据；
傅里叶反变换模块，用于对经过拓频的频率域单道地震数据进行反傅里叶变换。
6.根据权利要求5所述的拓宽地震数据频带的装置，其中所述滤波算子确定模块基于以下公式计算单道拓频滤波算子：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>Amp</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>Amp</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<mo>&times;</mo>
<msup>
<msub>
<mi>Amp</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中，f(w)表示所述单道拓频滤波算子，Amp'(w)表示所述拟合曲线，Ampmax'(w)表示所述拟合曲线最大值，ε表示拓频控制因子。
7.根据权利要求6所述的拓宽地震数据频带的装置，其中ε的取值范围在0.01至1之间。
8.根据权利要求1所述的拓宽地震数据频带的装置，其中所述拓频计算模块基于以下公式计算经过拓频的频率域单道地震数据：
F1(w)＝Fr1(w)+iFi1(w)
其中，
Fr1(w)＝Fr(w)×f(w)
Fi1(w)＝Fi(w)×f(w)
其中，F1(w)表示经过拓频的频率域单道地震数据，f(w)表示所述单道拓频滤波算子，Fr(w)表示所述频率域单道地震数据的实部，Fi(w)表示所述频率域单道地震数据的虚部，Fr1(w)表示经过拓频的频率域单道地震数据的实部，Fi1(w)表示经过拓频的频率域单道地震数据的虚部。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610167841.0A CN107219551B (zh) | 2016-03-22 | 2016-03-22 | 拓宽地震数据频带的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610167841.0A CN107219551B (zh) | 2016-03-22 | 2016-03-22 | 拓宽地震数据频带的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107219551A true CN107219551A (zh) | 2017-09-29 |
| CN107219551B CN107219551B (zh) | 2019-08-27 |
Family
ID=59928279
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610167841.0A Active CN107219551B (zh) | 2016-03-22 | 2016-03-22 | 拓宽地震数据频带的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107219551B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108363739A (zh) * | 2018-01-22 | 2018-08-03 | 中国海洋大学 | 一种基于稀疏采集的地震资料高低频拓展方法 |
| WO2019164606A1 (en) * | 2018-02-23 | 2019-08-29 | Exxonmobil Research And Engineering Company | Method of extracting intrinsic attenuation from seismic data |
| CN113721294A (zh) * | 2020-05-26 | 2021-11-30 | 中国石油化工股份有限公司 | 复数域最小二乘约束谱蓝化拓频方法 |
| CN115373035A (zh) * | 2022-08-15 | 2022-11-22 | 江苏南京地质工程勘察院 | 探地雷达低频信号的拓频方法及系统 |
| CN118465843A (zh) * | 2024-07-08 | 2024-08-09 | 中国地质科学院地球物理地球化学勘查研究所 | 保持相位特征的地震数据处理方法、装置、介质及产品 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176056A (zh) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | 基于时频分析的地震有效波多域能量补偿方法 |
| CN104122583A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种拓宽地震数据频谱的方法和装置 |
-
2016
- 2016-03-22 CN CN201610167841.0A patent/CN107219551B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176056A (zh) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | 基于时频分析的地震有效波多域能量补偿方法 |
| CN104122583A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种拓宽地震数据频谱的方法和装置 |
Non-Patent Citations (6)
| Title |
|---|
| BAI JUNYU ET AL.: "Nonlinear hybrid optimization algorithm for seismic impedance inversion", 《CPS/SEG BEIJING 2014 INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
| 张军华等: "地震低频信息缺失特征分析及拓频方法研究", 《APPLIED GEOPHYSICS（应用地球物理）》 * |
| 张在金等: "基于 GST 时变子波谱模拟提高分辨率方法研究", 《中国地球科学联合学术年会 2015》 * |
| 张盼等: "基于被动源多窗谱方法的主动源地震低频数据重构", 《APPLIED GEOPHYSICS（应用地球物理）》 * |
| 王小青等: "CPU/GPU协同计算在频率域二维全波形反演中的应用", 《CT 理论与应用研究》 * |
| 陈双全等: "应用傅里叶尺度变换提高地震资料分辨率", 《石油地球物理勘探》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108363739A (zh) * | 2018-01-22 | 2018-08-03 | 中国海洋大学 | 一种基于稀疏采集的地震资料高低频拓展方法 |
| CN108363739B (zh) * | 2018-01-22 | 2021-07-30 | 中国海洋大学 | 一种基于稀疏采集的地震资料高低频拓展方法 |
| WO2019164606A1 (en) * | 2018-02-23 | 2019-08-29 | Exxonmobil Research And Engineering Company | Method of extracting intrinsic attenuation from seismic data |
| CN113721294A (zh) * | 2020-05-26 | 2021-11-30 | 中国石油化工股份有限公司 | 复数域最小二乘约束谱蓝化拓频方法 |
| CN113721294B (zh) * | 2020-05-26 | 2023-06-23 | 中国石油化工股份有限公司 | 复数域最小二乘约束谱蓝化拓频方法 |
| CN115373035A (zh) * | 2022-08-15 | 2022-11-22 | 江苏南京地质工程勘察院 | 探地雷达低频信号的拓频方法及系统 |
| CN115373035B (zh) * | 2022-08-15 | 2024-08-23 | 江苏南京地质工程勘察院 | 探地雷达低频信号的拓频方法及系统 |
| CN118465843A (zh) * | 2024-07-08 | 2024-08-09 | 中国地质科学院地球物理地球化学勘查研究所 | 保持相位特征的地震数据处理方法、装置、介质及产品 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107219551B (zh) | 2019-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107219551B (zh) | 拓宽地震数据频带的方法及装置 | |
| CN106597532B (zh) | 一种结合井资料与层位资料的叠前地震数据频带拓展方法 | |
| DE69700084T2 (de) | Verfahren zur Transformierung eines periodischen Signales unter Verwendung eines geplätteten Spectrogrammes, Verfahren zur Transformierung von Schall bei Verwendung von Phasenkomponenten und Verfahren zur Analyse eines Signales unter Verwendung einer optimalen Interpolationsfunktion | |
| CN108009347B (zh) | 基于同步压缩联合改进广义s变换的时频分析方法 | |
| Czarnecki | The instantaneous frequency rate spectrogram | |
| CN110187388B (zh) | 一种基于变分模态分解的稳定地震品质因子q估计方法 | |
| CN106405645B (zh) | 一种基于资料品质分析的信噪比可控的地震拓频处理方法 | |
| CN110336587B (zh) | 一种多跳频信号侦察中获取组合时频分布的方法 | |
| CN105760347A (zh) | 一种基于数据/极值联合对称延拓的hht端点效应抑制方法 | |
| CN103728660A (zh) | 基于地震数据的多道匹配追踪方法 | |
| CN103258543B (zh) | 一种人工语音带宽扩展的方法 | |
| CN106483563A (zh) | 基于互补集合经验模态分解的地震能量补偿方法 | |
| CN106707334A (zh) | 一种提高地震资料分辨率的方法 | |
| CN106556865A (zh) | 一种串联型地震信号优化时频变换方法 | |
| Meena et al. | Cross-term suppression in the Wigner-Ville distribution using variational mode decomposition | |
| CN111198357A (zh) | 一种基于可调窗函数的s变换时频分析方法 | |
| CN113281809B (zh) | 一种地震信号的谱分析方法 | |
| CN107526103A (zh) | 地震资料处理方法及其阙值和有效信号频率的求取方法 | |
| CN111474582B (zh) | 生成高精度时频谱的精准s变换方法 | |
| CN109975867B (zh) | 一种频率域信号混叠的地震数据拓频方法 | |
| CN104658547A (zh) | 一种人工语音带宽扩展的方法 | |
| Ochiai et al. | A recursive adaptive method of impulse response measurement with constant SNR over target frequency band | |
| CN105093325B (zh) | 一种定量的提频方法 | |
| CN103885044A (zh) | 一种基于clean算法的窄带雷达回波杂噪抑制方法 | |
| CN109884705B (zh) | 双重约束时频域子波谱提高地震分辨率处理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |