CN110967756B - 基于正态分布的微地震定位精度评价方法及系统 - Google Patents
基于正态分布的微地震定位精度评价方法及系统 Download PDFInfo
- Publication number
- CN110967756B CN110967756B CN201811158023.XA CN201811158023A CN110967756B CN 110967756 B CN110967756 B CN 110967756B CN 201811158023 A CN201811158023 A CN 201811158023A CN 110967756 B CN110967756 B CN 110967756B
- Authority
- CN
- China
- Prior art keywords
- microseism
- offset
- normal distribution
- average energy
- positioning
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000006870 function Effects 0.000 claims description 48
- 238000005070 sampling Methods 0.000 claims description 10
- 230000008901 benefit Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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)计算所述等偏移距线上的震源的叠加能量:
其中,E(j)为震源的叠加能量,gk(j)为微地震的监测道记录,n为微地震监测道数,k为道序号,j为微地震监测记录的采样点序号,j=1,2,……,L,L为每道的总采样点数。
优选地,根据所述平均能量与所述偏移距的对应关系,拟合符合正态分布的所述平均能量特征函数。
优选地,所述平均能量特征函数的置信参数为:
其中,σ为置信参数,Ei为偏移距为i时的平均能量,d为偏移距最大值,i为偏移距。
优选地,根据所述置信参数,分析微地震定位精度包括:σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
根据本发明的另一方面,提出了一种基于正态分布的微地震定位精度评价系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据微地震的监测记录,确定微地震定位的定位位置;根据所述定位位置,求取不同偏移距处的平均能量,进而获得平均能量特征函数,所述平均能量特征函数服从正态分布;计算所述平均能量特征函数的置信参数,根据所述置信参数,分析微地震定位精度。
优选地,根据所述定位位置,求取不同偏移距处的平均能量包括:以所述定位位置为圆心,以所述偏移距为半径画圆,获得等偏移距线;计算所述等偏移距线上多个震源的叠加能量;计算多个震源的叠加能量的平均值,即为所述偏移距处的平均能量。
优选地,通过公式(1)计算所述等偏移距线上的震源的叠加能量:
其中,E(j)为震源的叠加能量,gk(j)为微地震的监测道记录,n为微地震监测道数,k为道序号,j为微地震监测记录的采样点序号,j=1,2,……,L,L为每道的总采样点数。
优选地,根据所述平均能量与所述偏移距的对应关系,拟合符合正态分布的所述平均能量特征函数。
优选地,所述平均能量特征函数的置信参数为:
其中,σ为置信参数,Ei为偏移距为i时的平均能量,d为偏移距最大值,i为偏移距。
优选地,根据所述置信参数,分析微地震定位精度包括:σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
其有益效果在于:
(1)运算的速度快,效率高;
(2)相对简单易于实现;
(3)通过微地震定位结果能量与偏移距构建正态分布函数,计算尺度函数来求取质控评价参数,对每个事件进行定位精度评价;通过阈值来筛选定位精度低的事件,从而实现对定位结果精度评价。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的基于正态分布的微地震定位精度评价方法的步骤的流程图。
图2示出了根据本发明的一个实施例的微地震的监测记录的示意图。
图3示出了根据本发明的一个实施例的平均能量特征函数的示意图。
图4示出了根据本发明的一个实施例的平均能量特征函数服从正态分布的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的基于正态分布的微地震定位精度评价方法的步骤的流程图。
在该实施例中,根据本发明的基于正态分布的微地震定位精度评价方法可以包括:步骤101,根据微地震的监测记录,确定微地震定位的定位位置;步骤102,根据定位位置,求取不同偏移距处的平均能量,进而获得平均能量特征函数,平均能量特征函数服从正态分布;步骤103,计算平均能量特征函数的置信参数,根据置信参数,分析微地震定位精度。
在一个示例中,根据定位位置,求取不同偏移距处的平均能量包括:以定位位置为圆心,以偏移距为半径画圆,获得等偏移距线;计算等偏移距线上多个震源的叠加能量;计算多个震源的叠加能量的平均值,即为偏移距处的平均能量。
在一个示例中,通过公式(1)计算等偏移距线上的震源的叠加能量:
其中,E(j)为震源的叠加能量,gk(j)为微地震的监测道记录,n为微地震监测道数,k为道序号,j为微地震监测记录的采样点序号,j=1,2,……,L,L为每道的总采样点数。
在一个示例中,根据平均能量与偏移距的对应关系,拟合符合正态分布的平均能量特征函数。
在一个示例中,平均能量特征函数的置信参数为:
其中,σ为置信参数,Ei为偏移距为i时的平均能量,d为偏移距最大值,i为偏移距。
在一个示例中,根据置信参数,分析微地震定位精度包括:σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
具体地,在微地震定位是微地震处理中最重要的环节,定位的精度也是验证微地震处理结果直接参数。微地震定位方法较多,不同监测方式其定位方法也不同,根据监测方式和资料的特点,选取适宜的定位方法求取震源的真实解。
从微地震事件能量叠加原理可知,一般越接近震源的真实位置,其叠加的能量值最大,距离真实位置越远,其叠加能量值越小,如果微地震事件定位精度高,那么叠加能量越集中,周围能量越小;定位精度越低,中心叠加能量与周围叠加能量越接近;因此可以根据震源扫描能量叠加原理进行微地震事件定位精度的评价。
根据本发明的基于正态分布的微地震定位精度评价方法可以包括:
根据微地震的监测记录,确定微地震定位的定位位置;以定位位置为圆心,以偏移距为半径画圆,获得等偏移距线;计算等偏移距线上多个震源的叠加能量,可以在等偏移距线上以定位位置为中心做出等角度线,例如十字线、60度线等,将等偏移距线平均分为4份、6份等,计算等偏移距线与等角度线的交点上的震源的叠加能量;计算多个震源的叠加能量的平均值,即为偏移距处的平均能量,根据平均能量与偏移距的对应关系,拟合符合正态分布的平均能量特征函数。
正态分布具有两个参数μ和σ2的连续型随机变量的分布,第一参数μ是服从正态分布的随机变量的均值,是正态分布的位置参数,描述正态分布的集中趋势位置。概率规律为取与μ邻近的值的概率大,而取离μ越远的值的概率越小。第二个参数σ2是此随机变量的方差,其中置信参数σ描述正态分布资料数据分布的离散程度,σ越大,数据分布越分散,σ越小,数据分布越集中,也称为是正态分布的形状参数,σ越大,曲线越扁平,反之,σ越小,曲线越瘦高。
若各个偏移距的平均能量特征函数E近似服从一个数学期望为μ,尺度参数为σ的概率分布,则其概率密度函数为:
则各个偏移距的平均叠加能量E函数为正态随机变量,服从正态分布,这里的d为偏移距,那么可以将尺度参数σ作为微地震事件定位的精度参数。
根据微地震定位的定位位置和随偏移距d能量叠加值序列E1、E2、……、Ed,通过公式(2)计算平均能量特征函数的置信参数,根据置信参数,分析微地震定位精度:σ越大,平均能量特征函数分布越分散,能量叠加值不突出,曲线越扁,微地震定位精度越低,σ越小,平均能量特征函数越集中,能量叠加值突出,曲线越瘦高,微地震定位精度越高。还可以通过设定阈值来筛选定位精度低的事件,从而实现对定位结果精度评价。
本方法通过平均能量特征函数服从正态分布的规律,实现微地震定位精度评价,运算的速度快,效率高,相对简单易于实现。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
以一口地面水力压裂微地震监测资料为实例,对其中一个微地震事件进行震源定位,其中,图2所示的微地震事件的监测记录是通过射孔时差动校正之后的记录,以便于进行微地震事件的识别。
定位方法选取震源扫描能量叠加方法,求取最大叠加能量其所对应的位置为震源的位置。
根据定位位置为圆心,以偏移距50米为半径画圆,获得60度线等偏移距线;计算等偏移距线上多个震源的叠加能量;形成二维平面能量叠加谱,如图3所示,由于选取的数据是地面监测的方式,水平方向定位精度高,选取x,y两个方向为变量,深度为一固定值,计算不同偏移距的选取的各个震源的能量值。
计算同一偏移距上多个震源的叠加能量的平均值Ed,即为偏移距d处的平均能量。根据微地震位置和随偏移距d能量叠加值序列E1、E2、……、Ed平均能量与偏移距的对应关系,拟合符合正态分布的平均能量特征函数,如图4所示。通过公式(2)计算平均能量特征函数的置信函数σ,根据σ的大小,分析微地震定位精度:σ越大,平均能量特征函数分布越分散,能量叠加值不突出,曲线越扁,微地震定位精度越低,σ越小,平均能量特征函数越集中,能量叠加值突出,曲线越瘦高,微地震定位精度越高。
综上所述,本发明通过平均能量特征函数服从正态分布的规律,实现微地震定位精度评价,运算的速度快,效率高,相对简单易于实现。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种基于正态分布的微地震定位精度评价系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据微地震的监测记录,确定微地震定位的定位位置;根据定位位置,求取不同偏移距处的平均能量,进而获得平均能量特征函数,平均能量特征函数服从正态分布;计算平均能量特征函数的置信参数,根据置信参数,分析微地震定位精度。
在一个示例中,根据定位位置,求取不同偏移距处的平均能量包括:以定位位置为圆心,以偏移距为半径画圆,获得等偏移距线;计算等偏移距线上多个震源的叠加能量;计算多个震源的叠加能量的平均值,即为偏移距处的平均能量。
在一个示例中,通过公式(1)计算等偏移距线上的震源的叠加能量:
其中,E(j)为震源的叠加能量,gk(j)为微地震的监测道记录,n为微地震监测道数,k为道序号,j为微地震监测记录的采样点序号,j=1,2,……,L,L为每道的总采样点数。
在一个示例中,根据平均能量与偏移距的对应关系,拟合符合正态分布的平均能量特征函数。
在一个示例中,平均能量特征函数的置信参数为:
其中,σ为置信参数,Ei为偏移距为i时的平均能量,d为偏移距最大值,i为偏移距。
在一个示例中,根据置信参数,分析微地震定位精度包括:σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
本系统通过平均能量特征函数服从正态分布的规律,实现微地震定位精度评价,运算的速度快,效率高,相对简单易于实现。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种基于正态分布的微地震定位精度评价方法,其特征在于,包括:
根据微地震的监测记录,确定微地震定位的定位位置;
根据所述定位位置,求取不同偏移距处的平均能量,进而获得平均能量特征函数,所述平均能量特征函数服从正态分布;
计算所述平均能量特征函数的置信参数,根据所述置信参数,分析微地震定位精度。
2.根据权利要求1所述的基于正态分布的微地震定位精度评价方法,其中,根据所述定位位置,求取不同偏移距处的平均能量包括:
以所述定位位置为圆心,以所述偏移距为半径画圆,获得等偏移距线;
计算所述等偏移距线上多个震源的叠加能量;
计算多个震源的叠加能量的平均值,即为所述偏移距处的平均能量。
4.根据权利要求1所述的基于正态分布的微地震定位精度评价方法,其中,根据所述平均能量与所述偏移距的对应关系,拟合符合正态分布的所述平均能量特征函数。
6.根据权利要求1所述的基于正态分布的微地震定位精度评价方法,其中,根据所述置信参数,分析微地震定位精度包括:
σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
7.一种基于正态分布的微地震定位精度评价系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
根据微地震的监测记录,确定微地震定位的定位位置;
根据所述定位位置,求取不同偏移距处的平均能量,进而获得平均能量特征函数,所述平均能量特征函数服从正态分布;
计算所述平均能量特征函数的置信参数,根据所述置信参数,分析微地震定位精度。
8.根据权利要求7所述的基于正态分布的微地震定位精度评价系统,其中,根据所述定位位置,求取不同偏移距处的平均能量包括:
以所述定位位置为圆心,以所述偏移距为半径画圆,获得等偏移距线;
计算所述等偏移距线上多个震源的叠加能量;
计算多个震源的叠加能量的平均值,即为所述偏移距处的平均能量。
10.根据权利要求7所述的基于正态分布的微地震定位精度评价系统,其中,根据所述置信参数,分析微地震定位精度包括:
σ越大表示微地震定位精度越低,σ越小表示微地震定位精度越高。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811158023.XA CN110967756B (zh) | 2018-09-30 | 2018-09-30 | 基于正态分布的微地震定位精度评价方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811158023.XA CN110967756B (zh) | 2018-09-30 | 2018-09-30 | 基于正态分布的微地震定位精度评价方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110967756A CN110967756A (zh) | 2020-04-07 |
CN110967756B true CN110967756B (zh) | 2021-09-17 |
Family
ID=70028863
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811158023.XA Active CN110967756B (zh) | 2018-09-30 | 2018-09-30 | 基于正态分布的微地震定位精度评价方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110967756B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4330872A (en) * | 1974-09-23 | 1982-05-18 | Mobil Oil Corporation | Common-offset-distance seismic trace filtering |
GB2357097A (en) * | 1999-12-08 | 2001-06-13 | Norske Stats Oljeselskap | Method of assessing positional uncertainty in drilling a well |
CN102879801B (zh) * | 2012-08-30 | 2015-07-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于射孔约束的EnKF微地震事件位置反演方法 |
CN103913774B (zh) * | 2014-04-02 | 2017-03-01 | 西南石油大学 | 基于微地震事件的储层地质力学参数反演方法 |
CN105093297B (zh) * | 2015-07-10 | 2017-06-16 | 中联煤层气有限责任公司 | 一种微地震定位精度的质量控制方法 |
CN106772577B (zh) * | 2016-06-29 | 2019-04-26 | 中国石油大学(华东) | 基于微地震数据和spsa优化算法的震源反演方法 |
-
2018
- 2018-09-30 CN CN201811158023.XA patent/CN110967756B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110967756A (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210041596A1 (en) | Petrophysical Inversion With Machine Learning-Based Geologic Priors | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN104502997B (zh) | 一种利用裂缝密度曲线预测裂缝密度体的方法 | |
US10809416B2 (en) | Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells | |
CN104280775B (zh) | 一种基于全波形矢量偏移叠加的微地震监测定位方法 | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN108333627B (zh) | 火成岩地区断裂真伪的识别方法及装置 | |
WO2016008105A1 (zh) | 一种基于柯西分布的叠后波阻抗反演方法 | |
CN104200115A (zh) | 一种基于地质统计学模拟的全地层速度建模方法 | |
CN107132578A (zh) | 一种微地震地面监测速度模型校正算法 | |
CN103869362A (zh) | 体曲率获取方法和设备 | |
CN102879815B (zh) | 基于空间型灰度伴随矩阵的结构性地震属性提取方法 | |
CN108957554A (zh) | 一种地球物理勘探中地震反演方法 | |
CN110967756B (zh) | 基于正态分布的微地震定位精度评价方法及系统 | |
CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
CN109425889B (zh) | 一种用于刻画古岩溶暗河的方法 | |
CN110068861B (zh) | 基于对称观测的vsp缝洞识别技术方法 | |
CN107942373A (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
CN103064113A (zh) | 确定地下速度结构的双参数展开方法及装置 | |
CN112987091A (zh) | 储层检测方法、装置、电子设备和存储介质 | |
CN111474580A (zh) | 一种基于炮检距矢量片的方位角道集提取方法和系统 | |
CN113296147B (zh) | 一种识别地震有限断层破裂参数的方法及系统 | |
CN109581521A (zh) | Tti各向异性的局部层析方法及系统 | |
CN109655906A (zh) | 地震相干属性计算方法及系统 | |
CN118050806A (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 |