CN110488385A - 三维断溶体模型的制作方法 - Google Patents
三维断溶体模型的制作方法 Download PDFInfo
- Publication number
- CN110488385A CN110488385A CN201910894466.3A CN201910894466A CN110488385A CN 110488385 A CN110488385 A CN 110488385A CN 201910894466 A CN201910894466 A CN 201910894466A CN 110488385 A CN110488385 A CN 110488385A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- ray
- face
- model
- disconnected
- 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.)
- Pending
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 9
- 235000019738 Limestone Nutrition 0.000 claims description 6
- 239000006028 limestone Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005260 corrosion Methods 0.000 claims description 3
- 230000007797 corrosion Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 2
- 239000000155 melt Substances 0.000 abstract 1
- 238000011161 development Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- 238000003325 tomography Methods 0.000 description 4
- 208000010392 Bone Fractures Diseases 0.000 description 2
- 206010017076 Fracture Diseases 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012800 visualization 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
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
三维断溶体模型的制作方法,这种方法的核心是用户选择的特定波形才能在合成地震图中建模,该方法是先对断溶体的各种地质体的速度,厚度,以及比例尺进行采集和统计,再进行断溶体二维平面速度模型的绘制,等间距绘制完成几个二维的平面速度模型在二维面上找一个控制点,由这点建立一条射线,与下一个剖面建立映射的关系;再用同样的方法与下一面建立同样的映射,利用样条插值的方法,对面与面之间的空白部分进行填充,从而形成三维模型,对于断熔体的内外部构造描述有了很大提高,提高了对断熔体精确描述的准确性。
Description
技术领域
本发明涉及涉及石油,页岩气储层地质学领域的一种对断裂带储层的建模方法,具体技术为基于小波对垂直反射率对数进行卷积的建模方法,可以有效提高对油田开发阶段断溶体的结构的识别与预测的精确性以及对断溶体内部表征的准确表现。
背景技术
目前,随着油田的勘探开发的深入及保持油田产量平稳的要求,且现在的勘探开发技术也得到了长足的进步,油田勘探开发的目标已逐渐从以往的背斜油气藏、构造油气藏过渡到岩性油气藏、隐蔽油气藏勘探,而断块型油气藏就是一类重要的隐蔽油气藏,其油气资源和生产能力,所占比重很大,勘探前景比较广阔,对于此类油气藏,断层控制着其形成及开发,因此识别断层的类型及空间展布就具有很重要的作用,在这其中,对断层及断层破碎带的认识程度更是影响勘探效果的重要影响因素。
断层破碎带外带以裂缝发育为特征,通常以单一的组裂缝为主, 裂缝发育程度较内带明显降低,地层连续性好,缺少角砾岩与碎裂岩,在地震剖面上通常呈现同相轴连续,振幅减弱,频率降低的特征,破碎带外带与内带、围岩呈渐变关系,常规地震剖面通常较难识别外带的分布,需要精细的裂缝与储层预测方法判识。
2012年Paton提出了现有方法的结果:通过数据量使局部数据统计适应地震表达的变化,该方法将地理体的手动交互式3维编辑与不透明度阈值相结合,在仅有数据驱动技术不足以解决地质目标的区域,对于我们的案例研究,我们提取了使用相似性与方位角属性的交叉图分析获得的干扰地质体,一些与SDZ方向相似的噪声仍然可以抵抗主要的交叉图选择。使用低相似值分布跟踪的SDZ已经被渲染并提取为单个地理体,这些表示已经使用低相干阈值(基于颜色不透明度值) 预定义的SDZ的体积可视化,然后可以用相应的原始地震信号特性或插值方法来覆盖或填充所得的三维地质体的完整属性。
发明内容
本发明的主要目的是提供一种对断溶体的三维结构进行的建模的方法,该方法基于小波对垂直反射率对数进行卷积的建模方法,简单的一维反射方法和简单的一维反射方法之间的良好折衷要求苛刻的FD方法是基于射线的方法,射线追踪方法使用射线理论来计算非均匀介质中沿射线路径的地震波传播时间和振幅,遵循波动方程的高频近似,作为高频或低波长近似,射线理论不能精确地应用于薄层内的结构主波长,该方法也不产生完整的波场;只有用户选择的特定波形才能在合成地震图中建模,此外,射线理论在临界偏移附近不准确,即使可以进行一些修正,对于水平分层介质中的深反射器,反射角通常很小,因此可以安全地使用射线理论,而在比较浅的深度,它可能不合适,然而,射线理论被广泛使用,特别是对于2维分层模型,因为该方法的某些版本相对简单,计算时间快,并且可以从事件旅行中识别合成地震图上所有事件的射线路径,如果不总是信号幅度,该方法对计算行程时间非常有用,一些基于射线的方法已经开发出以便能够处理PSDM图像中的3维横向分辨率和照明这种3维空间卷积方法模拟PSDM点散射响应和反射器),其计算时间接近1维卷积并且比 FD方法更易于访问,这种方法比较容易为地质学家学习,最近已应用于地质模型,包括断层显示出有缺陷的三维断层结果成像。
该方法包括以下步骤:
(1)先对断溶体的各种地质体的速度,厚度,以及比例尺进行采集和统计,再进行断溶体二维平面速度模型的绘制;
(2)等间距绘制完成几个二维的平面速度模型;
(3)在二维面上找一个控制点,由这点建立一条射线,与下一个剖面建立映射的关系(需判断该射线与下一个剖面的每一边交点个数的奇偶性,如果是奇数点的个数则表明另一面与该射线交点在该面以内,偶数则交点在该面以外);再用同样的方法与下一面建立同样的映射。
(4)利用样条插值的方法,对面与面之间的空白部分进行填充,从而形成三维模型。
样条插值的实现算法步骤:
假定有n+1个数据的节点分别为:
(x0,y0),(x1,y1),(x2,y2),……,(xn,yn)
a.计算步长hi=xi+1-xi(i=0,1,…,n-1);
b.将数据节点和指定的首位端点条件带入矩阵方程;
c.解矩阵方程,求得二次微分值mi。(该矩阵为三对角矩阵);
d.计算样条曲线的系数:
ai=yi;
其中i=0,1,2,3.....n-1
e.在每个子区间xi≤x≤xi+1中,创建方程:
gi(x)=ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3.
附图说明:
图1为该方法的实行路线图。
具体实施方式:
1.三维断溶体的模型制作的方法依次包括以下的详细步骤:
步骤一:分别得到上覆地层1和上覆地层2,目的地上覆地层,断裂溶蚀带,溶蚀缝洞系统,致密灰岩,致密灰岩1,致密灰岩2的层速度,根据加德纳公式计算出密度,由速度以及密度计算出波阻抗。
加德纳公式:ρ=0.31*v1/4;
ρ为该层岩体的密度,v为该层岩体的波速度。
波阻抗公式:δ=ρ*v;
δ为该层岩体的波阻抗。
步骤二:涉及的地震参数:子波频率为30Hz。
步骤三:根据示意图,先画出二维的断溶体速度,密度,波阻抗模型图。
步骤四:等间距绘制几个二维的平面速度模型;
步骤五:在二维面上找一个控制点,由这点建立一条射线与下一个剖面建立映射的关系(需判断该射线与下一个剖面每一边的交点个数奇偶性,如果是交点个数为奇数则表明另一面与该射线交点在该面以内,交点为偶数则交点在该面以外)再用同样的方法与下一面建立同样的映射;
步骤六:利用样条插值的方法,对面与面之间的空白部分进行填充,形成三维的模型。
样条插值的实现算法步骤:
假定有n+1个数据的节点分别为:
(x0,y0),(x1,y1),(x2,y2),……,(xn,yn)
a.计算步长hi=xi+1-xi(i=0,1,…,n-1);
b.将数据节点和指定的首位端点条件带入矩阵方程;
c.解矩阵方程,求得二次微分值mi。(该矩阵为三对角矩阵);
d.计算样条曲线的系数:
ai=yi;
其中i=0,1,2,3.....n-1
e.在每个子区间xi≤x≤xi+1中,创建方程:
gi(x)=ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3.
最后,本申请的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围内。
Claims (1)
1.三维断溶体模型的制作方法依次包括以下的详细步骤:
步骤一:分别得到上覆地层1和上覆地层2,目的地上覆地层,断裂溶蚀带,溶蚀缝洞系统,致密灰岩,致密灰岩1,致密灰岩2的速度,根据加德纳公式计算出密度,由速度以及密度计算出波阻抗;
加德纳公式:ρ=0.31*v1/4;
波阻抗公式:δ=ρ*v;
步骤二:涉及的地震参数:子波频率为30Hz;
步骤三:根据示意图,先画出二维的断溶体速度,密度,波阻抗模型图;
步骤四:等间距绘制几个二维的平面速度模型;
步骤五:在二维面上找一个控制点,由这点建立一条射线与下一个剖面建立映射的关系(需判断该射线与下一个剖面每一边的交点个数奇偶性,如果是交点个数为奇数则表明另一面与该射线交点在该面以内,交点为偶数则交点在该面以外)再用同样的方法与下一面建立同样的映射;
步骤六:利用样条插值的方法,对面与面之间的空白部分进行填充,形成三维的模型;
样条插值的实现算法步骤:
假定有n+1个数据的节点分别为:
(x0,y0),(x1,y1),(x2,y2),……,(xn,yn)
a.计算步长hi=xi+1-xi(i=0,1,…,n-1);
b.将数据节点和指定的首位端点条件带入矩阵方程;
c.解矩阵方程,求得二次微分值mi(该矩阵为三对角矩阵);
d.计算样条曲线的系数:
ai=yi;
其中i=0,1,2,3…..n-1
e.在每个子区间xi≤x≤xi+1中,创建方程:
gi(x)=ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910894466.3A CN110488385A (zh) | 2019-09-20 | 2019-09-20 | 三维断溶体模型的制作方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910894466.3A CN110488385A (zh) | 2019-09-20 | 2019-09-20 | 三维断溶体模型的制作方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110488385A true CN110488385A (zh) | 2019-11-22 |
Family
ID=68558878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910894466.3A Pending CN110488385A (zh) | 2019-09-20 | 2019-09-20 | 三维断溶体模型的制作方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110488385A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484909A (zh) * | 2021-09-07 | 2021-10-08 | 西南石油大学 | 一种基于几何网格化和参数分配的缝洞型储层建立方法 |
CN115576007A (zh) * | 2022-11-22 | 2023-01-06 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
RU2815465C1 (ru) * | 2022-11-22 | 2024-03-18 | Саутвест Петролеум Юниверсити | Система полудетерминированного моделирования трещиноватости на основе беспорядочной матрицы |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597496A (zh) * | 2015-01-30 | 2015-05-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种二维地震资料中速度的三维空间归位方法 |
CN106526670A (zh) * | 2016-09-21 | 2017-03-22 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
CN107240152A (zh) * | 2016-03-28 | 2017-10-10 | 中国石油化工股份有限公司 | 三维地质模型的建模方法及系统 |
CN107367757A (zh) * | 2017-07-20 | 2017-11-21 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体的刻画方法 |
CN107390264A (zh) * | 2017-07-20 | 2017-11-24 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体内部结构的表征方法 |
-
2019
- 2019-09-20 CN CN201910894466.3A patent/CN110488385A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597496A (zh) * | 2015-01-30 | 2015-05-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种二维地震资料中速度的三维空间归位方法 |
CN107240152A (zh) * | 2016-03-28 | 2017-10-10 | 中国石油化工股份有限公司 | 三维地质模型的建模方法及系统 |
CN106526670A (zh) * | 2016-09-21 | 2017-03-22 | 中石化石油工程技术服务有限公司 | 一种碎屑岩储层中地震属性砂体空间分布描述及评价的方法 |
CN107367757A (zh) * | 2017-07-20 | 2017-11-21 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体的刻画方法 |
CN107390264A (zh) * | 2017-07-20 | 2017-11-24 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体内部结构的表征方法 |
Non-Patent Citations (1)
Title |
---|
邓飞: "剖面三维地质建模与高斯射线束正演的研究与实现", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484909A (zh) * | 2021-09-07 | 2021-10-08 | 西南石油大学 | 一种基于几何网格化和参数分配的缝洞型储层建立方法 |
CN113484909B (zh) * | 2021-09-07 | 2021-11-19 | 西南石油大学 | 一种基于几何网格化和参数分配的缝洞型储层建立方法 |
CN115576007A (zh) * | 2022-11-22 | 2023-01-06 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
CN115576007B (zh) * | 2022-11-22 | 2023-03-14 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
RU2815465C1 (ru) * | 2022-11-22 | 2024-03-18 | Саутвест Петролеум Юниверсити | Система полудетерминированного моделирования трещиноватости на основе беспорядочной матрицы |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109388817B (zh) | 一种储层裂缝三维建模方法 | |
US9759826B2 (en) | System and method for generating an implicit model of geological horizons | |
Philip et al. | A precise method for determining contoured surfaces | |
US8838425B2 (en) | Generating facies probablity cubes | |
US8731875B2 (en) | System and method for providing data corresponding to physical objects | |
Wilson et al. | From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey | |
EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
US8731873B2 (en) | System and method for providing data corresponding to physical objects | |
CN110031896A (zh) | 基于多点地质统计学先验信息的地震随机反演方法及装置 | |
CN104459768A (zh) | 一种基于可视化的三维空间目标地质体追踪方法 | |
US20110218775A1 (en) | System and Method For Providing Data Corresponding To Physical Objects | |
EA031663B1 (ru) | Способ вычислительной геологии | |
CN104975851B (zh) | 用于振幅随炮检距变化道集分析的油藏模型优化方法 | |
CN110488385A (zh) | 三维断溶体模型的制作方法 | |
CN110007340A (zh) | 基于角度域直接包络反演的盐丘速度密度估计方法 | |
CN107884825A (zh) | 一种基于地震多属性的不确定性建模方法 | |
CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
CN108375794A (zh) | 基于对称观测的vsp缝洞绕射成像技术方法 | |
CN107942373A (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
CN115880455A (zh) | 基于深度学习的三维智能插值方法 | |
Lyu et al. | Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China | |
Fretwell et al. | A new approach to 3-D geological modeling of complex sand injectite reservoirs: The Alba field, United Kingdom central North Sea | |
CN111352154B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191122 |
|
WD01 | Invention patent application deemed withdrawn after publication |