CN109521470B - 分析地质构造对地震反演裂缝密度影响的方法 - Google Patents
分析地质构造对地震反演裂缝密度影响的方法 Download PDFInfo
- Publication number
- CN109521470B CN109521470B CN201811635126.0A CN201811635126A CN109521470B CN 109521470 B CN109521470 B CN 109521470B CN 201811635126 A CN201811635126 A CN 201811635126A CN 109521470 B CN109521470 B CN 109521470B
- Authority
- CN
- China
- Prior art keywords
- seismic
- fracture
- density
- seismic data
- gather
- 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 39
- 239000011159 matrix material Substances 0.000 claims abstract description 40
- 239000011435 rock Substances 0.000 claims abstract description 8
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 4
- 230000000694 effects Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000005755 formation reaction Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 2
- 230000005012 migration Effects 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 5
- 238000012937 correction Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral 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/30—Analysis
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
-
- 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/66—Subsurface modeling
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
本发明提供了分析地质构造对地震反演裂缝密度影响的方法，包括以下步骤，步骤S1，基于待研究目的区域地球物理资料和地质资料，构建含地质构造的地震数值模型；步骤S2，基于的地震数值模型，定义地震观测系统，记录生成合成地震数据；步骤S3，通过对步骤S2中的地震数据处理，以得到每个偏移距集合的平均振幅矩阵，近似形成目的层反射系数矩阵；步骤S4，构建岩石物理模型，以得到基函数系列和裂缝属性系列；步骤S5，基于步骤S3获取的目的层反射系数矩阵和步骤S4获取的基函数系列，合成地震数据的裂缝属性矩阵序列；步骤S6，对比步骤S4获取的裂缝属性系列和步骤S5的地震数据的裂缝属性矩阵序列，以获取合成地震数据的裂缝密度。
Description
技术领域
本发明涉及地震反演裂缝分析技术领域，更具体地，涉及一种分析地质构造对地震反演裂缝密度影响的方法。
背景技术
纵波在各向异性介质中传播，纵波的属性有各种变化规律，利用这些规律，可反演介质各向异性的参数，裂缝介质是各向异性介质的一种特殊类型，利用纵波属性的变化规律可反演裂缝的密度和走向。这些纵波属性主要包括反射振幅、动校正速度、旅行时和衰减等。
其中，纵波振幅属性，平面波在两种各向同性介质分界面处反射系数和透射系数的经典解析表达式为Zoeppritz方程(以人名Zoeppritz命名的反射投射系数方程)，它是各向同性介质中AVO(Amplitude Versus Offset，振幅随偏移距变化)分析的理论基础。由于解析表达式中不方便的参数化，后来学者们已经开发了一系列涵盖了合理的入射范围的近似。在实际应用中，AVO响应中的特征可以帮助确定岩石特性；另外，裂缝介质中的AVO分析通常考虑方位角变化。各向同性与HTI(horizontal transverse isotropy，水平横向各向异性)分界面处的纵波反射系数随入射角和方位角变化，这种变化规律揭示了纵波反射系数和AVO梯度都表现出随方位角为椭圆变化规律。椭圆率的变化表示地震各向异性的强度。换句话说，假设应用该各向同性/HTI模型，则从地震数据提取的反射系数可用于反演地震各向异性(或相对裂缝密度)的主裂缝方向和强度。
而纵波旅行时和动校正速度是另外两个重要属性，它们通常相互关联。动校正速度(或慢度)与方位角之间存在椭圆关系，将该技术应用于纵波地震资料的HTI介质，可以反演对称平面方向(裂缝走向)，并通过拟合椭圆变化而计算相对裂缝密度。因此，在纵波资料应用中，通过分析地震道集中的动校正速度，结合椭圆拟合技术可反演裂缝相对密度和裂缝走向。
地震衰减是由透射介质或系统引起的地震波振幅或能量的减小。含油气层储层中的对齐裂纹或裂缝被认为是流体的优质储集场所或通道，导致了所含流体或流体流动的优先方向。该特征进一步导致依赖方位角的衰减，并且衰减与方位角之间有密切椭圆关系，这种椭圆可以用来指导地震资料的纵波衰减裂缝反演。对于裂缝介质中纵波衰减的方位角变化，对于某一偏移距，通常采用谱比法计算不同方位角处的品质因子，然后将椭圆拟合技术应用于品质因子以获得主裂缝走向(对应于拟合椭圆的长轴)和相对裂缝密度(对应于椭圆率)。
目前采用奇异值分解法分析纵波振幅并计算裂缝密度是常见的缝隙方向，其具体操作方式是对于各向异性模型，上部为各向同性介质和下部为裂缝(等效的HTI)介质的交界面的情形，Varela等人(2007)提出的奇异值分解法(SVD)方法，纵波反射系数矩阵的奇异值分解产生如下的基函数系列f和地震属性系列C:
R(fd,θ,φ)＝C1(fd,θ)f1(φ)+C2(fd,θ)f2(φ)+C3(fd,θ)f3(φ) (1)
其中，基函数序列f(包括f1，f2，f3)随方位φ而改变，地震裂缝属性矩阵系列C(包括矩阵C1，C2，C3)与裂缝密度fd和入射角有关联。基函数f仅依赖于方位角，而属性矩阵C依赖于入射角θ和裂缝密度。
在实际裂缝密度反演过程中，首先通过岩石物理模型建立反射系数矩阵R，之后基于SVD方法分解反射系数矩阵得到基函数序列f和属性矩阵序列C，再从实际地震数据中提取的一个新的反射系数矩阵R’，它将与矩阵R共享相同的基函数f，使用R’和基函数f计算新的裂缝属性矩阵序列C’，对比属性矩阵C和新属性矩阵C’即可得出实际地震数据中不同入射角的裂缝密度。
使用地震数据反演裂缝密度的前提是目的层相对平缓，但是实际情况的地质模型目的层很多时候不完全满足目的层相对平缓的条件，比喻说目的层含背斜、向斜和断层等起伏构造，这些构造对裂缝密度反演的影响还有待研究，因此，急需发明一种方法来评估地质起伏构造对裂缝密度反演的影响。
发明内容
本发明提供一种分析地质构造对地震反演裂缝密度影响的新方法。
根据本发明的一个方面，提供一种分析地质构造对地震反演裂缝密度影响的方法，包括以下步骤，
步骤S1，基于待研究目的区域地球物理资料和地质资料，构建含地质构造的地震数值模型；
步骤S2，基于所述的地震数值模型，定义地震观测系统，通过数值有限差分方法模拟合成地震数据，记录生成合成地震数据；
步骤S3，通过对步骤S2中的所述地震数据处理，获取共中心点方位超道集的子道集，基于子道集按不同的偏移距分成若干偏移距集合，以得到每个所述偏移距集合的平均振幅矩阵，近似形成目的层反射系数矩阵；
步骤S4，构建岩石物理模型，获取不同偏移距/不同方位角对应的纵波反射系数矩阵，使用SVD数值分解方法，得到基函数系列和裂缝属性系列；
步骤S5，基于步骤S3获取的所述目的层反射系数矩阵和步骤S4获取的所述基函数系列，合成地震数据的裂缝属性矩阵序列；
步骤S6，对比步骤S4获取的裂缝属性系列和步骤S5的地震数据的裂缝属性矩阵序列，以获取合成地震数据的裂缝密度。
在上述方案基础上优选，所述步骤S1中，所述的地震数值模型包括模型纵横波速度、模型纵横波密度、目的裂缝层异性参数。
在上述方案基础上优选，所述步骤S3中，通过对步骤S2中的所述地震数据进行抽道集、速度分析和叠加与偏移，以获取共中心点方位超道集，并按方位将所述共中心点方位超道集分解成子道集。
在上述方案基础上优选，所述步骤S4中，通过采用SVD方法的公式(1)分解所述不同偏移距/不同方位角对应的纵波反射系数矩，以得到所述基函数系列和裂缝属性系列。
与现有技术相比较，本发明具备以下优点：
本发明的一种分析地质构造对地震反演裂缝密度影响的方法。传统的方法使用Ruger的椭圆拟合算法直接反演，精度低，而且反演的结果是裂缝密度的相对值，不是绝对值相比较；本发明使用SVD方法，反演得到的裂缝密度是绝对值，是真实的裂缝密度信息，更精确的反应大地构造对裂缝反演的影响。
附图说明
图1a表示本发明的模型1的各向异性速度模型；
图1b表示本发明的模型2的各向异性速度模型；
图1c表示本发明的模型3的各向同性速度模型；
图2a表示本发明的模型1某中心点道集对应的反射系数矩阵；
图2b表示本发明的模型2某中心点道集对应的反射系数矩阵；
图2c表示本发明的模型3某中心点道集对应的反射系数矩阵；
图3a表示使用SVD方法构建的基函数系列f1，f2，f3；
图3b表示使用SVD方法构建的地震裂缝属性C1；
图3c表示使用SVD方法构建的地震裂缝属性C2；
图3d表示使用SVD方法构建的地震裂缝属性C3；
图4a表示本发明的模型1的各共中心点位置反演得到的裂缝密度；
图4b表示本发明的模型2的各共中心点位置反演得到的裂缝密度图4c表示本发明的模型3的各共中心点位置反演得到的裂缝密度图5本发明的基于合成地震数据分析地质构造对地震反演裂缝密度的影响的方法流程框图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
请参阅图1所示，本发明提供了一个实际例子说明本发明方法的实施具体方式。本实施例中依据一个实际油田的参数来设计的目的裂缝层非常薄，只有110米，这可能会影响基于SVD方法的反演结果。这里我们设计对应5步实施。
(1)裂缝密度反演时，为了定量评估背斜的影响以及可能的薄层效应(或调谐效果)，我们设计3个模型：模型1(图1a)是直接取自速度剖面(图1c)和包含一个背斜和一层薄裂缝层；模型2(图1b)，由模型1中的第二层和第三层合并得到，参数来自模型1的裂缝层；模型3(图1c)由模型1中的第二层和第三层合并得到，但与模型1的第三层的参数一致。图1a至图1c中不同的颜色表示不同的层，是示意图，其中，白色、绿色和蓝色表示各向同性层，红色表示裂缝层；这三个模型都是三维模型，模型属性在y方向不变，显示的是x-z剖面；x和z都是长度网格个数，单位是10米。这三个模型是三维模型，模型属性在方向不变，显示的是-剖面。在的模型1中，穿过背斜，如果存在反演裂缝密度的变化，可以表明背斜的影响。通过比较模型1和模型2的反演的裂缝密度，可以发现可能的薄层效应。对于模型3，反演的裂缝密度，或假裂缝密度(因为没有裂缝)，如果存在，也显示了背斜的作用。预设裂缝层(红色对应层)密度为0.1，这里使用Chapman岩石物理模型计算对应的裂缝模拟的Thomason参数。在确定以上模型后，定义地震观测系统，使用有限差分方法模拟合成地震数据，记录生成标准的地震数据格式，对于不同的模型，使用相同的观测系统，和模拟方式，便于后续的数据对比。
(2)在生成合成数据后，我们实施地震数据常规处理，完成从炮集到共中心点方位超道集的转换，并变换成方位道集和偏移距集合，每个集合计算出平均振幅，得到了和偏移距/方位角相关的平均振幅矩阵，近似成目的层反射系数矩阵R’，如图2a、图2b和图2c是三个模型某位置对应的共中心点道集目的层反射系数R’例子中，颜色表示反射系数的大小，白色空格表示无效反射系数(或未能抽取到的反射系数)，三个图中，横轴X刻度1-9表示0-180度方位角，每20度计算一个方位角道集，纵轴Y刻度表示5-45度每个1度对应一个入射角格子，白色格子部分表示该位置没有数据道，没有可供计算的振幅。
(3)依据构建岩石物理模型，从测井数据取得各向异性裂缝层与各向同性盖层的弹性参数/Thomasen参数，计算不同入射角和方位角对应的反射系数矩阵R，使用SVD方法按照前面公式(1)，将R分解成f基函数序列和裂缝属性C矩阵序列，如图3。图3是分解后的基函数序列f(f1,f2,f3)和裂缝属性矩阵序列C(C1,C2,C3)。
(4)使用合成地震数据反射系数矩阵R’和上一步通过岩石物理模型得到的基函数序列f，计算合成地震数据的裂缝属性矩阵序列C’，对比合成地震数据的裂缝属性矩阵序列C’和第三步得到的裂缝属性矩阵序列C，其中，裂缝属性一般提取C1，C2和C3，每个参数都可以反演一个裂缝密度fd，综合考虑得到的裂缝密度fd的相对稳定可靠，确定成最终合成地震数据的裂缝密度fd，图4a、图4b和图4c表示三个模型的各中心点位置反演得到的裂缝密度，每个图中左边子图表示不同共中心点位置(x轴)，不同入射角(y轴)对应的反演得到的裂缝密度(颜色值)，右边子图是左子图对应的平均值(x轴表示不同共中心点位置，y轴是裂缝密度，是左边子图不同入射角对应不同裂缝密度的平均值)，也就是实际反演所得到的不同中心点位置最终裂缝密度值，以4a为例，一二三行表示C1，C2和C3三个属性分别对应的裂缝密度反演结果；左列表示，依据各CDP位置和各个入射角位置的C反演的裂缝密度矩阵，右列表示左边图的裂缝密度在不同入射角对应的平均值。
(5)讨论模型图1中背斜对裂缝反演得到裂缝密度的影响。对比图4a和4b可以看出，裂缝层的厚度反演得到的C1，C2和C3没什么大的区别，基本一样，这表明，模型1a中薄层厚度对裂缝反演密度没什么影响；对比图4b和4c我们可以看出，裂缝反演的结果基本验证了预设的裂缝值，特别是图4b的C2反演的结果为0.1左右，和预设的裂缝密度0.1接近，图4c的C2反演的结果为0.0左右，和预设的裂缝密度0.0接近，这说明SVD方法本身可靠；最后一点，也是本发明的创新之处，我们从图4a和4b可以看出，裂缝密度反演在构造背斜的中间的裂缝密度接近真实值0.1，在模型背斜的两侧，反演的裂缝密度大于0.1，这说明本例的设计模型中，背斜的存在导致反演裂缝密度的过高估计(高于模型设置裂缝密度的0.1)。
最后，本申请的方法仅为较佳的实施方案，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种分析地质构造对地震反演裂缝密度影响的方法，其特征在于，包括以下步骤，
步骤S1，基于待研究目的区域地球物理资料和地质资料，构建含地质构造的地震数值模型；
步骤S2，基于所述的地震数值模型，定义地震观测系统，通过数值有限差分方法模拟合成地震数据，记录生成合成地震数据；
步骤S3，通过对步骤S2中的所述合成地震数据处理，获取共中心点方位超道集的子道集，基于子道集按不同的偏移距分成若干偏移距集合，以得到每个所述偏移距集合的平均振幅矩阵，近似形成目的层反射系数矩阵；
步骤S4，构建岩石物理模型，获取不同偏移距/不同方位角对应的纵波反射系数矩阵，以得到基函数系列和裂缝属性系列；
步骤S5，基于步骤S3获取的所述目的层反射系数矩阵和步骤S4获取的所述基函数系列，合成地震数据的裂缝属性矩阵序列；
步骤S6，对比步骤S4获取的裂缝属性系列和步骤S5的合成地震数据的裂缝属性矩阵序列，以获取合成地震数据的裂缝密度。
2.如权利要求1所述的一种分析地质构造对地震反演裂缝密度影响的方法，其特征在于，所述步骤S1中，所述的地震数值模型包括模型纵横波速度、模型纵横波密度、目的裂缝层异性参数。
3.如权利要求1所述的一种分析地质构造对地震反演裂缝密度影响的方法，其特征在于，所述步骤S3中，通过对步骤S2中的所述合成地震数据进行抽道集、速度分析和叠加与偏移，以获取共中心点方位超道集，并按方位将所述共中心点方位超道集分解成子道集。
4.如权利要求1所述的一种分析地质构造对地震反演裂缝密度影响的方法，其特征在于，所述步骤S4中，通过采用SVD方法分解所述不同偏移距/不同方位角对应的纵波反射系数矩，以得到所述基函数系列和裂缝属性系列。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811635126.0A CN109521470B (zh) | 2018-12-29 | 2018-12-29 | 分析地质构造对地震反演裂缝密度影响的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811635126.0A CN109521470B (zh) | 2018-12-29 | 2018-12-29 | 分析地质构造对地震反演裂缝密度影响的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109521470A CN109521470A (zh) | 2019-03-26 |
| CN109521470B true CN109521470B (zh) | 2020-04-14 |
Family
ID=65798181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811635126.0A Active CN109521470B (zh) | 2018-12-29 | 2018-12-29 | 分析地质构造对地震反演裂缝密度影响的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109521470B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113985480B (zh) * | 2021-11-08 | 2024-01-26 | 西南石油大学 | 一种基于角度校正的avo反演方法及装置 |
| CN115576007B (zh) * | 2022-11-22 | 2023-03-14 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102169190A (zh) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
| CN105005079A (zh) * | 2015-07-14 | 2015-10-28 | 北京博达瑞恒科技有限公司 | 一种井曲线反演方法 |
| CN105629303A (zh) * | 2015-12-28 | 2016-06-01 | 中国石油大学(北京) | 基于岩石物理的叠前裂缝定量预测方法及系统 |
| CN106324665A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 反演裂缝密度的方法及系统 |
| CN109307884A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 一种显著反射界面自动提取方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7679993B2 (en) * | 2005-06-17 | 2010-03-16 | Schlumberger Technology Corporation | Method of characterizing a fractured reservoir using seismic reflection amplitudes |
-
2018
- 2018-12-29 CN CN201811635126.0A patent/CN109521470B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102169190A (zh) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
| CN106324665A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 反演裂缝密度的方法及系统 |
| CN105005079A (zh) * | 2015-07-14 | 2015-10-28 | 北京博达瑞恒科技有限公司 | 一种井曲线反演方法 |
| CN105629303A (zh) * | 2015-12-28 | 2016-06-01 | 中国石油大学(北京) | 基于岩石物理的叠前裂缝定量预测方法及系统 |
| CN109307884A (zh) * | 2017-07-28 | 2019-02-05 | 中国石油化工股份有限公司 | 一种显著反射界面自动提取方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109521470A (zh) | 2019-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| CN108802813B (zh) | 一种多分量地震资料偏移成像方法及系统 | |
| CN109490963B (zh) | 裂缝储层岩石物理建模方法及系统 | |
| CN108139498A (zh) | 具有振幅保持的fwi模型域角度叠加 | |
| CN106556861B (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
| CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
| US11733413B2 (en) | Method and system for super resolution least-squares reverse time migration | |
| CN109143351B (zh) | 叠前各向异性特征参数反演方法及计算机可读存储介质 | |
| CN113031068B (zh) | 一种基于反射系数精确式的基追踪叠前地震反演方法 | |
| CN104975851B (zh) | 用于振幅随炮检距变化道集分析的油藏模型优化方法 | |
| CN113740901A (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
| CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
| CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
| CN106125133B (zh) | 一种基于气云区约束下的精细速度建模方法 | |
| CN109143362B (zh) | 基于共散射角道集的散射波分离方法 | |
| CN107479091B (zh) | 一种提取逆时偏移角道集的方法 | |
| US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
| CN109188522B (zh) | 速度场构建方法及装置 | |
| CN110488385A (zh) | 三维断溶体模型的制作方法 | |
| CN108508486A (zh) | 一种基于构造演化地层应变的裂缝预测分析方法 | |
| Yang et al. | Separation of split shear waves based on a hodogram analysis of HTI media | |
| Wang et al. | Rock Fracture Monitoring Based on High‐Precision Microseismic Event Location Using 3D Multiscale Waveform Inversion | |
| CN109581521B (zh) | Tti各向异性的局部层析方法及系统 | |
| CN104316959B (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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190326 Assignee: CHENGDU JINGSHI PETROLEUM TECHNOLOGY CO.,LTD. Assignor: SOUTHWEST PETROLEUM University Contract record no.: X2023980035129 Denomination of invention: Method for analyzing the influence of tectonics on seismic inversion fracture density Granted publication date: 20200414 License type: Exclusive License Record date: 20230504 |