CN114114411B - 三维断层疏导性的定量确定方法及装置 - Google Patents
三维断层疏导性的定量确定方法及装置 Download PDFInfo
- Publication number
- CN114114411B CN114114411B CN202010903256.9A CN202010903256A CN114114411B CN 114114411 B CN114114411 B CN 114114411B CN 202010903256 A CN202010903256 A CN 202010903256A CN 114114411 B CN114114411 B CN 114114411B
- Authority
- CN
- China
- Prior art keywords
- section
- data
- model
- dimensional
- fault
- 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
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 17
- 238000011160 research Methods 0.000 claims abstract description 86
- 238000006073 displacement reaction Methods 0.000 claims abstract description 76
- 239000011435 rock Substances 0.000 claims abstract description 59
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000009826 distribution Methods 0.000 claims abstract description 47
- 238000004590 computer program Methods 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000003370 grooming effect Effects 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 12
- 238000005520 cutting process Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 8
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 230000036244 malformation Effects 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 241001442234 Cosa Species 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000008398 formation water Substances 0.000 claims description 4
- 239000003921 oil Substances 0.000 description 25
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 13
- 238000005553 drilling Methods 0.000 description 12
- 238000007789 sealing Methods 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000011835 investigation Methods 0.000 description 6
- 239000004927 clay Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 244000089409 Erythrina poeppigiana Species 0.000 description 1
- 235000009776 Rathbunia alamosensis Nutrition 0.000 description 1
- 101100510671 Rattus norvegicus Lnpep gene Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
技术领域
本发明涉及油气勘探开发技术领域,特别涉及一种三维断层疏导性的定量确定方法及装置。
背景技术
在断块油气藏勘探开发实践中,断层既可以起到油气的封闭作用也可以起到疏导油气的作用,不封闭的断层即可使油气穿越而过形成油气的侧向运移成藏。
断层的封闭性研究是断层疏导性研究的对立面和关键点,断层封闭性研究一直是困扰油气地质家们的一道难题。断层封闭具双向性,即垂向封闭和侧向封闭。对断层的封闭性进行评价,就是要在熟知研究地区断层特征的基础上,筛选对本区断层封闭或开启影响较大的几个变量,紧密结合正确的断层封闭机理,利用掌握的各种准确的地质和物探资料,选用合适的评价方法进行计算分析,使评价的效果最大程度地反映本区断层封闭的客观实际。许多学者认为断层封闭的主要依据有:①储层与非储层的岩性对置;②断裂带中填充物的物性。早期断层封闭性机制的研究假设沿断层面富粘土层的涂抹能在两盘储层产生低渗的断层泥。最新的研究表明相似的粘土涂抹或断层泥在近岸深水区遮挡油气聚集的断层面是普遍存在的。而近年来,用粘土涂抹方法,尤其是断层泥比率(SGR)法,评价预测厚层碎屑岩层序地区的断层封闭性,国内外许多学者做了许多工作,取得了较大的进展;至此之后,国外学者大多致力于研究断层本身特征,很少见定量评价断层封闭性方法的。
综上,考虑到现有的断层疏导性的确定方法,通常是先分析断层的封堵性,封堵的断层不能疏导油气,不封堵的断层可以疏导油气,形成油气渗流通道,使油气跨越断层运移到断层对盘,在断层分析过程中,有两个方面的问题。一是:现有的封闭性的确定方法,大多是通过SGR断层泥比率法确定封闭性。而SGR法具体实施时,往往只涉及断层一盘的泥岩厚度与断距,并没考虑到断层两盘的岩性影响、断面压实和闭合程度以及泥岩实际涂抹位移变化等因素,导致存在确定封闭性的准确度不高。二是欲评价断层疏导能力,首先得分析封闭性,过程复杂,没有直接的定量化的评价断层疏导能力的参数和计算方法,评价效率低下。计算结果容易出现偏差和误导,造成对断块圈闭封闭性分析不全面而引起勘探失利。
因此,目前缺乏一种断层疏导性的定量确定方法。
发明内容
本发明实施例提出一种三维断层疏导性的定量确定方法,用以定量地分析三维断层疏导性,准确度高,该方法包括:
根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
对所述三维地质模型中的断面进行网格化处理,获得网格模型;
获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。
本发明实施例提出一种三维断层疏导性的定量确定装置,用以定量地分析三维断层疏导性,准确度高,该装置包括:
数据获得模块,用于根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
三维地质模型建立模块,用于根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
网格模型获得模块,用于对所述三维地质模型中的断面进行网格化处理,获得网格模型;
第一计算模块,用于获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
第二计算模块,用于根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
疏导性确定模块,用于在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。
本发明实施例还提出了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述三维断层疏导性的定量确定方法。
本发明实施例还提出了一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述三维断层疏导性的定量确定方法的计算机程序。
在本发明实施例中,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;对所述三维地质模型中的断面进行网格化处理,获得网格模型;获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。在上述过程中,综合考虑了研究区叠后三维地震数据、测井信息,最后建立了三维地质模型,并在获得网格模型后,定量获得了每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,进而定量计算了三维断面疏导系数,使得根据三维断面疏导系数判断研究区三维断层是否疏导的过程更加准确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中三维断层疏导性的定量确定方法的流程图;
图2为本发明实施例中断面上实际泥岩涂抹的位移的示意图;
图3为本发明实施例中三维断层疏导性的定量确定方法的详细流程图;
图4为本发明实施例中三维地质模型的示意图;
图5为本发明实施例中岩性曲线图版的示意图;
图6为本发明实施例中三维断面疏导系数与深度的关系图版;
图7为发明实施例中三维断层疏导性的定量确定装置的示意图;
图8为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
在本说明书的描述中,所使用的“包含”、“包括”、“具有”、“含有”等,均为开放性的用语,即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施,其中的步骤顺序不作限定,可根据需要作适当调整。
图1为本发明实施例中三维断层疏导性的定量确定方法的流程图,如图1所示,该方法包括:
步骤101,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
步骤102,根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
步骤103,对所述三维地质模型中的断面进行网格化处理,获得网格模型;
步骤104,获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
步骤105,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
步骤106,在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。
在本发明实施例提出的方法中,综合考虑了研究区叠后三维地震数据、测井信息,最后建立了三维地质模型,并在获得网格模型后,定量获得了每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,进而定量计算了三维断面疏导系数,使得根据三维断面疏导系数判断研究区三维断层是否疏导的过程更加准确。
具体实施时,在步骤101中,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据,具体包括研究区叠后三维地震数据加载进专业地震解释软件Landmark或者Geoeast中,解释出断面和目的地层,得到断面的总断距、倾角、走向等要素,得到地层的层厚、倾角等要素。现有的方法获得断面数据和层位数据,进而确定断层封闭性来分析疏导性时,大多使用的是钻井资料,要求且钻井的位置应靠近目标断层。又因为钻井数据一般具有二维性,因此现有的方法利用钻井资料分析处理时,常常只能进行简单的二维数据计算分析,不能进行三维分析。相对地,分析处理范围有限,不能很好地满足实际的钻探需要。而本发明实施例中使用的三维地震数据满足了三维分析的需要。
在步骤102中,根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型,在一实施例中,根据所述断面数据、层位数据和测井信息,建立三维地质模型,包括:
根据层位数据,建立层位模型;
根据断面数据,建立断层模型;
基于测井信息,校正层位模型和断层模型;
根据校正后的层位模型和校正后的断面模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型。
具体实施时,为了根据解释得到的层位数据建立各个目的层层面的构造层面模型,具体实行时,在构造层面建模中,建立关键面的构造模型(关键面是指数据信息比较多的,能准确建立构造面的层面,如研究区的油组)。
在一实施例中,根据层位数据,建立层位模型,包括:
S1:根据层位数据,构造层位构造面模型;
将准备好的层位数据加载到专业的建模软件如petrol或者GPTmodel中,按照克里金差值、最小二乘等方法构造层位构造面模型;其中,可以根据层面的接触关系处理地层剥蚀、等地质情况;可以使用地震解释或者其他软件的成果数据做为约束条件,可以先生成地层体模型,同时在关键层面的约束下建立小层的构造层面模型;然后将各个小层根据需要进行垂向网格细分。
S2:在三维空间中分析层位构造面模型中上下相邻层位的关系,剔除引起串层或交叉的层位畸形点,获得新的层位构造面模型。
S3:对所述新的层位构造面模型进行平滑处理,获得层位模型,其中,平滑次数可以根据研究精度需要进行设定。
具体实施时,断层模型构建的关键是对各个断层之间拓扑关系,几何特征,属性特征进行可视化表达,将断层视为地层单元的边界,构建断层模型;利用剖面数据采用断面法构建断层模型,基于轮廓线重构的断层面构建发,油断层面的倾向和倾角构建小的断层面,有许多断层面构建断层模型,对模型的空间长沙县,空间计算,空间分析等,实现一种能解决上述问题的断层数据模型。具体实施时,可按照以下步骤进行:
在一实施例中,根据断面数据,建立断层模型,包括:
S1:剔除断面数据中的畸形数据点,获得断楞数据,所述断楞数据包括断层线数据、断层点构成数据;
将解释好的断层数据加载进专业的建模软件如petrol或者GPTmodel中,在三维视窗下分析断层数据,选定需要建模的范围,剔除畸形数据点,具体的数据格式有Irap、Zmap、CPS-3等。
S2:根据断楞数据,生成断层面模型;
具体可以采用克里金差值等方法,生成断层面模型。
S3:分析断层面模型中断层面的交接关系,对断层面模型中不符合预设要求的部分矫正,获得断层模型,其中,预设要求包括地质规律和特征,可以通过编辑的方式进行矫正。
在一实施例中,基于测井信息,校正层位模型和断层模型,包括:
S1:根据测井信息,获得井上分层体系数据;
可以根据区域地质认识和测井信息,建立井上分层体系,得到井上分层体系数据,其中,分层体系名称和地震上层位名称保持一致。
S2:利用井上分层体系数据,对层位模型进行校正,得到校正后的层位模型,使得层位模型和井上分层体系趋于一致,层位构造形态分布更加符合真实地质情况。
S3:利用校正后的层位模型,对断层模型进行校正,校正时,可以同步编辑,使得层位模型和断层模型更好的匹配,从而提高建模精度和质量。
在步骤103中,对所述三维地质模型中的断面进行网格化处理,获得网格模型,可以根据研究精度需要设置网格大小,网格化断面模型,具体可以为:
S1:根据研究需求和断面-层面交接关系,初步选定需要网格化的断面范围。
S2:根据断面几何形态,选定网格的类型,如三角网格或者是矩形网格。
S3:将所述目标断层的空间断面展布进行网格化,其中网格化步长和网格大小等参数可以根据实际研究需要进行设定,如4m×4m。
特别的,为了提高网格质量,在网格化过程中,网格主方向选择断面走势方向,断面局部区域为了精度要求,可以进行网格加密。
在一实施例中,所述方法还包括:
根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移。
图2为本发明实施例中断面上实际泥岩涂抹的位移的示意图,可以采用如下公式,根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移:
Si=Di/sina
其中,Si为第i个网格点内的断面上实际泥岩涂抹的位移;Di为第i个网格点内的断距,单位为cm;a为第i个网格点内的倾角。
在一实施例中,所述方法还包括:
根据研究区的测井信息获得研究区的地层岩性数据;
根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度。
在一实施例中,测井信息包括泥质含量数据、伽马数据、声波数据和电阻率数据中的至少一种。
在上述实施例中,地层岩性数据包括地层岩性,可以先研究区的测井信息标定研究区的地层岩性,具体包括将las格式测井曲线加载进专业的测井解释软件Geolog或者其他软件中,其中,las格式测井曲线包括泥质含量曲线、伽马曲线、声波曲线和电阻率曲线。
在一个实施例中,可以根据泥质含量曲线、伽马曲线,标定研究区的地层岩性,具体包括:根据研究区砂泥岩响应特征,确定砂泥岩泥质含量区分值和伽马区分值,研究区中大于泥质含量曲线中泥质含量数据大于泥质含量区分值且伽马曲线中伽马数据大于伽马区分值的部分表示泥质含量高,标定为泥岩;否则标定为砂岩,以此得到整个研究区地层的岩性信息。
在其他实施例中,也可以利用其它曲线进行岩性标定获得岩性信息。例如,利用伽马曲线和电阻率曲线进行标定,根据研究区砂泥岩响应特征,将研究区中伽马曲线中伽马数据小于伽马区分值且电阻率曲线中电阻率数据大于电阻率区分的部分标定为砂岩,否则标定为泥岩,从而得到整个研究区地层的岩性信息。
之后,根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度,特别的,当网格点内地层岩性数据加载后,如果泥岩层厚度大于网格大小,说明整个网格点都是泥岩,泥岩层厚度既为纵向上断距大小,如果泥岩层厚度小于纵向上断距大小,既按实际泥岩层厚度参与运算。
而每个网格点内的深度数据可以根据网格化后的目标断层的空间断层展布,从带深度的断面数据中直接读出。
在一实施例中,获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,包括:
采用如下公式,根据网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度:
Li=Hi/Si
其中,Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Hi为第i个网格点内的泥岩层厚度;Si为断面上实际泥岩涂抹的位移;
采用如下公式,根据每个网格点内的深度数据,计算每个网格点内的断面压力:
Pi=10-3Mi(ρl-ρw)·g·cosa
其中,Pi为第i个网格点内的断面压力;Mi为第i个网格点内的深度数据;ρl为上覆地层平均密度;ρw为地层水的平均密度;g为重力加速度;a为第i个网格点内的倾角。
在上述实施例中,由于断层之上静水压力对断层不产生作用,断面所受的促使其愈合的压力主要是上覆地层岩石骨架重力产生的压力,上覆地层平均密度由断裂带附近测井数据得到,地层水的平均密度由测试数据得到。
在一实施例中,采用如下公式,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数:
其中,C为三维断面疏导系数;Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Pi为第i个网格点内的断面压力。
获得了整个断面的三维断面疏导系数,还可以获得每个网格点的疏导系数,从而形成疏导系数图版,便于观察。
在步骤106中,在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导,这里需要用到三维断面疏导系数阈值。
在一个实施例中,三维断面疏导系数阈值根据已钻断块圈闭的油藏解剖确定,即分别计算多个已钻断块圈闭断面的不同深度点的三维断面疏导系数,得到三维断面疏导系数与深度的关系图版,断层两盘联通的断层对应的三维断面疏导系数即为三维断面疏导系数阈值。断层两盘联通的断层的判断方法为:如果断层两盘地层流体性质如地层水、地层原油等的性质趋于一致,既断层两盘地层极可能联通;还可以通过断层两盘地层压力来判断,如果压力想同或者趋于一致,既断层两盘地层联通,等等,以上方法不一而足。这里断块圈闭指的是,在断层控制下的形成的圈闭。相应地,已钻断块圈闭就是该圈闭已经钻井,根据钻井情况就可确定相应断层的联通疏导性能。
在本发明实施例中,根据钻井资料,即已钻断块圈闭资料,可以发现:在已钻断块圈闭中,三维断面疏导系数大于等于5.8/Mpa.m的区域疏导性好,小于5.8/Mpa.m的区域疏导性不好;进而可以将5.8/Mpa.m确定为三维断面疏导系数阈值,用于确定其他断层的网格点的连通疏导性。例如,对于断层的其他网格点,如果该网格点的三维断面疏导系数大于等于5.8/Mpa.m,则确定该网格点疏导性好,达到油气运移要求;如果该网格点的三维断面疏导系数小于5.8/Mpa.m,则确定该网格点疏导性差,未达到疏导油气运移要求。
在一实施例中,所述方法还包括:在研究区三维断层疏导时,确定对所述研究区的断块圈闭进行钻井。
本发明实施例提出的方法还可以预测油气沿着一系列顺向断块或者反向断块的运移路径和轨迹,分析油气成藏情况;还可以判断断层控制的注水层或者油层是否联通,从而制定科学合理的工作制度,从而提高断块圈闭的钻探成功率和开发经济效益。
基于上述实施例,本发明提出如下一个实施例来说明三维断层疏导性的定量确定方法的详细流程,图3为本发明实施例中三维断层疏导性的定量确定方法的详细流程图,如图3所示:
步骤301,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
步骤302,根据层位数据,构造层位构造面模型;
步骤303,在三维空间中分析层位构造面模型中上下相邻层位的关系,剔除引起串层或交叉的层位畸形点,获得新的层位构造面模型;
步骤304,对所述新的层位构造面模型进行平滑处理,获得层位模型;
步骤305,剔除断面数据中的畸形数据点,获得断楞数据,所述断楞数据包括断层线数据、断层点构成数据;
步骤306,根据断楞数据,生成断层面模型;
步骤307,分析断层面模型中断层面的交接关系,对断层面模型中不符合预设要求的部分矫正,获得断层模型;
步骤308,根据测井信息,获得井上分层体系数据;
步骤309,利用井上分层体系数据,对层位模型进行校正,得到校正后的层位模型;
步骤310,利用校正后的层位模型,对断层模型进行校正;
步骤311,根据校正后的层位模型和校正后的断面模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型;
步骤312,对所述三维地质模型中的断面进行网格化处理,获得网格模型;
步骤313,根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移;
步骤314,根据研究区的测井信息获得研究区的地层岩性数据;
步骤315,根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度;
步骤316,获取网格模型的断面中每个网格点内的深度数据;
步骤317,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
步骤318,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
步骤319,在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。
当然,可以理解的是,上述详细流程还可以有其他变化例,相关变化例均应落入本发明的保护范围。
下面给出一具体实施例,来说明三维断层疏导性的定量确定方法的具体应用。
S11:将目标区叠后地震数据加载进解释软件Landmark中,获得断层数据和层位数据。
S12:将断层数据和层位数据加载进专业的建模软件petrel中,分别构建层位模型和断层模型,得到三维地质模型,基于测井信息,校正层位模型和断层模型,根据校正后的层位模型和校正后的断面模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型,图4为本发明实施例中三维地质模型的示意图。
S13:对所述三维地质模型中的断面进行网格化处理,获得网格模型,具体用于网格化的网格的大小根据研究需要自定义,此次实施例中采用25×25m的网格步长大小,然后获取每个网格点的倾角和断距。具体可以参阅图2,图2中的P、Q为断面上任意两点,根据图2可知上述两点的断面倾角显然不同。
S14:获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,其中,每个网格点内的泥岩层厚度根据研究区的地层岩性数据确定,图5为本发明实施例中岩性曲线图版的示意图,通过将las格式测井曲线加载进专业的测井解释软件Geolog获得。
S15:根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数。
S16,在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导,其中,三维断面疏导系数阈值根据已钻两个断块圈闭确定,根据钻井资料,两个已钻断块圈闭关键成藏期为白垩系,均为干井,且均显示了良好的油气,表明该区油源充足。图6为本发明实施例中三维断面疏导系数与深度的关系图版,其中长20、长19和长2为三口钻井,从图版分析可知,发现高部位含油的控圈断层其三维断面疏导系数在5.85Mp.s以上;而失利目的层高部位含油的控圈断层其三维断面疏导系数在5.85Mp.s以下。即可以确定高部位含油的控圈断层其三维断面疏导系数阈值在5.85Mp.s。由于该三维断面疏导系数阈值来源于已钻油藏,已经得到了钻井证实,因此,可以代表该目标区的地质特征。
S17,将目标区三维断面疏导系数C大于15Mp.s的分为疏导性优的I类断层,将三维断面疏导系数C大于10Mp.s而小于等于15Mp.s的分为疏导性好的II类断层,将三维断面疏导系数C小于10Mp.s大于5.8Mp.s的断层分为疏导性一般的III类断层,将三维断面疏导系数C小于5.8Mpa.s的断层疏导性分为不疏导的封闭的IV类断层,其中,I类断层控制疏导形成的高部位圈闭,油气易跨越断层侧向运移向上运移形成聚集成藏,为首选部署圈闭;II类断层控制的圈闭次之;III类断层控制的圈闭为风险性较高的圈闭;IV类断层控制的圈闭为不疏导的断层,油气不易跨越断层形成聚集,其高部位不会形成圈闭,不建议钻探。详见表1的疏导性定量评价标准,按照上述思路部署井位从而避免了风险,为勘探部署和钻井提供可靠的指导。
表1
综上所述,在本发明实施例提出的方法中,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;对所述三维地质模型中的断面进行网格化处理,获得网格模型;获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。在上述过程中,综合考虑了研究区叠后三维地震数据、测井信息,最后建立了三维地质模型,并在获得网格模型后,定量获得了每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,进而定量计算了三维断面疏导系数,使得根据三维断面疏导系数判断研究区三维断层是否疏导的过程更加准确。
本发明实施例还提出一种三维断层疏导性的定量确定装置其原理与三维断层疏导性的定量确定方法类似,这里不再赘述。
图7为发明实施例中三维断层疏导性的定量确定装置的示意图,如图7所述,该装置包括:
数据获得模块701,用于根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
三维地质模型建立模块702,用于根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
网格模型获得模块703,用于对所述三维地质模型中的断面进行网格化处理,获得网格模型;
第一计算模块704,用于获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
第二计算模块705,用于根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
疏导性确定模块706,用于在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。
在一实施例中,测井信息包括泥质含量数据、伽马数据、声波数据和电阻率数据中的至少一种。
在一实施例中,三维地质模型建立模块702包括:
层位模型建立模块7021,用于根据层位数据,建立层位模型;
断面模型建立模块7022,用于根据断面数据,建立断层模型;
校正模块7023,用于基于测井信息,校正层位模型和断层模型;
生成模块7024,用于根据校正后的层位模型和校正后的断面模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型。
在一实施例中,层位模型建立模块7021具体用于:
根据层位数据,构造层位构造面模型;
在三维空间中分析层位构造面模型中上下相邻层位的关系,剔除引起串层或交叉的层位畸形点,获得新的层位构造面模型;
对所述新的层位构造面模型进行平滑处理,获得层位模型。
在一实施例中,断面模型建立模块7022具体用于:
剔除断面数据中的畸形数据点,获得断楞数据,所述断楞数据包括断层线数据、断层点构成数据;
根据断楞数据,生成断层面模型;
分析断层面模型中断层面的交接关系,对断层面模型中不符合预设要求的部分矫正,获得断层模型。
在一实施例中,校正模块7023具体用于:
根据测井信息,获得井上分层体系数据;
利用井上分层体系数据,对层位模型进行校正,得到校正后的层位模型;
利用校正后的层位模型,对断层模型进行校正。
在一实施例中,所述装置还包括第三计算模块707,用于:
根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移。
在一实施例中,所述装置还包括泥岩层厚度获得模块708,用于:
根据研究区的测井信息获得研究区的地层岩性数据;
根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度。
在一实施例中,第一计算模块704具体用于:
采用如下公式,根据网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度:
Li=Hi/Si
其中,Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Hi为第i个网格点内的泥岩层厚度;Si为断面上实际泥岩涂抹的位移;
采用如下公式,根据每个网格点内的深度数据,计算每个网格点内的断面压力:
Pi=10-3Mi(ρl-ρw)·g·cosa
其中,Pi为第i个网格点内的断面压力;Mi为第i个网格点内的深度数据;ρl为上覆地层平均密度;ρw为地层水的平均密度;g为重力加速度;a为第i个网格点内的倾角。
在一实施例中,第二计算模块705具体用于:
采用如下公式,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数:
其中,C为三维断面疏导系数;Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Pi为第i个网格点内的断面压力。
综上所述,在本发明实施例提出的装置中,根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;对所述三维地质模型中的断面进行网格化处理,获得网格模型;获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导。在上述过程中,综合考虑了研究区叠后三维地震数据、测井信息,最后建立了三维地质模型,并在获得网格模型后,定量获得了每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,进而定量计算了三维断面疏导系数,使得根据三维断面疏导系数判断研究区三维断层是否疏导的过程更加准确。
本申请的实施例还提供一种计算机设备,图8为本发明实施例中计算机设备的示意图,该计算机设备能够实现上述实施例中的三维断层疏导性的定量确定方法中全部步骤,所述计算机设备具体包括如下内容:
处理器(processor)801、存储器(memory)802、通信接口(CommunicationsInterface)803和通信总线804;
其中,所述处理器801、存储器802、通信接口803通过所述通信总线804完成相互间的通信;所述通信接口803用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输;
所述处理器801用于调用所述存储器802中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的三维断层疏导性的定量确定方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质,能够实现上述实施例中的三维断层疏导性的定量确定方法中全部步骤,所述计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述实施例中的三维断层疏导性的定量确定方法的全部步骤。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (18)
1.一种三维断层疏导性的定量确定方法,其特征在于,包括:
根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
对所述三维地质模型中的断面进行网格化处理,获得网格模型;
获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导;
获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,包括:
采用如下公式,根据网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度:
Li=Hi/Si
其中,Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Hi为第i个网格点内的泥岩层厚度;Si为断面上实际泥岩涂抹的位移;
采用如下公式,根据每个网格点内的深度数据,计算每个网格点内的断面压力:
Pi=10-3Mi(ρl-ρw)·g·cosa
其中,Pi为第i个网格点内的断面压力;Mi为第i个网格点内的深度数据;ρl为上覆地层平均密度;ρw为地层水的平均密度;g为重力加速度;a为第i个网格点内的倾角;
采用如下公式,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数:
其中,C为三维断面疏导系数;Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Pi为第i个网格点内的断面压力。
2.如权利要求1所述的三维断层疏导性的定量确定方法,其特征在于,测井信息包括泥质含量数据、伽马数据、声波数据和电阻率数据中的至少一种。
3.如权利要求1所述的三维断层疏导性的定量确定方法,其特征在于,根据所述断面数据、层位数据和测井信息,建立三维地质模型,包括:
根据层位数据,建立层位模型;
根据断面数据,建立断层模型;
基于测井信息,校正层位模型和断层模型;
根据校正后的层位模型和校正后的断层模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型。
4.如权利要求3所述的三维断层疏导性的定量确定方法,其特征在于,根据层位数据,建立层位模型,包括:
根据层位数据,构造层位构造面模型;
在三维空间中分析层位构造面模型中上下相邻层位的关系,剔除引起串层或交叉的层位畸形点,获得新的层位构造面模型;
对所述新的层位构造面模型进行平滑处理,获得层位模型。
5.如权利要求3所述的三维断层疏导性的定量确定方法,其特征在于,根据断面数据,建立断层模型,包括:
剔除断面数据中的畸形数据点,获得断楞数据,所述断楞数据包括断层线数据、断层点构成数据;
根据断楞数据,生成断层面模型;
分析断层面模型中断层面的交接关系,对断层面模型中不符合预设要求的部分矫正,获得断层模型。
6.如权利要求3所述的三维断层疏导性的定量确定方法,其特征在于,基于测井信息,校正层位模型和断层模型,包括:
根据测井信息,获得井上分层体系数据;
利用井上分层体系数据,对层位模型进行校正,得到校正后的层位模型;
利用校正后的层位模型,对断层模型进行校正。
7.如权利要求1所述的三维断层疏导性的定量确定方法,其特征在于,还包括:
根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移。
8.如权利要求1所述的三维断层疏导性的定量确定方法,其特征在于,还包括:
根据研究区的测井信息获得研究区的地层岩性数据;
根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度。
9.一种三维断层疏导性的定量确定装置,其特征在于,包括:
数据获得模块,用于根据研究区叠后三维地震数据,获得研究区的断面数据和层位数据;
三维地质模型建立模块,用于根据所述断面数据、层位数据和研究区的测井信息,建立三维地质模型;
网格模型获得模块,用于对所述三维地质模型中的断面进行网格化处理,获得网格模型;
第一计算模块,用于获取网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移和深度数据,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力;
第二计算模块,用于根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数;
疏导性确定模块,用于在所述三维断面疏导系数大于三维断面疏导系数阈值时,确定研究区三维断层疏导,否则确定研究区三维断层不疏导;
第一计算模块具体用于:
采用如下公式,根据网格模型的断面中每个网格点内的泥岩层厚度、断面上实际泥岩涂抹的位移,计算每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度:
Li=Hi/Si
其中,Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Hi为第i个网格点内的泥岩层厚度;Si为断面上实际泥岩涂抹的位移;
采用如下公式,根据每个网格点内的深度数据,计算每个网格点内的断面压力:
Pi=10-3Mi(ρl-ρw)·g·cosa
其中,Pi为第i个网格点内的断面压力;Mi为第i个网格点内的深度数据;ρl为上覆地层平均密度;ρw为地层水的平均密度;g为重力加速度;a为第i个网格点内的倾角;
第二计算模块具体用于:
采用如下公式,根据每个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度和断面压力,计算三维断面疏导系数:
其中,C为三维断面疏导系数;Li为第i个网格点内的泥岩层沿断面实际涂抹位移上的泥岩分布厚度;Pi为第i个网格点内的断面压力。
10.如权利要求9所述的三维断层疏导性的定量确定装置,其特征在于,测井信息包括泥质含量数据、伽马数据、声波数据和电阻率数据中的至少一种。
11.如权利要求9所述的三维断层疏导性的定量确定装置,其特征在于,三维地质模型建立模块包括:
层位模型建立模块,用于根据层位数据,建立层位模型;
断面模型建立模块,用于根据断面数据,建立断层模型;
校正模块,用于基于测井信息,校正层位模型和断层模型;
生成模块,用于根据校正后的层位模型和校正后的断层模型,处理层位和断面交接关系,形成包含层位和断面的三维地质模型。
12.如权利要求11所述的三维断层疏导性的定量确定装置,其特征在于,层位模型建立模块具体用于:
根据层位数据,构造层位构造面模型;
在三维空间中分析层位构造面模型中上下相邻层位的关系,剔除引起串层或交叉的层位畸形点,获得新的层位构造面模型;
对所述新的层位构造面模型进行平滑处理,获得层位模型。
13.如权利要求11所述的三维断层疏导性的定量确定装置,其特征在于,断面模型建立模块具体用于:
剔除断面数据中的畸形数据点,获得断楞数据,所述断楞数据包括断层线数据、断层点构成数据;
根据断楞数据,生成断层面模型;
分析断层面模型中断层面的交接关系,对断层面模型中不符合预设要求的部分矫正,获得断层模型。
14.如权利要求11所述的三维断层疏导性的定量确定装置,其特征在于,校正模块具体用于:
根据测井信息,获得井上分层体系数据;
利用井上分层体系数据,对层位模型进行校正,得到校正后的层位模型;
利用校正后的层位模型,对断层模型进行校正。
15.如权利要求9所述的三维断层疏导性的定量确定装置,其特征在于,还包括第三计算模块,用于:
根据每个网格点内的倾角和断距,计算每个网格点内的断面上实际泥岩涂抹的位移。
16.如权利要求9所述的三维断层疏导性的定量确定装置,其特征在于,还包括泥岩层厚度获得模块,用于:
根据研究区的测井信息获得研究区的地层岩性数据;
根据地层岩性数据,获得网格模型的断面中每个网格点内的泥岩层厚度。
17.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8任一项所述方法。
18.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至8任一项所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010903256.9A CN114114411B (zh) | 2020-09-01 | 2020-09-01 | 三维断层疏导性的定量确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010903256.9A CN114114411B (zh) | 2020-09-01 | 2020-09-01 | 三维断层疏导性的定量确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114114411A CN114114411A (zh) | 2022-03-01 |
CN114114411B true CN114114411B (zh) | 2023-08-22 |
Family
ID=80360685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010903256.9A Active CN114114411B (zh) | 2020-09-01 | 2020-09-01 | 三维断层疏导性的定量确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114114411B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114994282A (zh) * | 2022-06-01 | 2022-09-02 | 国家石油天然气管网集团有限公司 | 一种面向交变应力环境的断层滑移判别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194253A (zh) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | 一种面向三维地质层面结构的四面体网格生成方法 |
CN103267835A (zh) * | 2013-04-19 | 2013-08-28 | 山东大学 | 可预制断层的大型流-固耦合模型试验台架及试验方法 |
CN105089615A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种基于油藏模型的测井数据历史回归处理方法 |
CN108508484A (zh) * | 2018-02-23 | 2018-09-07 | 中国石油大学(华东) | 一种适用于张扭性断层泥岩涂抹效果的评价方法 |
CN110321595A (zh) * | 2019-06-06 | 2019-10-11 | 中国石油大学(华东) | 一种测井提取静态品质系数的断层封闭性评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6751558B2 (en) * | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
-
2020
- 2020-09-01 CN CN202010903256.9A patent/CN114114411B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194253A (zh) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | 一种面向三维地质层面结构的四面体网格生成方法 |
CN103267835A (zh) * | 2013-04-19 | 2013-08-28 | 山东大学 | 可预制断层的大型流-固耦合模型试验台架及试验方法 |
CN105089615A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种基于油藏模型的测井数据历史回归处理方法 |
CN108508484A (zh) * | 2018-02-23 | 2018-09-07 | 中国石油大学(华东) | 一种适用于张扭性断层泥岩涂抹效果的评价方法 |
CN110321595A (zh) * | 2019-06-06 | 2019-10-11 | 中国石油大学(华东) | 一种测井提取静态品质系数的断层封闭性评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114114411A (zh) | 2022-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11209561B2 (en) | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function | |
RU2661489C1 (ru) | Способ комплексирования исходных данных для уточнения фильтрационного строения неоднородных карбонатных коллекторов | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
CA2572449A1 (en) | Hydrodynamics-based gridding geologic modeling (hydro-gridding) | |
Aarre et al. | Seismic detection of subtle faults and fractures | |
CA2692425A1 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
CN110632658B (zh) | 断层的侧向封闭性分析方法及装置 | |
US20140278318A1 (en) | Natural Resource Reservoir Modeling | |
Trippetta et al. | Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy | |
Abeed et al. | Modeling of petrophysical properties and reserve estimation of Mishrif formation-Garraf oil field | |
Kurah et al. | Reservoir characterization and volumetric estimation of reservoir fluids using simulation and analytical methods: a case study of the coastal swamp depobelt, Niger Delta Basin, Nigeria | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN114114411B (zh) | 三维断层疏导性的定量确定方法及装置 | |
Suarez et al. | Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale | |
Alpak et al. | Shale-drape modeling for the geologically consistent simulation of clastic reservoirs | |
CN115857047B (zh) | 一种地震储层综合预测方法 | |
Emujakporue | Petrophysical properties distribution modelling of an onshore field, Niger Delta, Nigeria | |
CN113376692B (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
Qi et al. | Correlation of seismic attributes and geomechanical properties to the rate of penetration in the Mississippian Limestone, Oklahoma | |
Taura et al. | Structural and petrophysical controls on the remaining fluid distribution in the reservoirs of Gharif Formation before abandonment | |
Haris et al. | Reservoir Compartment Assessment: A Case Study Of Bangko And Bekasap Formation, Central Sumatra Basin Indonesia | |
CN113341480B (zh) | 分频rgb切片与多属性融合的砂质水合物储层预测方法 | |
Podnebesnykh et al. | New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field | |
Tamulonis | A methodology for integrating unconventional geologic and engineering data into a geocellular model | |
Kondratyev et al. | Submarine fan reservoir architecture and heterogeneity influence on hard-to-recover reserves. Achimov Fm |
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 |