CN109633745B - 一种三维构造图的制图方法及装置 - Google Patents
一种三维构造图的制图方法及装置 Download PDFInfo
- Publication number
- CN109633745B CN109633745B CN201710929513.4A CN201710929513A CN109633745B CN 109633745 B CN109633745 B CN 109633745B CN 201710929513 A CN201710929513 A CN 201710929513A CN 109633745 B CN109633745 B CN 109633745B
- Authority
- CN
- China
- Prior art keywords
- geological
- data
- dimensional
- map
- establishing
- 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 60
- 238000010586 diagram Methods 0.000 title claims description 29
- 238000010276 construction Methods 0.000 title claims description 5
- 238000005553 drilling Methods 0.000 claims abstract description 35
- 238000012545 processing Methods 0.000 claims description 30
- 238000013507 mapping Methods 0.000 claims description 22
- 238000009499 grossing Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims description 7
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims description 7
- 238000011084 recovery Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 abstract description 15
- 239000003208 petroleum Substances 0.000 abstract description 5
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 4
- 229930195733 hydrocarbon Natural products 0.000 abstract description 4
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 230000002349 favourable effect Effects 0.000 description 14
- 238000013508 migration Methods 0.000 description 9
- 230000005012 migration Effects 0.000 description 9
- 238000009825 accumulation Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response 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/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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种三维构造图的制图方法及装置,将待识别区的等时地层格架中地质层位的网格数据进行处理,分别建立了等值线图、流经线折线图,结合根据钻井资料建立的钻井位置图,将这些图层叠加显示,得到带流经线的三维古构造图,为石油地质工作者直观地三维立体显示某地质层位的古构造和现今构造空间特征以及古流体流经线、汇集区与方向,满足不同研究目的的需要,指导待识别区油气勘探。本发明可以用于实现含油气盆地某地质层位在某地质时期的古构造和流经线的三维立体显示,也可以用于实现含油气盆地某地质层位的现今构造和流经线的三维立体显示,指导油气勘探。
Description
技术领域
本发明属于油气勘探技术领域,具体涉及一种三维构造图的制图方法及装置。
背景技术
在油气勘探技术领域,油气运移、富集规律分析技术是预测有利圈闭及其含油气性和油气充满程度的关键技术,是精细勘探和提高钻井成功率的重要技术保证。根据“高效勘探、效益开发”对勘探技术的要求,该技术的核心任务之一是,必须更精准地刻画构造,同时更可靠地表征油气运移、富集规律。
王兴军在《低幅度构造变速成图技术》一文中,针对塔里木盆地低幅度构造区构造成图的难点,提出了相应的速度建场和变速成图方法:通过沿层速度分析技术提高叠加速度的横向分辨率;利用叠加速度或原始叠前道集作为输入反演层速度,前者不受水平层状介质及射线垂直入射等前提条件的限制,对复杂地层产状具有更好的适应性,后者克服了叠加速度精度的影响,其精度更高;利用叠加速度闭合差校正和均方根速度闭合差校正,得到相对精确的速度场,用于三维时深转换;通过研究深度误差与时间、速度、埋深等参数的相关性,找到深度误差分布规律,进而对构造图进行两次校正。与传统的直射线成图技术相比,此法可明显提高低幅度构造区深度域构造图的精度。
公布号为CN103885083、名称为《古构造图成图方法及装置》的中国专利公开了一种古构造图成图方法及装置,通过结合井校数据和多个地震属性,估算了井校点的深度误差和地震属性之间的响应关系,进而利用地震属性对于无井区域的深度进行修正,从而能够在井控数量较少的情况下生成更准确可靠的古构造图。
上述这些构造图成图方法虽然提高了对构造刻画和解释的精度,较可靠地揭示了油气富集有利区,但在用来分析和描述油气运移、富集规律时,均不能有效表征油气流体的运聚路径和有利聚集区与方向,制约了对有利圈闭及其含油气性和油气充满程度的较可靠预测。
发明内容
本发明的目的是提供一种三维构造图的制图方法及装置,用于解决现有技术不能有效表征油气流体的运聚路径和有利聚集区的问题。
为解决上述技术问题,本发明提出一种三维构造图的制图方法,包括以下解决方案:
方法方案一,包括以下步骤:
1)根据待识别区的钻井资料,获取井位坐标数据;根据待识别区的井震资料,建立待识别区的等时地层格架;
2)对所述等时地层格架中地质层位的数据采用网格化算法进行处理,得到地质层位的网格数据,将地质层位中任意两个地质层位的网格数据相减,得到两个地质层位之间地质单元的残余厚度网格数据;
3)对所述残余厚度网格数据,计算等值线的正交轨迹线,得到流经线数据;
4)根据所述井位坐标数据,建立钻井位置图;根据所述残余厚度网格数据,建立等值线图;根据所述流经线数据,建立流经线折线图,将钻井位置图、等值线图和流经线折线图进行三维立体叠合显示,得到两个地质层位中地质时代较老的地质层位的带流经线的三维古构造图。
方法方案二,在方法方案一的基础上,对所述两个地质层位中的任意一个地质层位的网格数据进行平滑处理、边界处理后,建立所述两个地质层位中地质时代较老的地质层位的构造等值线图和构造面三维立体图,该构造等值线图、构造面三维立体图和所述钻井位置图、流经线折线图进行三维立体叠合显示,得到所述两个地质层位中地质时代较老的地质层位的带流经线的三维现今构造图。
方法方案三,在方法方案一的基础上,步骤2)中所述处理包括真厚度恢复处理、平滑处理、边界处理。
方法方案四、五,分别在方法方案二、三的基础上,所述平滑处理为指定一个0~1之间的系数W,利用此系数用来平衡待光滑点Zf与其周围控制点的值Zc,进而得到光滑后的节点的值Z,其计算公式为Z=(Zf+W×Zc)/(1+W)。
方法方案六,在方法方案一的基础上,所述等时地层格架的建立步骤如下:
S1、根据待识别区的井震资料,进行骨干剖面地层划分,建立待识别区的骨干剖面等时地层格架;
S2、在骨干剖面等时地层格架的约束下,进行构造-地层联动解释,完成每一地质层位的全区追踪闭合,建立待识别区的等时地层格架。
方法方案七,在方法方案六的基础上,划分所述骨干剖面地层时,采取避开构造复杂地带,和分别连接设定个数的已钻井和取心井的原则。
方法方案八,在方法方案一的基础上,步骤2)中采用网格化算法进行处理为:
P1、对所述等时地层格架中每一地质层位的数据,将围绕待计算节点为中心、以搜索半径为半径的邻域平均分为八个扇区,然后在所分的每一个扇区中,找出距离中心点近的点;利用最小二乘法计算中心点的值的步骤为:当距离中心点距离最近的点的距离小于百分之一网格间距时,用所述距离最近的点的值代替中心点的值;当中心点邻域内没有数据点时,把所述中心点赋值为无效值;
P2、对所述等时地层格架中每一地质层位的数据,利用最小二乘法计算中心点的值。
为解决上述技术问题,本发明还提出一种三维构造图的制图装置,包括以下解决方案:
装置方案一,包括处理器,用于执行实现以下方法的指令:
1)根据待识别区的钻井资料,获取井位坐标数据;根据待识别区的井震资料,建立待识别区的等时地层格架;
2)对所述等时地层格架中地质层位的数据采用网格化算法进行处理,得到地质层位的网格数据,将地质层位中任意两个地质层位的网格数据相减,得到两个地质层位之间地质单元的残余厚度网格数据;
3)对所述残余厚度网格数据,计算等值线的正交轨迹线,得到流经线数据;
4)根据所述井位坐标数据,建立钻井位置图;根据所述残余厚度网格数据,建立等值线图;根据所述流经线数据,建立流经线折线图,将钻井位置图、等值线图和流经线折线图进行三维立体叠合显示,得到两个地质层位中地质时代较老的地质层位的带流经线的三维古构造图。
装置方案二,在装置方案一的基础上,对所述两个地质层位中的任意一个地质层位的网格数据进行平滑处理、边界处理后,建立所述两个地质层位中地质时代较老的地质层位的构造等值线图和构造面三维立体图,该构造等值线图、构造面三维立体图和所述钻井位置图、流经线折线图进行三维立体叠合显示,得到所述两个地质层位中地质时代较老的地质层位的带流经线的三维现今构造图。
装置方案三,在装置方案一的基础上,步骤2)中所述处理包括真厚度恢复处理、平滑处理、边界处理。
装置方案四、五,分别在装置方案二、三的基础上,所述平滑处理为指定一个0~1之间的系数W,利用此系数用来平衡待光滑点Zf与其周围控制点的值Zc,进而得到光滑后的节点的值Z,其计算公式为Z=(Zf+W×Zc)/(1+W)。
装置方案六,在装置方案一的基础上,所述等时地层格架的建立步骤如下:
S1、根据待识别区的井震资料,进行骨干剖面地层划分,建立待识别区的骨干剖面等时地层格架;
S2、在骨干剖面等时地层格架的约束下,进行构造-地层联动解释,完成每一地质层位的全区追踪闭合,建立待识别区的等时地层格架。
装置方案七,在装置方案六的基础上,划分所述骨干剖面地层时,采取避开构造复杂地带,和分别连接设定个数的已钻井和取心井的原则。
装置方案八,在装置方案一的基础上,步骤2)中采用网格化算法进行处理为:
P1、对所述等时地层格架中每一地质层位的数据,将围绕待计算节点为中心、以搜索半径为半径的邻域平均分为八个扇区,然后在所分的每一个扇区中,找出距离中心点近的点;利用最小二乘法计算中心点的值的步骤为:当距离中心点距离最近的点的距离小于百分之一网格间距时,用所述距离最近的点的值代替中心点的值;当中心点邻域内没有数据点时,把所述中心点赋值为无效值;
P2、对所述等时地层格架中每一地质层位的数据,利用最小二乘法计算中心点的值。
本发明的有益效果是:
本发明将待识别区的等时地层格架中地质层位的网格数据进行处理,分别建立了等值线图、流经线折线图,结合根据钻井资料建立的钻井位置图,将这些图层叠加显示,得到带流经线的三维古构造图,表达了其中一个地质时代较新层位相应地质时期、流体沿地质时代较老层位的运聚路径、聚集区与方向,以及流体流经线与该地质层位其时构造的空间关系,实现了流体流经线与构造的三维立体叠合成图。本发明根据勘探生产的研究目的需要,所生成图件可以更直观清晰地表达在某个地质时期如大量生排烃期油气流体沿某地质层位的运聚路径、有利聚集区和方向,提高对沿途有利圈闭及其含油气性的预测精度,提升油气钻探的成功率。
进一步,本发明针对其中地质时代较老层位的网格数据进行处理,建立该地质层位的构造等值线图和构造面三维立体图,结合钻井位置图和流经线折线图,进行三维立体叠合显示,得到该地质层位的带流经线的三维现今构造图,有效提升了钻井的成功率,为石油地质工作者更直观明了地刻画、描述油气流体的运聚路径和有利聚集区与方向提供一种有效可行的技术途径,能够满足含油气盆地高效勘探、效益开发的需要。
附图说明
图1是本发明技术方案流程框图;
图2是A盆地某单井合成地震记录标定地震地质层位图;
图3是A盆地中沿某条井震联合剖面的地震层序格架图;
图4是A盆地中沿某条地震剖面的等时地层格架图;
图5是A盆地某地质层位的三维古构造图;
图6是A盆地某地质层位的带流经线三维古构造图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
在目前常用的地质成图软件中,有的软件,如地震地质解释软件LandMark中的ZmapPlus地质绘图模块,可以很好地绘制地质平面图;有的软件,如Surfer,可以很好地实现地质图三维立体显示和生成矢量图,它们在地质成图方面有着各自的独特优势。但由于不同成图软件之间数据格式不兼容,从某一软件导出的成图数据,需要进行数据格式转换后,才能应用于另一软件进行地质成图,达到追求的地质图件效果。
随着油气勘探的深入,石油地质工作者持续加强精细勘探,对于石油地质图件所表达的地质信息和内容要求越来越丰富和直观,精度要求越来越高,单一软件成图已不能完全满足当前精细勘探地质成图的需要。充分发挥不同软件地质成图的各自优势,将具有快捷强大的地质平面图绘制功能的软件(如地震地质解释软件LandMark中的ZmapPlus地质绘图模块)与具有强大的三维立体地质成图功能的软件(如Surfer)有机地结合,通过数据格式转换,按照追求的地质图件效果,进行地质成图,可为石油地质工作者更直观明了地刻画、描述油气流体的运聚路径和有利聚集区以及提高有利圈闭预测成功率提供一种有效可行的技术途径。
因此,本发明利用地震地质解释软件LandMark中的ZmapPlus地质绘图模块和Surfer软件,以中国石油化工股份有限公司A盆地为例,结合附图,对本发明作进一步描述。
本发明的一种三维构造图制图方法的实施例:
1、A盆地资料收集与处理:
1.1区域地质资料收集;
1.2地震、钻井和岩心资料收集,要求覆盖整个A盆地;
1.3井震资料降噪保真处理;
1.4井震资料品质评价,地震资料品质满足波组特征清晰、断层断点清楚,钻井资料满足校深准确、岩心归位正确。
2、对由步骤1得到的达到资料品质要求的A盆地井震资料,在地层学理论指导下进行地层划分,建立A盆地等时地层格架:
2.1依据井震资料中的地层界面标志,识别不整合面,划分测井层序和地震层序;
2.2利用各单井剖面资料、井旁地震子波做合成地震记录,标定步骤2.1划分地震层序的地震地质层位,图2显示了A盆地某单井合成地震记录对井旁地震地质层位的标定结果;
2.3利用单井剖面资料绘制连井剖面和井震联合剖面,将步骤2.1划分的测井层序与井旁地震层序进行一致性匹配,建立起井震联合剖面地震层序,图3显示了A盆地中沿某条井震联合剖面的地震层序格架;
2.4在步骤1.4的地震资料中,按照一要穿越A盆地的主要部位、尽量连接较多已钻井、最好能通过一些重要取心井,二要尽可能避开构造复杂地带的原则,选取十字网格分布的若干骨干剖面,形成骨干剖面网,在步骤2.3确定的井震联合剖面地震层序的约束下,对步骤2.2标定的地震地质层位进行追踪闭合,建立A盆地骨干剖面网等时地层格架,图4显示了A盆地骨干剖面网中某条地震剖面的等时地层格架;
2.5在由步骤2.4得到的A盆地骨干剖面网等时地层格架约束下,对步骤2.2标定的某一地震地质层位,进行A盆地地震资料的全区构造-地层联动解释,完成该地质层位的全区追踪闭合和断层断点的闭合;
2.6重复步骤2.5,完成步骤2.2标定的每一地震地质层位的全区追踪闭合和断层断点的闭合,建立A盆地等时地层格架,图4同时显示了A盆地中沿某条地震剖面的等时地层格架。
3、对由步骤2.6得到的A盆地等时地层格架中的某两个地质层位,例如图4中的Sb1和Sb5,以(x,y,z)的数据格式进行ASCII输出,进行数据网格计算和网格数据处理:
3.1数据的网格计算包括两部分:数据的收集和数据的计算;
3.2步骤3.1所述的数据收集是,针对步骤3所述的Sb1地质层位的输出数据,将围绕待计算节点为中心、以搜索半径为半径的邻域平均分为八个扇区,然后在所分的每一个扇区中,找出距离中心点近的点;
3.3步骤3.1所述的数据计算,是针对步骤3所述的Sb1地质层位的输出数据,利用最小二乘法的算法计算中心点的值;如果距离中心点距离最近的点的距离小于百分之一网格间距时,就用该点的值代替中心点的值;如果中心点邻域内没有数据点,则把该中心点赋值为无效值;例如,由步骤3获得某一地质层位的n组数据(xi,yi,zi,ri),i=1,…,n;其中xi、yi、zi、ri分别为x坐标、y坐标、数据点值及到中心点的距离;①如果n=0,则中心点数据值为无效值ZNON;②如果n不等于0,若有最小距离点小于百分之一网格间距时,中心点数据值等于该点的数据值;否则,a、如果n<3,利用反距离法计算中心点数据值;b、如果n≥3,令wi=1/ri,用最小二乘法确定一个平面方程Ax+By+C=z,其中A、B、C为待定系数,该方程使得最小,求得A、B、C的值,进而得到中心点数据值z;
3.4经由步骤3.2和步骤3.3,完成Sb1地质层位的数据网格化;
3.5针对步骤3所述的Sb5地质层位的输出数据,重复步骤3.2和步骤3.3,完成Sb5地质层位的数据网格化;
3.6对由步骤3.4和步骤3.5得到的两个地质层位的网格数据进行相减计算,即Sb1地质层位的网格数据减Sb5地质层位的网格数据,得到Sb1地质层位和Sb5地质层位之间地质单元的残余厚度数据网格;
3.7对步骤3.6所述的地质单元的残余厚度数据网格进行真厚度恢复处理;
3.8对由步骤3.7得到的真厚度恢复处理后的网格数据进行光滑处理:数据光滑处理是指定一个0~1之间的系数W,利用此系数用来平衡待光滑点Zf与其周围控制点的值Zc,进而得到光滑后的节点的值Z,其计算公式为Z=(Zf+W×Zc)/(1+W);例如,待光滑点Zf周围有8个控制点,值分别为976、960、1005、912、1067、986、1004、994,则Zc=(976+960+1005+912+1067+986+1004+994)/8=988;Zf为1000,W为0.25,则Z=(Zf+W×Zc)/(1+W)=(1000+0.25×988)/(1+0.25)=997.6;
3.9对由步骤3.8得到的光滑后的网格数据进行边界处理,该数据能够按ASCII输出,建立与Surfer软件之间的数据接口。
4、对由步骤3.9得到的步骤3.6所述的地质单元的网格数据进行流经线数据计算:
4.1步骤4所述的流经线数据计算即通过计算等值线的正交线得到流经线数据;
4.2由步骤4.1得到的流径线数据按ASCII输出,建立与Surfer软件之间的数据接口。
5、根据步骤1.2钻井资料中的井位坐标数据、步骤3.9得到的网格数据和步骤4.2得到的流径线数据,应用Surfer软件进行三维立体叠合成图:
5.1根据步骤1.2钻井资料中的井位坐标数据,建立钻井位置图图层,用Surfer软件的Post map显示;
5.2针对步骤3.6所述的地质单元,根据步骤3.9得到的网格数据,建立等值线图层,用Surfer软件的Contour map显示;根据步骤4.2得到的流径线数据,建立流径线折线图图层,用Surfer软件的Base map显示;根据步骤3.9得到的网格数据,建立构造面三维立体图图层,用Surfer软件的3D Surface map显示;图5为A盆地中Sb1地质层位的构造等值线图和构造面三维立体图,属古构造图;
5.3对步骤3.4所述的Sb1地质层位的网格数据,先后进行光滑处理和边界处理,根据处理后得到的网格数据,建立Sb1地质层位的构造等值线图图层和构造面三维立体图图层,分别用Surfer软件的Contour map和3D Surface map显示;
5.4针对步骤3所述的Sb1和Sb5地质层位,将步骤5.1和步骤5.2得到的各图层,通过Surfer软件的Overlay map进行三维立体叠合显示,生成成果图件,得到Sb1地质层位对应于Sb5地质层位相应地质时期的带流经线三维构造图,属带流经线三维古构造图,图6为A盆地中Sb1地质层位的带流经线三维古构造图,表达了在Sb5地质层位相应地质时期,流体沿Sb1地质层位的运聚路径、聚集区与方向以及流体流经线与该地质层位其时构造的空间关系;
5.5将步骤5.1得到的钻井位置图图层和步骤5.2得到的流径线折线图图层,与步骤5.3得到的Sb1地质层位的构造等值线图图层和构造面三维立体图图层,通过Surfer软件的Overlay map进行三维立体叠合显示,生成成果图件,得到步骤3所述的Sb1地质层位的带流经线三维构造图,属带流经线三维现今构造图。
针对步骤3所述的Sb1和Sb5地质层位,由步骤5.4和步骤5.5得到的Sb1地质层位的带流经线三维构造图,可以为石油地质工作者直观地三维立体显示Sb1地质层位在Sb5地质层位对应地质时期的古构造和现今构造的空间特征,以及在Sb5地质层位对应地质时期古流体沿Sb1地质层位的流经线、汇集区与方向,满足不同勘探研究目的的需要,指导A盆地油气勘探。
本发明是在建立起待识别区等时地层格架基础上,实现流体流经线与构造的三维立体叠合成图,根据勘探生产的研究目的需要,所生成图件可以更直观清晰地表达在某个地质时期如大量生排烃期油气流体沿某地质层位的运聚路径、有利聚集区和方向,提高对沿途有利圈闭及其含油气性的预测精度,提升油气钻探的成功率,为石油地质工作者更直观明了地刻画、描述油气流体的运聚路径和有利聚集区与方向提供一种有效可行的技术途径,能够满足含油气盆地高效勘探、效益开发的需要。
本发明还提出了一种三维构造图的制图装置,包括处理器,用于执行实现以下方法的指令:
1)根据待识别区的钻井资料,获取井位坐标数据;根据待识别区的井震资料,建立待识别区的等时地层格架;
2)对所述等时地层格架中地质层位的数据采用网格化算法进行处理,得到地质层位的网格数据,将地质层位中任意两个地质层位的网格数据相减,得到两个地质层位之间地质单元的残余厚度网格数据;
3)对所述残余厚度网格数据,计算等值线的正交轨迹线,得到流经线数据;
4)根据所述井位坐标数据,建立钻井位置图;根据所述残余厚度网格数据,建立等值线图;根据所述流经线数据,建立流经线折线图,将钻井位置图、等值线图和流经线折线图进行三维立体叠合显示,得到两个地质层位中地质时代较老层位的带流经线的三维古构造图。
上述实施例中所指的构造图制图装置,实际上是基于本发明方法流程的一种计算机解决方案,即一种软件构架,可以应用到处理器中,上述装置即为与方法流程相对应的处理进程。由于对上述方法的介绍已经足够清楚完整,故不再详细进行描述。
以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (14)
1.一种三维构造图的制图方法,其特征在于,包括以下步骤:
1)根据待识别区的钻井资料,获取井位坐标数据;根据待识别区的井震资料,建立待识别区的等时地层格架;
2)对所述等时地层格架中地质层位的数据采用网格化算法进行处理,得到地质层位的网格数据,将地质层位中任意两个地质层位的网格数据相减,得到两个地质层位之间地质单元的残余厚度网格数据;
3)对所述残余厚度网格数据,计算等值线的正交轨迹线,得到流经线数据;
4)根据所述井位坐标数据,建立钻井位置图;根据所述残余厚度网格数据,建立等值线图;根据所述流经线数据,建立流经线折线图,将钻井位置图、等值线图和流经线折线图进行三维立体叠合显示,得到两个地质层位中地质时代较老的地质层位的带流经线的三维古构造图。
2.根据权利要求1所述的三维构造图的制图方法,其特征在于,对所述两个地质层位中的任意一个地质层位的网格数据进行平滑处理、边界处理后,建立所述两个地质层位中地质时代较老的地质层位的构造等值线图和构造面三维立体图,该构造等值线图、构造面三维立体图和所述钻井位置图、流经线折线图进行三维立体叠合显示,得到所述两个地质层位中地质时代较老的地质层位的带流经线的三维现今构造图。
3.根据权利要求1所述的三维构造图的制图方法,其特征在于,步骤2)中采用网格化算法处理后,进行真厚度恢复处理、平滑处理、边界处理。
4.根据权利要求2或3所述的三维构造图的制图方法,其特征在于,所述平滑处理为指定一个0~1之间的系数W,利用此系数用来平衡待光滑点Zf与其周围控制点的值Zc,进而得到光滑后的节点的值Z,其计算公式为Z=(Zf+W×Zc)/(1+W)。
5.根据权利要求1所述的三维构造图的制图方法,其特征在于,所述等时地层格架的建立步骤如下:
S1、根据待识别区的井震资料,进行骨干剖面地层划分,建立待识别区的骨干剖面等时地层格架;
S2、在骨干剖面等时地层格架的约束下,进行构造-地层联动解释,完成每一地质层位的全区追踪闭合,建立待识别区的等时地层格架。
6.根据权利要求5所述的三维构造图的制图方法,其特征在于,划分所述骨干剖面地层时,采取避开构造复杂地带,和分别连接设定个数的已钻井和取心井的原则。
7.根据权利要求1所述的三维构造图的制图方法,其特征在于,步骤2)中采用网格化算法进行处理为:
P1、对所述等时地层格架中每一地质层位的数据,将围绕待计算节点为中心、以搜索半径为半径的邻域平均分为八个扇区,然后在所分的每一个扇区中,找出距离中心点近的点;利用最小二乘法计算中心点的值的步骤为:当距离中心点距离最近的点的距离小于百分之一网格间距时,用所述距离最近的点的值代替中心点的值;当中心点邻域内没有数据点时,把所述中心点赋值为无效值;
P2、对所述等时地层格架中每一地质层位的数据,利用最小二乘法计算中心点的值。
8.一种三维构造图的制图装置,其特征在于,包括处理器,用于执行实现以下方法的指令:
1)根据待识别区的钻井资料,获取井位坐标数据;根据待识别区的井震资料,建立待识别区的等时地层格架;
2)对所述等时地层格架中地质层位的数据采用网格化算法进行处理,得到地质层位的网格数据,将地质层位中任意两个地质层位的网格数据相减,得到两个地质层位之间地质单元的残余厚度网格数据;
3)对所述残余厚度网格数据,计算等值线的正交轨迹线,得到流经线数据;
4)根据所述井位坐标数据,建立钻井位置图;根据所述残余厚度网格数据,建立等值线图;根据所述流经线数据,建立流经线折线图,将钻井位置图、等值线图和流经线折线图进行三维立体叠合显示,得到两个地质层位中地质时代较老的地质层位的带流经线的三维古构造图。
9.根据权利要求8所述的三维构造图的制图装置,其特征在于,对所述两个地质层位中的任意一个地质层位的网格数据进行平滑处理、边界处理后,建立所述两个地质层位中地质时代较老的地质层位的构造等值线图和构造面三维立体图,该构造等值线图、构造面三维立体图和所述钻井位置图、流经线折线图进行三维立体叠合显示,得到所述两个地质层位中地质时代较老的地质层位的带流经线的三维现今构造图。
10.根据权利要求8所述的三维构造图的制图装置,其特征在于,步骤2)中采用网格化算法处理后,进行真厚度恢复处理、平滑处理、边界处理。
11.根据权利要求9或10所述的三维构造图的制图装置,其特征在于,所述平滑处理为指定一个0~1之间的系数W,利用此系数用来平衡待光滑点Zf与其周围控制点的值Zc,进而得到光滑后的节点的值Z,其计算公式为Z=(Zf+W×Zc)/(1+W)。
12.根据权利要求8所述的三维构造图的制图装置,其特征在于,所述等时地层格架的建立步骤如下:
S1、根据待识别区的井震资料,进行骨干剖面地层划分,建立待识别区的骨干剖面等时地层格架;
S2、在骨干剖面等时地层格架的约束下,进行构造-地层联动解释,完成每一地质层位的全区追踪闭合,建立待识别区的等时地层格架。
13.根据权利要求12所述的三维构造图的制图装置,其特征在于,划分所述骨干剖面地层时,采取避开构造复杂地带,和分别连接设定个数的已钻井和取心井的原则。
14.根据权利要求8所述的三维构造图的制图装置,其特征在于,步骤2)中采用网格化算法进行处理为:
P1、对所述等时地层格架中每一地质层位的数据,将围绕待计算节点为中心、以搜索半径为半径的邻域平均分为八个扇区,然后在所分的每一个扇区中,找出距离中心点近的点;利用最小二乘法计算中心点的值的步骤为:当距离中心点距离最近的点的距离小于百分之一网格间距时,用所述距离最近的点的值代替中心点的值;当中心点邻域内没有数据点时,把所述中心点赋值为无效值;
P2、对所述等时地层格架中每一地质层位的数据,利用最小二乘法计算中心点的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710929513.4A CN109633745B (zh) | 2017-10-09 | 2017-10-09 | 一种三维构造图的制图方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710929513.4A CN109633745B (zh) | 2017-10-09 | 2017-10-09 | 一种三维构造图的制图方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633745A CN109633745A (zh) | 2019-04-16 |
CN109633745B true CN109633745B (zh) | 2020-09-11 |
Family
ID=66051242
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710929513.4A Active CN109633745B (zh) | 2017-10-09 | 2017-10-09 | 一种三维构造图的制图方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633745B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111751880B (zh) * | 2019-03-29 | 2023-07-25 | 中国石油天然气集团有限公司 | 井震数据联合显示方法及装置 |
CN110259430A (zh) * | 2019-06-27 | 2019-09-20 | 北京中油瑞飞信息技术有限责任公司 | 一种部署井位的方法及装置 |
CN110821496B (zh) * | 2019-10-17 | 2021-06-29 | 中国石油天然气集团有限公司 | 有机质页岩相模式建立方法及有机质页岩评价方法 |
CN113945993A (zh) * | 2020-07-16 | 2022-01-18 | 中国石油化工股份有限公司 | 一种滩相储层预测方法及装置 |
CN112415600B (zh) * | 2020-11-02 | 2024-04-30 | 中国石油天然气集团有限公司 | 一种即时互动表层构造恢复方法及装置 |
CN115393528B (zh) * | 2022-09-20 | 2023-04-18 | 中国地震局地球物理研究所 | 一种古地震探槽的三维建模方法及其系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9104697B2 (en) * | 2012-07-24 | 2015-08-11 | Navico Holdings As | Aquatic geographic information system |
CN103206207B (zh) * | 2013-03-01 | 2015-06-17 | 陕西延长石油(集团)有限责任公司研究院 | 基于生烃期古构造的油藏有利区的识别方法 |
CN105425316A (zh) * | 2015-10-28 | 2016-03-23 | 中国石油化工股份有限公司 | 一种基于印模法的古地貌恢复方法 |
CN106837319B (zh) * | 2016-12-27 | 2019-10-11 | 中国石油天然气股份有限公司 | 一种密井网条件下的地下单一古河道识别方法 |
-
2017
- 2017-10-09 CN CN201710929513.4A patent/CN109633745B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109633745A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633745B (zh) | 一种三维构造图的制图方法及装置 | |
Qi et al. | Volumetric aberrancy to map subtle faults and flexures | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
EA011781B1 (ru) | Способ геологического моделирования с помощью основанного на гидродинамике гридинга (гидросеток) | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
Benesh et al. | Geometry, kinematics, and displacement characteristics of tear-fault systems: An example from the deep-water Niger Delta | |
CN108005646B (zh) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 | |
NO315445B1 (no) | Fremgangsmåte for fremstilling av minst ±n borehullsbane i et rom-tid-domene | |
Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
CN107942379A (zh) | 一种提高复杂断块速度模型精度的方法 | |
CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
CN106569272A (zh) | 一种基于数据性质空间升维的地震属性融合方法 | |
Barthwal et al. | Passive seismic tomography using recorded microseismicity: Application to mining-induced seismicity | |
Shaban et al. | Comparison between curvature and 3D strain analysis methods for fracture predicting in the Gachsaran oil field (Iran) | |
US20210396897A1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
CN104749625A (zh) | 一种基于正则化技术的地震数据倾角估计方法及装置 | |
Al Rassas et al. | Application of 3d reservoir geological model on es1 formation, block nv32, shenvsi oilfield, China | |
Chen et al. | Geological risk mapping and prospect evaluation using multivariate and Bayesian statistical methods, western Sverdrup Basin of Canada | |
Gupta et al. | Seismic attribute driven integrated characterization of the Woodford Shale in west-central Oklahoma | |
CN109283574B (zh) | 低频模型构建方法及计算机可读存储介质 | |
CN115880455A (zh) | 基于深度学习的三维智能插值方法 | |
Shao et al. | 3D Geological Modeling under Extremely Complex Geological Conditions. | |
RU2681250C1 (ru) | Способ квазитрехмерного моделирования эффективных газо- и нефтенасыщенных толщин залежей углеводородов | |
Solum et al. | Increasing interpreter capability in structurally complex settings through combined fieldwork, interpretation, and geocellular modeling |
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 |