CN111856567B - 基于地震叠加速度的压实系数确定方法及装置 - Google Patents
基于地震叠加速度的压实系数确定方法及装置 Download PDFInfo
- Publication number
- CN111856567B CN111856567B CN201910363616.8A CN201910363616A CN111856567B CN 111856567 B CN111856567 B CN 111856567B CN 201910363616 A CN201910363616 A CN 201910363616A CN 111856567 B CN111856567 B CN 111856567B
- Authority
- CN
- China
- Prior art keywords
- well
- logging
- compaction
- speed
- seismic
- 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
- 238000005056 compaction Methods 0.000 title claims abstract description 186
- 238000000034 method Methods 0.000 title claims abstract description 112
- 238000004364 calculation method Methods 0.000 claims abstract description 106
- 230000008859 change Effects 0.000 claims abstract description 21
- 238000004590 computer program Methods 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 14
- 238000009933 burial Methods 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 20
- 230000033001 locomotion Effects 0.000 abstract description 8
- 239000007787 solid Substances 0.000 abstract description 7
- 230000008569 process Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 8
- 238000012937 correction Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 6
- 230000007547 defect Effects 0.000 description 6
- 238000011084 recovery Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 208000019888 Circadian rhythm sleep disease Diseases 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 208000001456 Jet Lag Syndrome Diseases 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000033915 jet lag type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition 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
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
技术领域
本发明涉及石油勘探领域,尤其是地震资料处理解释范畴,具体涉及一种基于地震叠加速度的地层压实系数确定方法及装置。
背景技术
地层压实系数是指地层岩石随着地层埋深深度的增加,岩石体积和孔隙受压减小的表征参数,是地层岩石弹性性质的重要参数,被广泛应用于油气藏动态地质储量评价中。
同一沉积相带地层,由于构造运动褶皱变形,地层埋藏差异很大,地层随埋藏深度变化,压实系数也随之越来越大,为了恢复沉积时地层厚度,需要求取地层压实系数,传统求取地层压实校正系数的方法是通过井点拟合求取平面上地层压实系数,但是,在无井区域或井少地区(例如勘探地区)地层压实校正系数只能通过井点外推的方法,该方法与使用外推插值方法有很大关系,而且存在一定的主关系,即不同的人可能得到不同压实校正系数平面变化规律。综上,目前缺乏一种行之有效的针对勘探地区确定压实系数的方法。
发明内容
针对现有技术中的问题,本发明能够建立一种在无井区域可以准确求取压实系数的方法。
为解决上述技术问题,本发明提供以下技术方案:
第一方面,本发明提供一种基于地震叠加速度的压实系数确定方法,包括:
根据地震叠加速度及待测区域有井地区的第一测井数据,计算待测区域无井地区的第二测井数据;
利用第二测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
一实施例中,基于地震叠加速度的压实系数确定方法还包括:
根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
一实施例中,根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型,包括:
根据多个第一测井数据及相应的深度数据,计算H-V速度计算公式的相变因子系数及压实系数;
利用H-V速度计算公式,根据相变因子系数及压实系数,生成压实系数计算模型。
一实施例中,根据地震叠加速度及待测区域有井地区的第一测井数据,计算待测区域无井地区的第二测井数据,包括:
将地震叠加速度及相应的第一测井数据进行拟合,获得测井数据计算模型;
利用测井数据计算模型,根据地震叠加速度计算待测区域无井地区的第二测井数据。
第二方面,本发明提供一种基于地震叠加速度的压实系数确定装置,该装置包括:
第二测井数据计算单元,用于根据地震叠加速度及待测区域有井地区的第一测井数据,计算待测区域无井地区的第二测井数据;
压实系数计算单元,用于利用第二测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
一实施例中,基于地震叠加速度的压实系数确定装置还包括:
压实系数计算模型生成单元,用于根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
一实施例中,压实系数计算模型生成单元包括:
系数计算模型,用于根据多个第一测井数据及相应的深度数据,计算H-V速度计算公式的相变因子系数及压实系数;
压实系数计算模型生成模型,用于利用H-V速度计算公式,根据相变因子系数及压实系数,生成压实系数计算模型。
一实施例中,第二测井数据计算单元包括:
测井数据计算模型获得模型,用于将地震叠加速度及相应的第一测井数据进行拟合,获得测井数据计算模型;
第二测井数据计算模型,用于利用测井数据计算模型,根据地震叠加速度计算待测区域无井地区的第二测井数据。
第三方面,本发明提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行程序时实现基于地震叠加速度的压实系数确定方法的步骤。
第四方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现基于地震叠加速度的压实系数确定方法的步骤。
从上述描述可知,本发明提供基于地震叠加速度的压实系数确定方法及装置,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的实施例中的基于地震叠加速度的压实系数确定方法的一流程示意图;
图2为本发明的实施例中的基于地震叠加速度的压实系数确定方法的另一流程示意图;
图3为本发明的实施例中步骤200的流程示意图;
图4为本发明的实施例中步骤100的流程示意图;
图5为本发明的具体应用实例中基于地震叠加速度的压实系数确定方法的流程示意图;
图6为本发明的具体应用实例中同一地层不同深度地震速度与测井速度比值图;
图7为本发明的具体应用实例中速度随深度变化图;
图8为本发明的实施例中的基于地震叠加速度的压实系数确定装置的结构示意图;
图9为本发明的实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的实施例提供一种基于地震叠加速度的压实系数确定方法的具体实施方式,参见图1,该方法具体包括如下内容:
步骤100:根据地震叠加速度及待测区域有井地区的测井数据,计算待测区域无井地区的测井数据。
在地震勘探中,地震叠加速度通常指多次覆盖技术中的地震叠加速度。根据共反射点时距曲线求得的速度叫做叠加速度。求取叠加速度有两种基本方法,相关方法和叠加方法。
可以理解的是,针对同一区块,其地层地质特征固定,地震叠加速度与测井数据必然存在某种联系,具体来说,地震叠加速度与测井数据中的地层速度或者利用测井数据来求取的地层速度存在特定关系,另外,有井地区及无井地区同属于待测区域。
步骤200:利用待测区域无井地区的测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
可以理解的是步骤200中的待测区域无井地区的测井数据,并不是通过测井仪器测量出的实际测井数据,而是通过地震叠加速度反推出的待测区域无井地区的预测测井数据。
从上述描述可知,本发明提供基于地震叠加速度的压实系数确定方法,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
一实施例中,参见图2,基于地震叠加速度的压实系数确定方法还包括:
步骤200a:根据多个待测区域有井地区的测井数据及相应的深度数据,利用H-V(High-Velocity)速度计算公式,生成压实系数计算模型。
可以理解的是,待测区域有井地区的测井数据是通过测井仪器实际测量出的测井数据。压实系数是指地层岩石随着地层埋深深度的增加,岩石体积和孔隙受压减小的表征参数,故可以选用H-V速度计算公式来生成压实系数计算模型。
一实施例中,参见图3,步骤200a包括:
步骤201:根据多个待测区域有井地区的测井数据及相应的深度数据,计算H-V速度计算公式的相变因子系数及压实系数。
可以理解的是,H-V速度计算公式是个通用公式,可以应用于不同地区,但针对不同地区时,其公式中的相变因子系数及压实系数不相同,故需要具体计算。
步骤202:利用H-V速度计算公式,根据相变因子系数及压实系数,生成压实系数计算模型。
将步骤201所计算出的相变因子系数及压实系数带入到H-V速度计算公式(通用公式)中,获得待测区域的压实系数计算模型。
一实施例中,参见图4,步骤100包括:
步骤101:将地震叠加速度及相应的待测区域有井地区的测井数据进行拟合,获得测井数据计算模型。
步骤102:利用测井数据计算模型,根据地震叠加速度计算待测区域无井地区的测井数据。
步骤101至步骤102中,通过待测地区的地震叠加速度及有井地区的测井数据建立起待测地区的地震叠加速度及测井数据的关系式,并通过该关系式及地震叠加速度计算待测地区的无井地区的测井数据,可以理解的是,地震叠加速度具有平面意义,即该数据覆盖待测地区,而测井数据相对来说是“点数据”,只能代表一定面积的地下地质特征。
为进一步地说明本方案,本发明以川东H井区为例,提供基于地震叠加速度的压实系数确定方法的具体应用实例,该具体应用实例具体包括如下内容,参见图5。
S0:获取川东H井区有井地区的测井数据及地震叠加速度。
步骤S0中的测井数据为常规测井曲线数据,包括:三岩性测井曲线数据、三孔隙度测井曲线数据及三电性测井曲线数据,一实施例中,测井曲线数据可以为:伽马曲线、井径曲线及自然电位曲线(三岩性测井曲线数据);中子孔隙度曲线、声波时差测井曲线、密度测井曲线(三孔隙度测井曲线数据);深侧向电阻率、浅侧向电阻率曲线及微球型聚焦电阻率(三电性测井曲线数据)。
S1:建立川东H井区有井地区的测井速度与地震叠加速度之间的关系式。
通过读取该地区4个位置地震叠加速度与测井速度关系(表1及图6),建立地震叠加速度与测井速度关系式:
V测井=(11500*V地震)/(H+12000) (1)
其中:V测井:井点处测井速度;V地震:井点处地震速度;H为地层深度。
表1川东H井区地震叠加速度与测井速度拟合关系
深度(m) | 地震速度(m/s) | 测井速度(m/s) | 比例因子 |
800 | 4000 | 3824 | 1.05 |
1400 | 4300 | 3992 | 1.08 |
2000 | 4750 | 4160 | 1.14 |
2500 | 5100 | 4300 | 1.19 |
可以理解的是,可以利用步骤S0中的声波时差测井曲线(或密度测井曲线)来计算测井速度。
S2:生成川东H井区的压实系数计算模型。
根据多个待测区域有井地区的测井速度及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
基于沉积岩学知识和岩石物理理论,可以近似认为同一地区的地层速度随埋深的增加而增加。
H-V速度计算公式:V=a*H+b (2)
V为地层速度,m/s;
H为深度,m;
a为压实系数,主要与埋深、岩性有关;
b为相变因子,主要与沉积相、流体势场有关(一个地区该值为定值)。
故a=(V-b)/H (3)
表2川东H井区测井速度与深度关系
井名 | 速度(m/s) | 地层深度(m) |
C8 | 3992 | 1400 |
C27 | 4050 | 1450 |
C34 | 4030 | 1625 |
C28 | 4100 | 1680 |
C23 | 4200 | 2100 |
通过以上5口井,拟合出速度随深度变化关系图(图7)
可以求取压实系数a:0.283相变因子b:3625
由于压实a值与地层埋深有关,无井区a可以通过地震速度求取;相变因子b值与沉积相有关,所以该地区b值确定为3625。
故川东H井区压实系数计算模型为:a=(V测井-3625)/H。 (5)
S3:计算川东H井区无井地区的压实系数。
结合公式(1)及公式(5),计算川东H井区无井地区的压实系数,可以理解的是,鉴于地震叠加速度数据覆盖川东H井区,故S3中可以计算川东H井区无井地区的“点”的压实系数,也可以计算川东H井区无井地区“面”的压实系数平面图分布图。
S4:约束川东H井区无井地区压实系数平面图。
利用川东H井区的有井地区的压实系数,并结合地震叠加速度约束压实系数平面图分布图,以最大限度恢复地层真实厚度。
从上述描述可知,本发明提供基于地震叠加速度的压实系数确定方法,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
基于同一发明构思,本申请实施例还提供了基于地震叠加速度的压实系数确定装置,可以用于实现上述实施例所描述的方法,如下面的实施例。由于基于地震叠加速度的压实系数确定装置解决问题的原理与基于地震叠加速度的压实系数确定方法相似,因此基于地震叠加速度的压实系数确定装置的实施可以参见基于地震叠加速度的压实系数确定方法实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
本发明的实施例提供一种能够实现基于地震叠加速度的压实系数确定方法的基于地震叠加速度的压实系数确定装置的具体实施方式,参见图8,基于地震叠加速度的压实系数确定装置具体包括如下内容:
第二测井数据计算单元10,用于根据地震叠加速度及待测区域有井地区的第一测井数据,计算待测区域无井地区的第二测井数据;
压实系数计算单元20,用于利用第二测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
一实施例中,基于地震叠加速度的压实系数确定装置还包括:
压实系数计算模型生成单元,用于根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
一实施例中,压实系数计算模型生成单元包括:
系数计算模型,用于根据多个第一测井数据及相应的深度数据,计算H-V速度计算公式的相变因子系数及压实系数;
压实系数计算模型生成模型,用于利用H-V速度计算公式,根据相变因子系数及压实系数,生成压实系数计算模型。
一实施例中,第二测井数据计算单元包括:
测井数据计算模型获得模型,用于将地震叠加速度及相应的第一测井数据进行拟合,获得测井数据计算模型;
第二测井数据计算模型,用于利用测井数据计算模型,根据地震叠加速度计算待测区域无井地区的第二测井数据。
从上述描述可知,本发明提供基于地震叠加速度的压实系数确定装置,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
本申请的实施例还提供能够实现上述实施例中的基于地震叠加速度的压实系数确定方法中全部步骤的一种电子设备的具体实施方式,参见图9,电子设备具体包括如下内容:
处理器(processor)1201、存储器(memory)1202、通信接口(CommunicationsInterface)1203和总线1204;
其中,处理器1201、存储器1202、通信接口1203通过总线1204完成相互间的通信;通信接口1203用于实现服务器端设备、测量设备以及用户端设备等相关设备之间的信息传输。
处理器1201用于调用存储器1202中的计算机程序,处理器执行计算机程序时实现上述实施例中的基于地震叠加速度的压实系数确定方法中的全部步骤,例如,处理器执行计算机程序时实现下述步骤:
步骤100:根据地震叠加速度及待测区域有井地区的测井数据,计算待测区域无井地区的测井数据。
步骤200:利用待测区域无井地区的测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
从上述描述可知,本申请实施例中的电子设备,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
本申请的实施例还提供能够实现上述实施例中的基于地震叠加速度的压实系数确定方法中全部步骤的一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的基于地震叠加速度的压实系数确定方法的全部步骤,例如,处理器执行计算机程序时实现下述步骤:
步骤100:根据地震叠加速度及待测区域有井地区的测井数据,计算待测区域无井地区的测井数据。
步骤200:利用待测区域无井地区的测井数据及预先建立的压实系数计算模型计算待测区域无井地区的压实系数。
从上述描述可知,本申请实施例中的计算机可读存储介质,通过有井地区的测井速度求取压实系数,并建立有井地区的地震叠加速度和测井速度关系,从而间接建立地震叠加速度和压实系数的关系,这样在无井区域可以间接通过地震叠加速度来求取压实校正系数,并通过无井区地震速度、有井地区压实系数来约束压实系数的计算结果,以最大限度恢复地层真实厚度。本方法适用于勘探初期阶段,井网控制程度比较低,难以通过有井地区准确求取地层压实系数,地层恢复不准确,从而避免了传统方法的必须应用于井网密集地区的缺点。综上,本方法对无井区地层压实系数求取这个核心关键过程进行了更为有益的改进,对压实系数平面规律分析过程进行了优化,使该发明方法能够适用于多期次构造运动变化的地震工区,且方法实现过程简单,易于编程实现。方法应用获得的成果具有较高的合理性和可靠性,为构造演化分析提供扎实的数据基础,对提高勘探初期油气成藏分析具有重要的应用价值。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于硬件+程序类实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书的实施例可提供为方法、系统或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已,并不用于限制本说明书实施例。对于本领域技术人员来说,本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本说明书实施例的权利要求范围之内。
Claims (10)
1.一种基于地震叠加速度的压实系数确定方法,其特征在于,包括:
根据地震叠加速度及待测区域有井地区的第一测井数据,计算所述待测区域无井地区的第二测井数据;
利用所述第二测井数据及预先建立的压实系数计算模型计算所述待测区域无井地区的压实系数;
所述根据地震叠加速度及待测区域有井地区的第一测井数据,计算所述待测区域无井地区的第二测井数据,包括:
建立有井地区的测井速度与地震叠加速度之间的关系式;利用声波时差测井曲线或密度测井曲线来计算测井速度,通过读取地区地震叠加速度与测井速度关系,
建立地震叠加速度与测井速度关系式:
V测井=(11500*V地震)/(H+12000) (1)
其中:V测井:井点处测井速度;V地震:井点处地震速度;H为地层深度;
利用所述第二测井数据及预先建立的压实系数计算模型计算所述待测区域无井地区的压实系数,包括:
根据多个待测区域有井地区的测井速度及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型:
H-V速度计算公式:V=a*H+b (2)
V为地层速度,m/s;
H为深度,m;
a为压实系数,与埋深、岩性有关;
b为相变因子,与沉积相、流体势场有关,一个地区为定值;
a=(V-b)/H (3)
结合公式(1)及公式(3),计算无井地区的压实系数。
2.如权利要求1所述的压实系数确定方法,其特征在于,还包括:根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
3.如权利要求2所述的压实系数确定方法,其特征在于,所述根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型,包括:
根据多个第一测井数据及相应的深度数据,计算所述H-V速度计算公式的相变因子系数及压实系数;
利用所述H-V速度计算公式,根据所述相变因子系数及压实系数,生成压实系数计算模型。
4.如权利要求1所述的压实系数确定方法,其特征在于,所述根据地震叠加速度及待测区域有井地区的第一测井数据,计算所述待测区域无井地区的第二测井数据,包括:
将所述地震叠加速度及相应的所述第一测井数据进行拟合,获得测井数据计算模型;
利用所述测井数据计算模型,根据所述地震叠加速度计算所述待测区域无井地区的第二测井数据。
5.一种基于地震叠加速度的压实系数确定装置,其特征在于,包括:
第二测井数据计算单元,用于根据地震叠加速度及待测区域有井地区的第一测井数据,计算所述待测区域无井地区的第二测井数据;
压实系数计算单元,用于利用所述第二测井数据及预先建立的压实系数计算模型计算所述待测区域无井地区的压实系数;
所述根据地震叠加速度及待测区域有井地区的第一测井数据,计算所述待测区域无井地区的第二测井数据,包括:
建立有井地区的测井速度与地震叠加速度之间的关系式;利用声波时差测井曲线或密度测井曲线来计算测井速度,通过读取地区地震叠加速度与测井速度关系,
建立地震叠加速度与测井速度关系式:
V测井=(11500*V地震)/(H+12000) (1)
其中:V测井:井点处测井速度;V地震:井点处地震速度;H为地层深度;
利用所述第二测井数据及预先建立的压实系数计算模型计算所述待测区域无井地区的压实系数,包括:
根据多个待测区域有井地区的测井速度及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型:
H-V速度计算公式:V=a*H+b(2)
V为地层速度,m/s;
H为深度,m;
a为压实系数,主要与埋深、岩性有关;
b为相变因子,与沉积相、流体势场有关,一个地区为定值;
a=(V-b)/H (3)
结合公式(1)及公式(3),计算无井地区的压实系数。
6.如权利要求5所述的压实系数确定装置,其特征在于,还包括:
压实系数计算模型生成单元,用于根据多个第一测井数据及相应的深度数据,利用H-V速度计算公式,生成压实系数计算模型。
7.如权利要求6所述的压实系数确定装置,其特征在于,包括:所述压实系数计算模型生成单元包括:
系数计算模型,用于根据多个第一测井数据及相应的深度数据,计算所述H-V速度计算公式的相变因子系数及压实系数;
压实系数计算模型生成模型,用于利用所述H-V速度计算公式,根据所述相变因子系数及压实系数,生成压实系数计算模型。
8.如权利要求5所述的压实系数确定装置,其特征在于,所述第二测井数据计算单元包括:
测井数据计算模型获得模型,用于将所述地震叠加速度及相应的所述第一测井数据进行拟合,获得测井数据计算模型;
第二测井数据计算模型,用于利用所述测井数据计算模型,根据所述地震叠加速度计算所述待测区域无井地区的第二测井数据。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4任一项所述基于地震叠加速度的压实系数确定方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4任一项所述基于地震叠加速度的压实系数确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910363616.8A CN111856567B (zh) | 2019-04-30 | 2019-04-30 | 基于地震叠加速度的压实系数确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910363616.8A CN111856567B (zh) | 2019-04-30 | 2019-04-30 | 基于地震叠加速度的压实系数确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111856567A CN111856567A (zh) | 2020-10-30 |
CN111856567B true CN111856567B (zh) | 2023-09-26 |
Family
ID=72965778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910363616.8A Active CN111856567B (zh) | 2019-04-30 | 2019-04-30 | 基于地震叠加速度的压实系数确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111856567B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0397313A2 (en) * | 1989-03-23 | 1990-11-14 | Kenneth Rayvenor Lusty Armitage | Comprehensive system for evaluation of velocities of acoustical energy of sedimentary rocks |
US5475589A (en) * | 1992-07-08 | 1995-12-12 | Spiral Holding, Ltd. | System for evaluating seismic sequence lithology and property, and for evaluating risk associated with predicting potential hydrocarbon reservoir, seal, trap or source |
US6278948B1 (en) * | 1999-04-02 | 2001-08-21 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data |
CN101533103A (zh) * | 2009-04-13 | 2009-09-16 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
CN104678440A (zh) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | 井约束二维地震变速度场非线性误差校正方法 |
WO2015196779A1 (zh) * | 2014-06-24 | 2015-12-30 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN107783187A (zh) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | 一种将测井速度和地震速度结合建立三维速度场的方法 |
CN107843927A (zh) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | 基于井震联合速度的页岩地层压力预测方法及装置 |
CN107942384A (zh) * | 2017-11-14 | 2018-04-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用三维地震资料实现定向井智能导向的预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10120092B2 (en) * | 2016-08-22 | 2018-11-06 | Chevron U.S.A. Inc. | System and method for analysis of depositional settings of subsurface reservoirs |
-
2019
- 2019-04-30 CN CN201910363616.8A patent/CN111856567B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0397313A2 (en) * | 1989-03-23 | 1990-11-14 | Kenneth Rayvenor Lusty Armitage | Comprehensive system for evaluation of velocities of acoustical energy of sedimentary rocks |
US5475589A (en) * | 1992-07-08 | 1995-12-12 | Spiral Holding, Ltd. | System for evaluating seismic sequence lithology and property, and for evaluating risk associated with predicting potential hydrocarbon reservoir, seal, trap or source |
US6278948B1 (en) * | 1999-04-02 | 2001-08-21 | Conoco Inc. | Method for gravity and magnetic data inversion using vector and tensor data |
CN101533103A (zh) * | 2009-04-13 | 2009-09-16 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
WO2015196779A1 (zh) * | 2014-06-24 | 2015-12-30 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN104678440A (zh) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | 井约束二维地震变速度场非线性误差校正方法 |
CN107783187A (zh) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | 一种将测井速度和地震速度结合建立三维速度场的方法 |
CN107843927A (zh) * | 2016-09-20 | 2018-03-27 | 中国石油化工股份有限公司 | 基于井震联合速度的页岩地层压力预测方法及装置 |
CN107942384A (zh) * | 2017-11-14 | 2018-04-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用三维地震资料实现定向井智能导向的预测方法 |
Non-Patent Citations (12)
Title |
---|
4d time strain and the seismic signature of geomechanical compaction at genesis;James Rickett等;《Leading Edge》;第26卷(第5期);第644-647页 * |
Quantitative estimation of compaction and velocity changes using 4d impedance and traveltime changes;Landr M等;《Geophysics》;第69卷(第4期);第949-957页 * |
一个深水区深层沉积相早期预测方法――以南海北部深水区D凹陷渐新统为例;陈宇航等;《天然气地球科学》;第26卷(第11期);第2122-2130页 * |
利用Eaton法计算地层孔隙压力的不确定性分析;臧艳彬等;《西南石油大学学报(自然科学版)》;第34卷(第04期);第55-61页 * |
北三台地区用地震资料进行泥岩压实和油气运移研究;刘福宁;《石油实验地质》;第15卷(第01期);第74-79页 * |
地震速度谱资料预测地层压力的陷阱分析;张建宁;《勘探地球物理进展》;第28卷(第01期);第51-56页 * |
异常地层压力的预测及其在地质风险评价中的应用;刘福宁;《石油勘探与开发》;第19卷(第06期);第9-14页 * |
沾车地区速度场的建立及应用;王楠;《油气地球物理》;第8卷(第04期);第15-18页 * |
波形指示叠前地震反演方法在致密含油薄砂层预测中的应用;章雄等;《物探与化探》;第42卷(第3期);第545-554页 * |
济阳拗陷地层压力预测方法;谭明友;《石油地球物理勘探》;第39卷(第03期);第314-318页 * |
精细压力预测技术研究;蒋涔等;《海洋石油》;第36卷(第03期);第26-29页 * |
缅甸中央盆地北部新生代古地貌重建;李任远等;《地质科技情报 》;第37卷(第3期);第126-133页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111856567A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang | Incorporating geological conceptual models and interpretations into reservoir modeling using multiple-point geostatistics | |
US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
EA020635B1 (ru) | Способ сейсмической разведки для поиска углеводородов с применением поля средних скоростей сейсмической волны | |
CN106707344B (zh) | 一种地层层序划分方法及装置 | |
CN103245972A (zh) | 一种在二维空间内确定复杂地质构造的方法 | |
CN106646603A (zh) | 一种实际时移地震资料处理差异的可靠性判断方法 | |
CN108873065A (zh) | 砂岩优质储层预测方法及装置 | |
CN111751870B (zh) | 叠后层间多次波压制方法及装置 | |
Sancevero et al. | Comparing deterministic and stochastic seismic inversion for thin-bed reservoir characterization in a turbidite synthetic reference model of Campos Basin, Brazil | |
Pochat et al. | A simple method of determining sand/shale ratios from seismic analysis of growth faults: an example from upper Oligocene to lower Miocene Niger Delta deposits | |
Zhou et al. | Bayesian time-lapse difference inversion based on the exact Zoeppritz equations with blockiness constraint | |
CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
CN111856567B (zh) | 基于地震叠加速度的压实系数确定方法及装置 | |
CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
Quinay et al. | Waveform Inversion for Modeling Three‐Dimensional Crust Structure with Topographic Effects | |
CN109991663B (zh) | 工区地震速度体校正方法和装置 | |
AlKawai et al. | Integrating statistical rock physics and pressure and thermal history modeling to map reservoir lithofacies in the deepwater Gulf of Mexico | |
CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
Kumar et al. | Appraisal of reservoir porosity using a machine learning approach: A study from the Eocene‐Miocene interval of Upper Assam Basin, NE, India | |
CN111562629A (zh) | 一种基于等效孔隙截面指数的饱和度确定方法及装置 | |
Song et al. | The impact of local velocity anomaly on structure confirming and the varying velocity mapping strategy | |
CN112835096B (zh) | 气层识别方法及装置 | |
Pereira et al. | Geostatistical seismic inversion for frontier exploration | |
CN112824940B (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 |