CN106094038B - 适用于tti介质的频率域有限元全吸收pml方法 - Google Patents
适用于tti介质的频率域有限元全吸收pml方法 Download PDFInfo
- Publication number
- CN106094038B CN106094038B CN201610565960.1A CN201610565960A CN106094038B CN 106094038 B CN106094038 B CN 106094038B CN 201610565960 A CN201610565960 A CN 201610565960A CN 106094038 B CN106094038 B CN 106094038B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mover
- pml
- 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.)
- Active
Links
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种适用于TTI介质的频率域有限元全吸收PML方法，构造了一种频域有限元全吸收PML边界条件，实现了TTI介质复杂井眼条件的高精度声波测井模拟，得到了与理论解吻合的数值模拟结果，为声波测井仪器制造及利用声波测井资料评价复杂油气藏地层提供坚实的数据支撑；本发明综合了C‑PML当中的X、Z向吸收边界的拉伸因子和M‑PML当中的X、Z向吸收边界的拉伸因子分别融合为完全吸收的X向、Z向吸收边界的倏逝波抑制因子和切向衰减因子。本发明的有益效果是：通过建立全吸收PML的边界条件，兼具C‑PML和M‑PML的优势，实现了频域有限元PML吸收边界在求解TTI介质时的高效性。能高效、精确的进行TTI页岩地层的声波测井数值模拟，为声波测井仪器制造及页岩气藏开发提供必不可少的理论依据。
Description
技术领域
本发明属于地球物理(测井)勘探方法，尤其涉及弹性波数值模拟领域。
背景技术
随着油气田勘探开发的进行，页岩气等非常规油气资源成为勘探开发的热点。利用数值模拟手段考查页岩地层中井眼声波传播规律对有效勘探开发页岩油气及研究仪器构造指导仪器开发制造，具有重要意义。在各种常用数值模拟方法中，有限元方法具有其独特的优势：精度高、适应复杂边界、有利于不同场的耦合以及易于处理各种类型声源，在弹性波波场模拟中具有重要的应用价值。页岩地层可以简化为VTI(垂向横向各向同性)介质。直接对VTI地层进行模拟比较复杂，通过Bond变换旋转参考坐标系，使得VTI介质转化为TTI(水平横向各向同性)介质。频域有限元方法是TTI介质中有效的数值模拟方法。
井孔声场数值模拟所模拟的地层是无限大的，受计算机容量等的限制，数值模拟只能在有限区域内进行。为了能在有限区域中模拟无限区域中的波动过程，需要在有限区域的边界处引入吸收边界。现有比较先进的两种有限元吸收边界为C-PML和M-PML吸收边界，C-PML在极端入射角入射的情况下精度较高，但是在某些极端各向异性介质中会出现数值不稳定。M-PML在极端各向异性介质中具有较高的稳定性，但需要设置较厚的PML吸收层，大大增大了计算量，且修正因子不够优化的情况下，准确度不够高。
发明内容：
本发明的目的是提出适用于TTI介质的频率域有限元全吸收PML方法的技术方案。在频域下进行,具有更好的稳定性、更高的计算效率、更好的吸收效果、声源易于处理等优点；采用该吸收边界技术可以在吸收边界很薄的情况下有较高的精度和更好的稳定性，从而提高计算效率，此外在某些极端各向异性介质中，常规PML方法无法使用，本方法依然有效。
为了实现上述目的，本发明的技术方案是：适用于TTI介质的频率域有限元全吸收PML方法，用于构造频域有限元全吸收PML边界条件，实现TTI介质复杂井眼条件的高精度声波测井模拟；
所述方法采用复拉伸坐标系的弹性波方程，所述复拉伸坐标系的弹性波方程在所述频域下表示形式为：
在PML区域内，加入拉伸函数，将上述方程表示为：
所述PML是完全匹配层，sx和sz为PML吸收边界内的拉伸函数，sx和sz分别构造x和z方向的两个拉伸函数分量，sx和sz的表达式如下：
上述式中：
为x方向的实数衰减因子，
为x方向复频移因子，
为x方向倏逝波抑制因子，
Lx为吸收层x方向的厚度；
m为切向衰减因子；
为z方向的实数衰减因子，
为z方向复频移因子，
为z方向倏逝波抑制因子，
Lz为吸收层z方向的厚度，
ω为角频率，为频率的2π倍，
i为虚数单位。
更进一步，所述实数衰减因子复频移因子α0＝2ω，倏逝波抑制因子中k0取值范围为1～20，切向衰减因子m取值范围为0.005～0.02。
本发明的主要优点是：频率域有限元的全吸收PML，兼具有普通PML，C-PML和M-PML的优点，实现了TTI介质复杂井眼条件的高精度声波测井模拟，获得了在TTI介质中进行声波测井有限元模拟时的高效性、稳定性，可以设置非常薄的吸收边界，节省计算时间和计算内存需求，在极端各向异性介质中也可保持良好的稳定性；可以准确的实现TTI介质声波测井数值模拟，获得复杂介质大斜度井中的高精度模拟结果；得到了与理论解吻合的数值模拟结果，为声波测井仪器制造及利用声波测井资料评价复杂油气藏地层提供坚实的数据支撑。
下面结合附图和实施实例对本发明作详细描述。
附图说明
图1是本发明一个各向同性均匀模型分别采用从左到右M-PML、C-PML和全吸收PML的波场传播效果图；
图2是本发明一个极端各向异性TTI介质中m＝0.000和m＝0.02时的波场快照图，左侧为m＝0.000，右侧为m＝0.02时的波场快照图；
图3是均匀各向同性模型的参考解的波形图；
图4是本发明当m＝0.000时的波形图；
图5是本发明当m＝0.02(即采用全吸收PML时)，吸收效果较好时的波形图；
具体实施方式：
下面对本发明做进一步的详细说明。
频域有限元法具有更稳定、计算效率更高、PML效果更好以及更容易处理各种类型的声源等优点，在地球物理正演模拟尤其是在页岩气地层声波测井模拟中具有重要的应用价值。频域有限元法的基本思路是通过有限元方法求解频域下的弹性波方程，获得频谱并通过傅里叶反变换获得时域下的波形为了同时满足稳定性和准确性的要求，本发明构造了一种全吸收PML吸收边界，采用最优化的参数，使得新构造的PML兼具C-PML和M-PML两种吸收边界的优势。
本发明主要在二维直角坐标系下完成。在该坐标系下弹性波在弹性介质中的传播满足弹性波方程，在时间域的表达式为：
其中ρ为弹性介质的密度，为位移向量，C为广义虎克矩阵，在二维直角坐标系下的表达式为：
弹性波在频域下的形式满足弹性波方程
在二维直角坐标系下，上式化为
其中各应力变量的表达式为：
引入复拉伸坐标系
可将弹性波方程在频域下表示为如下形式：
进行坐标系复拉伸后的应力变量为
上式即为带有PML的弹性波方程在频域下的表达式。在非PML区域，sx和sz的取值为1；在PML区域，sx和sz即为本发明构造的PML完全匹配层的x方向和z方向的拉伸函数。sx和sz的表达式如下：
上述式中：
为x方向的实数衰减因子，
为x方向复频移因子，
为x方向倏逝波抑制因子，
Lx为吸收层x方向的厚度；
m为切向衰减因子；
为z方向的实数衰减因子，
为z方向复频移因子，
为z方向倏逝波抑制因子，
Lz为吸收层z方向的厚度，
ω为角频率，为频率的2π倍，
i为虚数单位。
其中，实数衰减因子复频移因子α0＝2ω，倏逝波抑制因子中k0取值范围为1～20，切向衰减因子m取值范围为0.005～0.02。
为了便于推导弹性波方程在频域下的等效积分弱形式，将各等式的左右两边均乘以变量sxsz，可得
将上式与势函数在弹性介质区域ΩE作内积，并对公式两边代入格林公式，整理可得
其中LE代表弹性介质区域的边界，代表弹性介质边界的外法向方向。
下面求解流固耦合问题。声波方程在频域下的等效积分弱形式
声波场和位移场在流固边界处在时间域满足如下公式：
其中ρf为流体介质的密度，在频率域的表达式为：
故上式最终化为：
其中LA→E代表流体区域和弹性介质区域的交界面。
由于流体区域和弹性介质区域的交界面需要满足法向应力连续，切向应力为零的边界条件，故
最终化为：
其中LE→A代表弹性介质区域和流体区域的交界面。
大量数据模拟实例表明，相比于传统PML和C-PML，优化的α和k能够在无精度损失的前提下极大的提高吸收效果。为了实现在TTI介质的模拟中足够精确，令m取一个足够小的值(通常为0.005到0.02之间)，就可以在不损失精度的前提下达到最大的稳定性。
实施例一：
考虑一个TTI声波测井模型，其声源采用的是中心频率为3kHz的单极子声源。井外地层为将一个VTI地层旋转45°所得，其广义Hooke矩阵为在此问题中，吸收边界层只占8个网格的厚度，远远小于波长，这对数值模拟方法是一个非常大的挑战。
图1从左到右分别为采用M-PML，C-PML和本发明的完全吸收PML吸收边界的波场图。图中可以看到，针对该模型，M-PML和C-PML吸收边界下，随着波场的传播有非常明显的边界反射和不稳定，图中虚线所示区域，而采用本发明提出的完全吸收PML边界条件，有更好的吸收效果和更好的稳定性。
接下来考察m取值不同时的吸收效果对比。图2为前文所述TTI介质模型中m取值不同时的频域响应。左图为m＝0.000时，3000Hz频率的响应特征，在这种情况下，在左上角可观察到明显的虚假反射；而在右图m＝0.02的情况下，可看到虚假反射得到了很好的消除，具有更好的吸收效果。
另外考查接收器接收到的波形特征。其实现方法为以100Hz为间隔，共取100个频率点，得到从100Hz到10000Hz的频率响应特征，并通过傅里叶反变换计算得到时域波形。图3为时域理论参考解的波形。图4为当m＝0.000时的时域波形特征，此时可观察到全场存在明显的虚假反射，在6ms前后，八个接收器上都有不同的来自界面的反射信号；而在图5当m＝0.02时的时域波形特征中，全场无虚假反射，且与参考解吻合较好。说明该全吸收PML吸收边界的参数优化对于数值模拟的准确性是非常重要的。
本发明的方法将已被用于一个实际的页岩TTI地层大斜度井问题中，并取得了良好的效果。
Claims (2)
1.适用于TTI介质的频率域有限元全吸收PML方法，用于构造频域有限元全吸收PML边界条件，实现TTI介质复杂井眼条件的高精度声波测井模拟；其特征在于：
所述方法采用复拉伸坐标系的弹性波方程，所述复拉伸坐标系的弹性波方程在所述频域下表示形式为：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>&rho;</mi>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>&rho;</mi>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
在PML区域内，加入拉伸函数，将上述方程表示为：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>&rho;</mi>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>s</mi>
<mi>x</mi>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>s</mi>
<mi>z</mi>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>&rho;</mi>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mover>
<mi>w</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>s</mi>
<mi>x</mi>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>s</mi>
<mi>z</mi>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
所述PML是完全匹配层，sx和sz为PML吸收边界内的拉伸函数，sx和sz分别构造x和z方向的两个拉伸函数分量，sx和sz的表达式如下：
<mrow>
<msub>
<mi>s</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<msub>
<mi>md</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<mi>i</mi>
<mi>&omega;</mi>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>d</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<msub>
<mi>md</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<mi>i</mi>
<mi>&omega;</mi>
</mrow>
</mfrac>
</mrow>
上述式中：
为x方向的实数衰减因子，
为x方向复频移因子，
为x方向倏逝波抑制因子，
Lx为吸收层x方向的厚度；
m为切向衰减因子；
为z方向的实数衰减因子，
为z方向复频移因子，
为z方向倏逝波抑制因子，
Lz为吸收层z方向的厚度，
ω为角频率，为频率的2π倍，
i为虚数单位。
2.根据权利要求1所述的适用于TTI介质的频率域有限元全吸收PML方法，其特征在于，所述实数衰减因子复频移因子α0＝2ω，倏逝波抑制因子中k0取值范围为1～20，切向衰减因子m取值范围为0.005～0.02。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610565960.1A CN106094038B (zh) | 2016-07-18 | 2016-07-18 | 适用于tti介质的频率域有限元全吸收pml方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610565960.1A CN106094038B (zh) | 2016-07-18 | 2016-07-18 | 适用于tti介质的频率域有限元全吸收pml方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106094038A CN106094038A (zh) | 2016-11-09 |
| CN106094038B true CN106094038B (zh) | 2017-11-14 |
Family
ID=57220655
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610565960.1A Active CN106094038B (zh) | 2016-07-18 | 2016-07-18 | 适用于tti介质的频率域有限元全吸收pml方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094038B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108229000B (zh) * | 2017-12-29 | 2020-06-09 | 电子科技大学 | 利用混合的三棱柱—四面体网格实现dgtd中pml的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
| CN102722651A (zh) * | 2012-06-01 | 2012-10-10 | 西安理工大学 | 二维柱坐标完全匹配吸收边界的实现方法 |
| CN103616721A (zh) * | 2013-11-25 | 2014-03-05 | 中国石油天然气股份有限公司 | 基于二阶偏微分波动方程的pml吸收边界条件的方法 |
| CN104237944B (zh) * | 2014-10-09 | 2015-12-30 | 王兵 | 一种适用于交错网格有限差分的全吸收pml方法 |
-
2016
- 2016-07-18 CN CN201610565960.1A patent/CN106094038B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
| CN102722651A (zh) * | 2012-06-01 | 2012-10-10 | 西安理工大学 | 二维柱坐标完全匹配吸收边界的实现方法 |
| CN103616721A (zh) * | 2013-11-25 | 2014-03-05 | 中国石油天然气股份有限公司 | 基于二阶偏微分波动方程的pml吸收边界条件的方法 |
| CN104237944B (zh) * | 2014-10-09 | 2015-12-30 | 王兵 | 一种适用于交错网格有限差分的全吸收pml方法 |
Non-Patent Citations (1)
| Title |
|---|
| 多轴卷积完全匹配层吸收边界条件;田坤;《石油地球物理勘探》;20140228;第49卷(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094038A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101251604B (zh) | 二参数转换波速度分析及动校正方法 | |
| CN104570072B (zh) | 一种粘弹性介质中的球面pp波反射系数建模方法 | |
| CN105095631A (zh) | 一种页岩各向异性岩石物理建模方法 | |
| CN103149586B (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
| CN106019375B (zh) | 一种页岩气地层层理地球物理评价方法 | |
| CN102445708A (zh) | 三维等效富泥质砂岩速度预测模型 | |
| CN101201409B (zh) | 一种地震数据变相位校正方法 | |
| CN102176053A (zh) | 提升波动方程叠前深度偏移成像效果的方法 | |
| CN112987088B (zh) | 一种渗流介质地震横波数值模拟和成像方法 | |
| CN104749628A (zh) | 一种基于弥散黏滞性波动方程的吸收边界反射方法 | |
| CN109655890B (zh) | 一种深度域浅中深层联合层析反演速度建模方法及系统 | |
| CN102707313B (zh) | 一种基于脉冲耦合神经网络的拟声波曲线构建方法 | |
| CN102262243A (zh) | 一种滤波法可控震源地震数据谐波干扰压制方法 | |
| CN105954803A (zh) | 叠后地震反演方法及装置 | |
| CN106094038B (zh) | 适用于tti介质的频率域有限元全吸收pml方法 | |
| Gao et al. | Unsplit complex frequency shifted perfectly matched layer for second-order wave equation using auxiliary differential equations | |
| CN104237944B (zh) | 一种适用于交错网格有限差分的全吸收pml方法 | |
| CN106443773B (zh) | 一种压制地震剖面数据中煤层屏蔽效应的方法及装置 | |
| CN105447225A (zh) | 一种应用于声波有限差分数值模拟的组合吸收边界条件 | |
| CN105527648B (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及系统 | |
| CN104007465B (zh) | 基于弹性反演的纯纵波数据刻画储层发育细节方法 | |
| CN107315193B (zh) | 一种利用地震瞬时质心频率进行含气性检测的方法 | |
| CN102269821A (zh) | 一种wefox分裂法双向聚焦叠前地震成像方法 | |
| CN109358363A (zh) | 基于多次反射等效界面模型的三维鸣震追踪与预测方法 | |
| CN107024716B (zh) | 一种地震波场吸收补偿成像方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170620 Address after: 102249 Beijing city Changping District Road No. 18 Applicant after: China University of Petroleum (Beijing) Address before: College of information China University of Petroleum No. 18 Beijing city Changping District Road 102249 Applicant before: Wang Bing |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |