CN106772590B - 一种剧烈起伏自由地表有限差分正演模拟系统及方法 - Google Patents
一种剧烈起伏自由地表有限差分正演模拟系统及方法 Download PDFInfo
- Publication number
- CN106772590B CN106772590B CN201710160189.4A CN201710160189A CN106772590B CN 106772590 B CN106772590 B CN 106772590B CN 201710160189 A CN201710160189 A CN 201710160189A CN 106772590 B CN106772590 B CN 106772590B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- wave field
- equation
- under
- module
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000009795 derivation Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000000844 transformation Methods 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
本发明公开了一种剧烈起伏自由地表有限差分正演模拟系统及方法,属于石油勘探,本发明能够对剧烈起伏地表进行准确模拟的有限差分正演模拟,通过引入过渡坐标系和矩形网格坐标系,在三个坐标系之间进行转换,既能克服传统传统有限差分方法在处理起伏地表时的缺陷,又能保持有限差分方法的精度,而不明显增加计算量,开发剧烈起伏自由地表有限差分正演模拟方法,为研究剧烈起伏地表地区的波场传播规律提供了理论基础,同时,为剧烈起伏地表地区的采集数据进行后续成像和反演工作提供了高精度的地震波正演基础。
Description
技术领域
本发明属于石油勘探领域,具体涉及一种剧烈起伏自由地表有限差分正演模拟系统及方法。
背景技术
剧烈的起伏地表广泛地存在于中国地区。地震勘探的前缘正逐步向复杂地表地区转移,常规的处理和解释方法对上述复杂起伏地表已经不再适用。为了更好地解决勘探开发难题,须密切研究起伏地表条件下介质地震波正演模拟,对于起伏地表条件下的地质构造进行正演模拟,可从数值模拟结果中认识、分析地震波波场特征及地震波传播规律等,更好地指导后续的地震勘探。
有限差分方法因其计算速度快、占用内存低、容易实现等优点得到广泛应用,但有限差分方法是基于矩形网格剖分的,因此在处理起伏地表时存在很大的缺陷,在模拟地震波传播时会产生明显的绕射噪音及不稳定。一种常用的解决方法是对物理模型进行曲网格剖分,并将物理域的曲网格变换为矩形网格,但该方法也存在一定的局限性,该方法对稳定性的要求较高,无法处理剧烈的起伏地表。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种剧烈起伏自由地表有限差分正演模拟系统及方法,设计合理,克服了现有技术的不足,具有良好的效果。
为了实现上述目的,本发明采用如下技术方案:
一种剧烈起伏自由地表有限差分正演模拟系统,包括模型输入模块、正演模拟模块、速度场转换模块、波场计算模块、波场存储模块、波场反变换模块以及结果输出模块;
模型输入模块,被配置为用于输入速度场和起伏地表的高程;
正演模拟模块,被配置为用于进行剧烈起伏自由地表有限差分正演模拟;
速度场转换模块,被配置为用于将速度场变换到矩形网格坐标系下;
波场计算模块,被配置为用于在过渡坐标系下更新地下区域波场以及在过渡坐标系下采用自由地表边界条件更新地表波场;
波场存储模块,被配置为用于在矩形网格坐标系下存储波场;
波场反变换模块,被配置为用于将波场变换到笛卡尔坐标系下;
结果输出模块,被配置为用于输出最终的波场快照和炮记录。
此外,本发明还提到一种剧烈起伏自由地表有限差分正演模拟方法,该方法采用如上所示的一种剧烈起伏自由地表有限差分正演模拟系统,包括如下步骤:
步骤1:通过模型输入模块输入速度场和起伏地表的高程,建立观测系统;
步骤2:通过速度场转换模块采用如下所示的变换方程将速度场变换到矩形网格坐标系下;
其中,x和z为笛卡尔坐标下的空间坐标,ξ和η为矩形网格坐标系下的空间坐标,z0(ξ)为笛卡尔坐标系下的高程函数,ηmax为矩形网格坐标系下的最大采样点数;
步骤3:引入过渡坐标系,通过波场计算模块在过渡坐标系下采用如下所示的一阶速度-应力弹性波方程更新地下区域波场;
其中,ux和uz分别为水平分量和垂直分量的速度,τxx和τzz是正应力,τxz是切应力,t是时间;λ和μ是拉梅常数,ρ是密度;表示过渡坐标系的横纵坐标;系数A和B由下式求得:
步骤4:通过波场计算模块在过渡坐标系下采用如下所示的自由地表边界条件更新地表波场;
其中,C-J八个系数可由下式求得:
步骤5:通过波场存储模块在矩形网格坐标系下存储区域波场和地表波场;
步骤6:通过波场反变换模块将区域波场和地表波场变换到笛卡尔坐标系下;
步骤7:通过结果输出模块输出最终的波场快照和炮记录。
优选地,对步骤3中过渡坐标系下的一阶速度-应力弹性波方程的推导过程如下:笛卡尔坐标系下的一阶速度-应力弹性波方程为:
应用链式法则,(2)式可变换为
对于方程(3),我们需要计算和可由下式求得
方程(4)变换得到方程(5)
由映射方程
可得
将方程(7)代入方程(5)可得
将方程(8)带入方程(3),得到过渡坐标系下的一阶速度-应力方程如方程(9)所示。
优选地,对步骤4中过渡坐标系下的自由地表边界条件的推导过程如下:
笛卡尔坐标系下的自由地表边界条件为:
其中
表示过渡坐标系的横坐标,将方程(11)代入到方程(10)中可得
对方程(12)应用链式法则可得
l表示过渡坐标系的纵坐标,在起伏地表处
将方程(14)代入方程(13)可得
整理方程(15)可得
本发明所带来的有益技术效果:
本发明能够对剧烈起伏地表进行准确模拟的有限差分正演模拟,通过引入过渡坐标系和矩形网格坐标系,在三个坐标系之间进行转换,既能克服传统传统有限差分方法在处理起伏地表时的缺陷,又能保持有限差分方法的精度,而不明显增加计算量,开发剧烈起伏自由地表有限差分正演模拟方法,为研究剧烈起伏地表地区的波场传播规律提供了理论基础,同时,为剧烈起伏地表地区的采集数据进行后续成像和反演工作提供了高精度的地震波正演基础。
附图说明
图1为本发明一种剧烈起伏自由地表有限差分正演模拟方法的流程图。
图2为三个坐标系网格示意图,其中,图(a)为笛卡尔坐系网格示意图;图(b)为过渡坐标系网格示意图;图(c)为矩形网格坐标系网格示意图。
图3为加拿大逆掩断层速度模型示意图,其中,图(a)为纵波速度模型示意图;图(b)为横波速度模型示意图。
图4为起伏地表高程函数示意图。
图5为加拿大逆掩断层网格剖分图,其中,图(a)为笛卡尔坐标系网格剖分图;图(b)为过渡坐标系网格剖分图。
图6为采用本发明方法得到的波场快照图,其中,图(a)为0.625s水平分量的波场快照图;图(b)为0.625s垂直分量的波场快照图;图(c)为0.7s水平分量的波场快照图;图(d)为0.7s垂直分量的波场快照图。
图7为采用传统方法得到的波场快照图,其中,图(a)为0.625s水平分量的波场快照图;图(b)为0.625s垂直分量的波场快照图;图(c)为0.7s水平分的波场快照图量;图(d)为0.7s垂直分量的波场快照图。
图8为采用本发明方法得到的炮记录。
图9为本发明一种剧烈起伏自由地表有限差分正演模拟系统的结构示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
实施例1
如图9所示,一种剧烈起伏自由地表有限差分正演模拟系统,包括模型输入模块、正演模拟模块、速度场转换模块、波场计算模块、波场存储模块、波场反变换模块以及结果输出模块;
模型输入模块,被配置为用于输入速度场和起伏地表的高程;
正演模拟模块,被配置为用于进行剧烈起伏自由地表有限差分正演模拟;
速度场转换模块,被配置为用于将速度场变换到矩形网格坐标系下;
波场计算模块,被配置为用于在过渡坐标系下更新地下区域波场以及在过渡坐标系下采用自由地表边界条件更新地表波场;
波场存储模块,被配置为用于在矩形网格坐标系下存储波场;
波场反变换模块,被配置为用于将波场变换到笛卡尔坐标系下;
结果输出模块,被配置为用于输出最终的波场快照和炮记录。
实施例2
在上述实施例的基础上,本发明还提到一种剧烈起伏自由地表有限差分正演模拟方法,其流程图如图1所示,包括如下步骤:
步骤1:通过模型输入模块输入速度场和起伏地表的高程,建立观测系统;
步骤2:通过速度场转换模块采用如下所示的变换方程将速度场变换到矩形网格坐标系下;
其中,x和z为笛卡尔坐标下的空间坐标,ξ和η为矩形网格坐标系下的空间坐标,z0(ξ)为笛卡尔坐标系下的高程函数,ηmax为矩形网格坐标系下的最大采样点数;
步骤3:引入过渡坐标系,通过波场计算模块在过渡坐标系下采用如下所示的一阶速度-应力弹性波方程更新地下区域波场;
其中,ux和uz分别为水平分量和垂直分量的速度,τxx和τzz是正应力,τxz是切应力,t是时间;λ和μ是拉梅常数,ρ是密度;表示过渡坐标系的横纵坐标;系数A和B由下式求得:
步骤4:通过波场计算模块在过渡坐标系下采用如下所示的自由地表边界条件更新地表波场;
其中,C-J八个系数可由下式求得:
步骤5:通过波场存储模块在矩形网格坐标系下存储区域波场和地表波场;
步骤6:通过波场反变换模块将区域波场和地表波场变换到笛卡尔坐标系下;
步骤7:通过结果输出模块输出最终的波场快照和炮记录。
优选地,对步骤3中过渡坐标系下的一阶速度-应力弹性波方程的推导过程如下:笛卡尔坐标系下的一阶速度-应力弹性波方程为:
应用链式法则,(2)式可变换为
对于方程(3),我们需要计算和可由下式求得
方程(4)变换得到方程(5)
由映射方程
可得
将方程(7)代入方程(5)可得
将方程(8)带入方程(3),得到过渡坐标系下的一阶速度-应力方程如方程(9)所示。
优选地,对步骤4中过渡坐标系下的自由地表边界条件的推导过程如下:
笛卡尔坐标系下的自由地表边界条件为:
其中
表示过渡坐标系的横坐标,将方程(11)代入到方程(10)中可得
对方程(12)应用链式法则可得
l表示过渡坐标系的纵坐标,在起伏地表处
将方程(14)代入方程(13)可得
整理方程(15)可得
图2为三个坐标系网格示意图,图(a)笛卡尔坐系网格示意图;图(b)过渡坐标系网格示意图;图(c)矩形网格坐标系网格示意图;
本发明一种剧烈起伏自由地表有限差分正演模拟方法,应用于国际标准的起伏地表加拿大逆掩断层模型(如图3所示),取得了理想的计算效果。输入速度场和起伏地表的高程(如图4所示),建立观测系统;将速度场变换到矩形网格坐标系下(如图5b所示);引入过渡坐标系,在过渡坐标系下(如图5a所示)更新地下区域波场;在过渡坐标系下采用自由地表边界条件更新地表波场(如图6所示);在矩形网格坐标系下存储波场;将波场变换到笛卡尔坐标系下;输出最终的波场快照和炮记录(如图8所示)。本发明的剧烈起伏自由地表有限差分正演模拟方法所得的波场快照(如图6所示),波形清晰,没有出现绕射噪音和不稳定现象。而传统有限差分方法得到的波场快照(如图7所示)出现了明显的不稳定。
本发明能够对剧烈起伏地表进行准确模拟的有限差分正演模拟,通过引入过渡坐标系和矩形网格坐标系,在三个坐标系之间进行转换,既能克服传统传统有限差分方法在处理起伏地表时的缺陷,又能保持有限差分方法的精度,而不明显增加计算量,开发剧烈起伏自由地表有限差分正演模拟方法,为研究剧烈起伏地表地区的波场传播规律提供了理论基础,同时,为剧烈起伏地表地区的采集数据进行后续成像和反演工作提供了高精度的地震波正演基础。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (4)
1.一种剧烈起伏自由地表有限差分正演模拟系统,其特征在于,包括模型输入模块、正演模拟模块、速度场转换模块、波场计算模块、波场存储模块、波场反变换模块以及结果输出模块;
模型输入模块,被配置为用于输入速度场和起伏地表的高程;
正演模拟模块,被配置为用于进行剧烈起伏自由地表有限差分正演模拟;
速度场转换模块,被配置为用于将速度场变换到矩形网格坐标系下;
波场计算模块,被配置为用于在过渡坐标系下更新地下区域波场以及在过渡坐标系下采用自由地表边界条件更新地表波场;
波场存储模块,被配置为用于在矩形网格坐标系下存储波场;
波场反变换模块,被配置为用于将波场变换到笛卡尔坐标系下;
结果输出模块,被配置为用于输出最终的波场快照和炮记录。
2.一种剧烈起伏自由地表有限差分正演模拟方法,其特征在于,采用如权利要求1所述的一种剧烈起伏自由地表有限差分正演模拟系统,包括如下步骤:
步骤1:通过模型输入模块输入速度场和起伏地表的高程,建立观测系统;
步骤2:通过速度场转换模块采用如下所示的变换方程将速度场变换到矩形网格坐标系下;
其中,x和z为笛卡尔坐标下的空间坐标,ξ和η为矩形网格坐标系下的空间坐标,z0(ξ)为笛卡尔坐标系下的高程函数,ηmax为矩形网格坐标系下的最大采样点数;
步骤3:引入过渡坐标系,通过波场计算模块在过渡坐标系下采用如下所示的一阶速度-应力弹性波方程更新地下区域波场;
其中,ux和uz分别为水平分量和垂直分量的速度,τxx和τzz是正应力,τxz是切应力,t是时间;λ和μ是拉梅常数,ρ是密度;表示过渡坐标系的横纵坐标;系数A和B由下式求得:
步骤4:通过波场计算模块在过渡坐标系下采用如下所示的自由地表边界条件更新地表波场;
其中,C-J八个系数可由下式求得:
步骤5:通过波场存储模块在矩形网格坐标系下存储区域波场和地表波场;
步骤6:通过波场反变换模块将区域波场和地表波场变换到笛卡尔坐标系下;
步骤7:通过结果输出模块输出最终的波场快照和炮记录。
3.根据权利要求2所述的剧烈起伏自由地表有限差分正演模拟方法,其特征在于,对步骤3中过渡坐标系下的一阶速度-应力弹性波方程的推导过程如下:
笛卡尔坐标系下的一阶速度-应力弹性波方程为:
应用链式法则,(2)式可变换为
对于方程(3),我们需要计算和可由下式求得
方程(4)变换得到方程(5)
由映射方程
可得
将方程(7)代入方程(5)可得
将方程(8)带入方程(3),得到过渡坐标系下的一阶速度-应力方程如方程(9)所示。
4.根据权利要求2所述的剧烈起伏自由地表有限差分正演模拟方法,其特征在于,对步骤4中过渡坐标系下的自由地表边界条件的推导过程如下:
笛卡尔坐标系下的自由地表边界条件为:
其中
表示过渡坐标系的横坐标,将方程(11)代入到方程(10)中可得
对方程(12)应用链式法则可得
表示过渡坐标系的纵坐标,在起伏地表处
将方程(14)代入方程(13)可得
整理方程(15)可得
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710160189.4A CN106772590B (zh) | 2017-03-17 | 2017-03-17 | 一种剧烈起伏自由地表有限差分正演模拟系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710160189.4A CN106772590B (zh) | 2017-03-17 | 2017-03-17 | 一种剧烈起伏自由地表有限差分正演模拟系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772590A CN106772590A (zh) | 2017-05-31 |
CN106772590B true CN106772590B (zh) | 2018-10-12 |
Family
ID=58966907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710160189.4A Active CN106772590B (zh) | 2017-03-17 | 2017-03-17 | 一种剧烈起伏自由地表有限差分正演模拟系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772590B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188512B (zh) * | 2018-09-17 | 2020-01-14 | 中国石油大学(华东) | 基于非规则扇形网格的起伏隧道空间正演模拟系统及方法 |
CN112684499B (zh) * | 2020-12-03 | 2022-12-06 | 南方科技大学 | 一种地震波仿真分析方法、存储介质及设备 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100545679C (zh) * | 2007-07-11 | 2009-09-30 | 成都理工大学 | 起伏地表的地震叠后正演方法 |
CN102062875A (zh) * | 2010-11-30 | 2011-05-18 | 中国石油集团川庆钻探工程有限公司 | 起伏地表弹性波波动方程正演方法 |
CN102565856B (zh) * | 2010-12-29 | 2013-11-13 | 中国石油天然气集团公司 | 一种基于波动方程正演的近地表噪音压制方法 |
CN102749643B (zh) * | 2011-04-22 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种面波地震记录的频散响应获取方法及其装置 |
CN102353988B (zh) * | 2011-07-08 | 2013-07-24 | 中国科学院地质与地球物理研究所 | 基于图形处理器计算起伏地表直接叠前逆时偏移的方法 |
CN102590857A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 真地表起伏叠前深度域双程波成像方法 |
CN104123449B (zh) * | 2014-07-16 | 2017-02-15 | 吉林大学 | 复杂山地区域的分区局部变加密不等距双重网格剖分方法 |
WO2016108896A1 (en) * | 2014-12-31 | 2016-07-07 | Landmark Graphics Corporation | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid |
-
2017
- 2017-03-17 CN CN201710160189.4A patent/CN106772590B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106772590A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104122585B (zh) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 | |
Kirchner et al. | Capabilities and limitations of numerical ice sheet models: a discussion for Earth-scientists and modelers | |
CN108646293B (zh) | 基于黏声拟微分方程的黏声起伏地表正演模拟系统及方法 | |
CN108508482A (zh) | 一种地下裂缝地震散射响应特征模拟方法 | |
CN105911584B (zh) | 一种隐式交错网格有限差分弹性波数值模拟方法及装置 | |
CN106772590B (zh) | 一种剧烈起伏自由地表有限差分正演模拟系统及方法 | |
Tayfun et al. | Distribution of sea-surface elevations in intermediate and shallow water depths | |
CN110244355B (zh) | 一种基于震源断层模型的脉冲地震动模拟方法 | |
CN106227957A (zh) | 等效裂缝建模的方法 | |
CN112464520B (zh) | 局部重力异常深度反演方法及装置 | |
CN102062875A (zh) | 起伏地表弹性波波动方程正演方法 | |
CN108647457B (zh) | 一种sccrf生成方法及系统 | |
CN111123374A (zh) | 一种基于匹配滤波的探地雷达全波形反演方法 | |
CN109459789A (zh) | 基于振幅衰减与线性插值的时间域全波形反演方法 | |
Long et al. | Data assimilation into a numerical equatorial ocean model. I. The model and the assimilation algorithm | |
Keylock | Hypothesis testing for nonlinear phenomena in the geosciences using synthetic, surrogate data | |
CN105807317B (zh) | 基于切比雪夫伪谱法的各向异性衰减面波模拟方法 | |
CN111399074A (zh) | 一种重力和重力梯度模量联合三维反演方法 | |
Krylov et al. | Numerical modeling of nonlinear response of seafloor porous saturated soil deposits to SH-wave propagation | |
CN113627657A (zh) | 一种使用机器学习模型的砂岩型铀成矿有利区预测方法 | |
CN105093265A (zh) | 一种模拟地震波在ti介质中传播规律的方法 | |
Kuznetsov et al. | Features of wind waves at the southeastern coast of Sakhalin according to bottom-pressure measurements | |
Zhou et al. | Numerical simulation of sea surface directional wave spectra under typhoon wind forcing | |
CN111126439A (zh) | 一种煤矿井下二维矿震成像训练数据集生成方法 | |
CN105319595A (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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 065000 No. 84 Jinguang Road, Guangyang District, Langfang City, Hebei Province Co-patentee after: China Petroleum University (East China) Patentee after: Institute of Geophysical and Geochemical Exploration under China Academy of Geos Address before: 065000 No. 84 Jinguangdao, Guangyang District, Handan City, Hebei Province Co-patentee before: China Petroleum University (East China) Patentee before: Institute of Geophysical and Geochemical Exploration under China Academy of Geos |