CN111815769B - 逆冲推覆构造带构造的建模方法、计算设备及存储介质 - Google Patents
逆冲推覆构造带构造的建模方法、计算设备及存储介质 Download PDFInfo
- Publication number
- CN111815769B CN111815769B CN201910292415.3A CN201910292415A CN111815769B CN 111815769 B CN111815769 B CN 111815769B CN 201910292415 A CN201910292415 A CN 201910292415A CN 111815769 B CN111815769 B CN 111815769B
- Authority
- CN
- China
- Prior art keywords
- target
- thrust
- fault
- data
- model
- 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 104
- 238000003860 storage Methods 0.000 title claims abstract description 11
- 238000010276 construction Methods 0.000 title claims description 80
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 49
- 238000005253 cladding Methods 0.000 claims abstract description 45
- 238000004458 analytical method Methods 0.000 claims abstract description 38
- 238000009826 distribution Methods 0.000 claims abstract description 21
- 208000010392 Bone Fractures Diseases 0.000 claims description 29
- 206010017076 Fracture Diseases 0.000 claims description 29
- 238000005516 engineering process Methods 0.000 claims description 27
- 230000000877 morphologic effect Effects 0.000 claims description 23
- 238000011161 development Methods 0.000 claims description 21
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 238000005755 formation reaction Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000011160 research Methods 0.000 description 9
- 238000005553 drilling Methods 0.000 description 6
- 238000004445 quantitative analysis Methods 0.000 description 5
- 238000005520 cutting process Methods 0.000 description 4
- 238000007405 data analysis Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000011248 coating agent Substances 0.000 description 3
- 238000000576 coating method Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000011010 flushing procedure Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 210000001331 nose Anatomy 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 230000005389 magnetism Effects 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
- 238000002360 preparation method Methods 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种逆冲推覆构造带构造的建模方法、装置、计算设备及存储介质,属于计算机技术领域,所述方法包括:根据反映逆冲推覆构造带构造的各地层的形态以及地理位置分布的形态数据,和,二维声学波动方程建立目标构造模型,对该目标构造模型展开正演,得到目标量版,根据目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型,根据目标多层面差值算法对第一层面网格模型进行优化,得到第二层面网格模型,可见本方案能够解决角点网格方法建立的三维模型,在逆冲推覆构造带的断层部分发生变形,不能真实反映逆冲推覆构造带的三维结构的问题,提高了用于反映逆冲推覆构造带的三维模型的准确度。
Description
技术领域
本申请实施例涉及石油工业油气田开发领域,特别涉及一种逆冲推覆构造带构造的建模方法、计算设备及存储介质。
背景技术
在石油工业的地质勘探中,技术人员通常采用建模的方法,对地质构造解释。同时,建模也是对地质结构进行三维再现的重要手段。
在油气藏勘探和开采领域中,逆冲推覆构造带是存储油气的有利构造。因此,对逆冲推覆构造带进行三维再现,是开采逆冲推覆构造带中的油气藏的必要准备。在一些技术中,技术人员利用角点网格方法对逆冲推覆构造带进行建模。
然而,角点网格方法建立的三维模型在逆冲推覆构造带的断层部分发生变形,不能真实反映逆冲推覆构造带的三维结构。
发明内容
本申请实施例提供了一种逆冲推覆构造带构造的建模方法、装置、计算设备及存储介质,可以解决角点网格方法建立的三维模型在逆冲推覆构造带的断层部分发生变形,不能真实反映逆冲推覆构造带的三维结构的问题。所述技术方案如下:
根据本申请的一方面内容,提供了一种逆冲推覆构造带构造的建模方法,所述方法包括:
根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;
对所述目标构造模型开展正演,获取目标量版,所述目标量版包括指示至少两种类型且放大n倍的断块构造样式,n大于1;
根据所述目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型;
根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
根据本申请的另一方面内容,提供了一种逆冲推覆构造带构造的建模装置,所述装置包括:
第一建模模块,用于根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;
量版获取模块,用于对所述目标构造模型开展正演,获取目标量版,所述目标量版包括至少两种类型且放大n倍的断块构造样式,n大于1;
第二建模模块,用于根据所述目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型;
第三建模模块,用于根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
根据本申请的另一方面内容,提供了一种计算设备,所述计算设备包括处理器和存储器,所述存储器中存储有至少一条指令,所述指令由所述处理器加载并执行以实现如本申请实施提供的逆冲推覆构造带构造的建模方法。
根据本申请的另一方面内容,提供了一种计算机可读存储介质,所述存储介质中存储有至少一条指令,所述指令由处理器加载并执行以实现如本申请实施提供的逆冲推覆构造带构造的建模方法。
本申请实施例提供的技术方案带来的有益效果可以包括:
根据形态数据和二维声学波动方程建立目标构造模型,其中,形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布,对目标构造模型开展正演,获取目标量版,该目标量版包括指示至少两种类型且放大n倍的断块构造样式,根据目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型,根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二网格模型,该第二层面网格模型的精度高于第一层面网格模型,且第二层面网格模型是逆冲推覆构造带构造的三维模型。由于通过本申请能够利用目标构造模型的正演获得目标量版,并在此基础上对不同序级的断层进行展布分析获得第一层面网格模型,再对该第一层面网格模型进行优化,得到第二层面网格模型,可见,本实施例能够获得较之角点网格方法更为精确的模型,从而提高了用于反映逆冲推覆构造带构造的模型的精准度,有利于后续的地质结构分析和油气开采工作。
附图说明
为了更清楚地介绍本申请实施例中的技术方案,下面将对本申请实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本申请一个示例性实施例提供的计算设备的结构框图;
图2是本申请一个示例性实施例提供的逆冲推覆构造带构造的建模方法的流程图;
图3是本申请另一个示例性实施例提供的逆冲推覆构造带构造的建模方法流程图;
图4是本申请实施例提供的一种复杂推覆构造三维构造建模技术的实现示意图;
图5是本申请提供的一种C气藏所在的构造区带的沉积微相模型剖面图;
图6是本申请提供的一种2D声学波动方程正演模拟Y型断层模板;
图7是本申请提供的一种逆冲推覆构造地区三维区内构造类型特征图;
图8是本申请提供的一种复杂构造带inline1528线不同断层级别、组合关系时深转换地震剖面图;
图9是本申请实施例提供的逆冲复杂构造带非结构化网格建模技术反映复杂构造的示意图;
图10是本申请实施例提供的一种不同低序级断层平面网格分布图;
图11是本申请实施例提供的一种旋回差异算法确定单井最小垂厚对比图;
图12是本申请实施例提供的一种逆冲推覆构造层面中的算法模拟适应性分析图;
图13是本申请实施例提供的一种复杂逆冲推覆构造带三维构造模型图;
图14是本申请一个示例性实施例提供的逆冲推覆构造带构造的建模装置的结构框图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请实施方式作进一步地详细描述。
下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在本申请的描述中,需要理解的是,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性。在本申请的描述中,需要说明的是,除非另有明确的规定和限定,术语“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本申请中的具体含义。此外,在本申请的描述中,除非另有说明,“多个”是指两个或两个以上。“和/或”,描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
在本申请实施例中,构造建模是石油地质行业常用的构造解释、三维地质建模必备的过程和手段。在构造建模中,主要的技术手段包括地质测控手段、地震测量手段和钻井采样手段等。技术人员能够通过上述多种技术手段采集到原始的地址数据,再利用层位追踪、蚂蚁追踪、三维地址建模软件等行为构造描述和断层预测,最终为钻井实施及油气藏认识提供三维构造模型。
通常目前构造建模技术是针对构造平缓、断裂发育特征清楚的气藏应用。针对复杂逆冲推覆的构造建模技术尚处于技术攻关阶段,业内没有成熟的方案可以参考借鉴。
鉴于本申请提供的一种逆冲推覆构造带构造的建模方法,计算设备能够针对逆冲推覆构造碳酸盐岩礁滩气藏开展精细的三维构造建模工作,建立比较接近地下真实情况的高精度三维构造模型,可以对气藏强非均质性特征进行精细表征,同时有利于开展后期的气藏数值模拟研究,为气藏的储量评价和后续开发提供技术支撑。
另外,本申请实施例提供的逆冲推覆构造带构造的建模方法,能够用于碳酸盐气藏的三维建模过程中,能有效解决复杂逆冲推覆构造气藏断裂复杂、断块多等在三维构造建模中难于表征的问题。
示例性地,本申请实施例所示的逆冲推覆构造带构造的建模方法,可以应用在计算设备中,该计算设备具备显示屏且具备逆冲推覆构造带构造的建模的功能。计算设备可以包括台式电脑、电脑一体机、服务器、工作站和超级计算机等设备。
请参考图1,图1是本申请一个示例性实施例提供的计算设备的结构框图,如图1所示,该计算设备包括处理器120和存储器140,所述存储器140中存储有至少一条指令,所述指令由所述处理器120加载并执行以实现如本申请各个方法实施例所述的逆冲推覆构造带构造的建模方法。
在本申请中,计算设备100是具备海量数据处理功能的电子设备。当计算设备100运行数据分析软件时,计算设备100能够根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;对所述目标构造模型开展正演,获取目标量版,所述目标量版包括指示至少两种类型且放大n倍的断块构造样式,n大于1;根据所述目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型;根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
处理器120可以包括一个或者多个处理核心。处理器120利用各种接口和线路连接整个计算设备100内的各个部分,通过运行或执行存储在存储器140内的指令、程序、代码集或指令集,以及调用存储在存储器140内的数据,执行计算设备100的各种功能和处理数据。可选的,处理器120可以采用数字信号处理(Digital Signal Processing,DSP)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)、可编程逻辑阵列(ProgrammableLogic Array,PLA)中的至少一种硬件形式来实现。处理器120可集成中央处理器(CentralProcessing Unit,CPU)、图像处理器(Graphics Processing Unit,GPU)和调制解调器等中的一种或几种的组合。其中,CPU主要处理操作***、用户界面和应用程序等;GPU用于负责显示屏所需要显示的内容的渲染和绘制;调制解调器用于处理无线通信。可以理解的是,上述调制解调器也可以不集成到处理器120中,单独通过一块芯片进行实现。
存储器140可以包括随机存储器(Random Access Memory,RAM),也可以包括只读存储器(Read-Only Memory)。可选的,该存储器140包括非瞬时性计算机可读介质(non-transitory computer-readable storage medium)。存储器140可用于存储指令、程序、代码、代码集或指令集。存储器140可包括存储程序区和存储数据区,其中,存储程序区可存储用于实现操作***的指令、用于至少一个功能的指令(比如触控功能、声音播放功能、图像播放功能等)、用于实现下述各个方法实施例的指令等;存储数据区可存储下面各个方法实施例中涉及到的数据等。
请参考图2,图2是本申请一个示例性实施例提供的逆冲推覆构造带构造的建模方法的流程图。该逆冲推覆构造带构造的建模方法可以应用在上述所示的计算设备中。在图2中,逆冲推覆构造带构造的建模方法包括:
步骤210,根据形态数据和二维声学波动方程建立目标构造模型,形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布。
在本申请实施例中,计算设备可以是能够进行海量数据处理的电子设备。例如,超级计算机和大型运算设备能够对地质勘探领域中采集到的海量数据进行处理。
在本申请实施例中,计算设备能够获取形态数据。形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布。在一种可能的实现方式中,形态数据包括地震数据、钻井信息数据和测井数据中至少一种。需要说明的是,钻井信息数据可以实现为岩心数据。在一种实际应用的场景中,岩心数据、地震数据和测井数据是两两互为独立来源的数据。本申请中的计算设备能够将上述形态数据进行汇总,组成基础研究数据库。
可选地,在一种可能的实现方式中,形态数据还包括历史地质数据。
可选地,在一种可能的实现方式中,形态数据还包括生产动态数据。
在本申请实施例中,计算设备将根据形态数据,并结合二维声学波动方程建立目标构造模型。
步骤220,对目标构造模型开展正演,获取目标量版,目标量版包括指示至少两种类型且放大n倍的断块构造样式,n大于1。
在本申请实施例中,计算设备将对目标构造模型开展正演,获取目标量版,该目标量版包括至少两种类型且放大n倍的断块构造样式,n大于1。
需要说明的是,正演又称地球物理正演(英文:geophyrsical forwardcalculation),是指在地球物理资料解释理论中,由地质体的赋存状态和物性参数计算该地质体引起的场异常或效应的过程。其中,地质体的赋存状态包括形状、产状和空间位置中至少一种,物性参数包括密度、磁性、电性、弹性和速度中至少一种。
可选地,已知地质体的赋存状态和物性可统称为模型。
可选地,目标量版是经过正演获取的理想化的模板,该目标量版用于指导实际的断层解释。
步骤230,根据目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型。
在本申请实施例中,计算设备能够对至少两个不同序级的断层进行展布分析,建立第一层面网格模型。
可选地,计算设备能够在确定断层影响界限及断层规模的基础上,利用非结构网格法(英文:delaunay)建模技术和小波分解技术,开展不同序级断层展布分析,确定平面网格大小及方位,建立高精度层面网格模型,也即,建立第一层面网格模型。
需要说明的是,平面网格大小和方位是建立第一层面网格模型所需要的参数。非结构网格法建模技术和小波分解技术是用于建模的算法。断层影响界限及断层规模是地质模型数据。当计算设备同时获取上述参数、算法和地质模型数据时,计算设备能够建立第一层面网格模型。
步骤240,根据目标多层面差值算法对第一层面网格模型进行优化,建立第二层面网格模型,第二层面网格模型的精度高于第一层面网格模型,且第二层面网格模型是逆冲推覆构造带构造的三维模型。
在本申请实施例中,计算设备能够根据目标多层面差值算法对第一层面网格模型进行优化,建立第二层面网格模型。需要说明的是,目标多层面差值算法用于提高层面网格模型的精度。其中,第二层面网格模型的精度高于第一层面网格模型,且第二层面网格模型是逆冲推覆构造带构造的三维模型。
综上所述,本实施例提供的逆冲推覆构造带构造的建模方法,根据形态数据和二维声学波动方程建立目标构造模型,其中,形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布,对目标构造模型开展正演,获取目标量版,该目标量版包括指示至少两种类型且放大n倍的断块构造样式,根据目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型,根据目标多层面差值算法对第一层面网格模型进行优化,建立第二网格模型,该第二层面网格模型的精度高于第一层面网格模型,且第二层面网格模型是逆冲推覆构造带构造的三维模型。由于通过本申请能够利用目标构造模型的正演获得目标量版,并在此基础上对不同序级的断层进行展布分析获得第一层面网格模型,再对该第一层面网格模型进行优化,得到第二层面网格模型,可见,本实施例能够获得较之角点网格方法更为精确的模型,从而提高了用于反映逆冲推覆构造带构造的模型的精准度,有利于后续的地质结构分析和油气开采工作。
基于上一个实施例所公开的方案,计算设备还能够在建立模型之前获取目标勘探数据,通过分析目标勘探数据获取逆冲推覆构造带构造的形态数据,并在建立第二层面网格模型时,使用中间层面网格模型作为生成模型的过渡,提高了逆冲推覆构造带构造三维模型的精确度。详细地,请参考如下实施例。
请参见图3,图3是本申请另一个示例性实施例提供的逆冲推覆构造带构造的建模方法流程图。该逆冲推覆构造带构造的建模方法可以应用在上述所示的计算设备中。在图3中,该逆冲推覆构造带构造的建模方法包括:
步骤310,获取目标勘探数据,目标勘探数据是指示逆冲推覆构造带构造在被勘探时收集的地质数据。
在本申请实施例中,计算设备能够获取目标勘探数据,该目标勘探数据是工作人员在针对逆冲推覆构造带构造,进行勘探时收集的地质数据。计算设备可以通过人工输入或者磁盘拷贝等方法获得该目标勘探数据。
在一种可能的实现方式中,计算设备还能够通过获取目标单井的单井数据,目标单井是逆冲推覆构造带中的单井,以获取目标勘探数据。
需要说明的是,单井数据可以包括测井数据、井点分层数据和微相数据中至少一种。可选地,在获取目标单井的单井数据时,本实施例还可以对单井数据进行校正。校正操作包括取奇异值和岩心校正等工作。
步骤320,分析目标勘探数据,获取逆冲推覆构造带构造的形态数据。
在本申请实施例中,计算设备能够通过三维地震和井精细标定技术分析目标勘探数据,获取逆冲推覆构造带构造的形态数据。
在一种可能的实现方式中,计算设备还能够通过小波分解技术处理单井数据,获取逆冲推覆构造带的形态数据,其中,形态数据包括气藏地界数据,以及,目标单井的断层参数。
步骤330,根据形态数据和二维声学波动方程建立目标构造模型。
步骤330的执行过程和步骤210的执行过程相同,此处不再赘述。
步骤340,对目标构造模型开展正演,获取目标量版。
步骤340的执行过程和步骤220的执行过程相同,此处不再赘述。
步骤350,根据目标量版,确定逆冲推覆构造带构造的断层影响界限和断层规模。
在本申请实施例中,计算设备能够根据目标量版,确定逆冲推覆构造带构造的断层影响界限和断层规模。需要说明的是,该操作可以是地质模型参数中的定量分析。
步骤360,根据断层影响界限和断层规模,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型。
在本申请实施例中,计算设备能够在断层影响界限和断层规模的基础上,利用非结构网络法(英文:delaunay)建模技术和小波分解技术,开展不同序级断层展布分析。
步骤370,根据目标多层面差值算法对第一层面网格模型进行优化,得到中间层面网格模型。
在本申请实施例中,计算设备还能够从候选多层面差值算法中确定出目标多层面差值算法。其中,候选多层面差值算法包括克里金算法、最小曲率法和收敛法中至少两种算法。
步骤380,对中间层面网格模型进行井分层数据内插,得到第二层面网格模型。
在本申请实施例中,计算设备将对中间层面网格模型进行井分层数据内插,优化地质断层截切关系,建立层面与断层有效接触关系,从而构建更加复杂精确的断层模型。
综上所述,本申请提供的逆冲推覆构造带构造的建模方法,能够利用目标构造模型的正演获得目标量版,并在此基础上对不同序级的断层进行展布分析获得第一层面网格模型,再对该第一层面网格模型进行优化,得到第二层面网格模型,可见,本实施例能够获得较之角点网格方法更为精确的模型,从而提高了用于反映逆冲推覆构造带构造的模型的精准度,有利于后续的地质结构分析和油气开采工作。
作为本申请另一种能够实现的方式,本申请提供如下执行流程,该执行流程可以应用在如图1所示的计算设备中,该流程包括如下步骤:
本申请可以采用非结构网格算法建立一体化模型,形成“模型正反演、组合定量分析、非结构化网格建立”为核心的复杂推覆构造三维构造建模方法。
步骤(1),收集气藏地震数据、钻井信息数据、测井数据、生产动态数据及地质背景相关资料,运用三维地震与井精细标定技术,对复杂逆冲构造带构造区带进行规律分析及识别划分。
步骤(2),利用步骤(1)中的构造带划分结果、断层产状、构造应力分析参数等数据,采用2D声学波动方程建立构造模型开展正演,最终形成不同类型断块构造样式放大倍数的量版。
步骤(3),在步骤(2)建成的不同类型断块构造样式放大倍数的量版基础上,开展断层级别定量划分,确定断层影响界限及断层规模。
步骤(4)在步骤(3)确定断层影响界限及断层规模的基础上,利用非结构网格法(delaunay)建模技术、小波分解技术,开展不同序级断层展布分析,确定平面网格大小及方位,建立高精度层面网格模型。
步骤(5)在步骤(4)建立的高精度层面网格模型的基础上,开展克里金算法、最小曲率法、收敛法多层面差值算法优选,最终优选收敛法建立层面网格模型。
步骤(6)在步骤(5)在层面网格模型的基础上,开展井分层数据内插,优化断层截切关系,建立层面与断层有效接触关系,构建更加精确的复杂断层模型。
本申请实施例通过上述步骤(1)至步骤(6)组成的三维构造建模方法,能够通过“模型正反演、组合定量分析、非结构化网格建立”为核心的复杂推覆构造三维构造建模技术,突出构造解释与地质建模工作界限,以复杂构造体为研究单元。常规平面网格参数选取以计算机运算能力为标准,本申请证实复杂构造带构造建模需要充分考虑断层组合关系及断层影响界限,明确最低序级断层界限对复杂构造带平面网格大小及方位具有明显控制作用,为此提出以“不同序级断层组合分析”为确定层面网格参数的原则。此次提出的平面模型网格大小及方位受最低序级断层控制的观点,奠定了复杂构造带构造建模的技术基础。
可选地,本申请提供的这种复杂逆冲推覆构造三维构造建模方法从地震、单井数据出发,结合构造正反演技术利用非结构化网格建摸方法精细雕刻出构造模型,提高了构造模型的精度。
可选地,本申请通过应用提出了常规三维构造建模容易忽视的平面网格方位、大小、平面差值算法优选等参数及执行操作,促进了逆冲复杂构造地质建模的技术发展。
作为本申请另一种能够实现的方式,本申请还提供有如下应用场景。请参见图4,图4是本申请实施例提供的一种复杂推覆构造三维构造建模技术的实现示意图。在图4中,基础研究数据库400的数据来源可以包括岩心数据411,测井数据412、地震数据413和历史数据414。需要说明的是,历史数据414又可称之为前人成果。
在图4中,计算设备能够按照复杂构造带构造样式、断层精细划分和成因分析技术三个流程分别从基础研究数据库400中提取数据进行分析。
在复杂构造带构造样式流程中,包括构造样式42a1,模型正演42a2,构造形态42a3和构造分类42a4四个步骤。
在断层精细划分流程中,包括断层划分42b1,期次识别42b2,组合样式42b3和断裂级别42b4四个步骤。
在成因分析技术流程中,包括成因分析42c1,组合样式42c2、应力分析42c3和综合分析42c4四个步骤。
在基础研究数据库400中的数据经过上述复杂构造带构造样式、断层精细划分和成因分析技术三个流程的处理后,能够被后续地质模型参数的处理步骤处理,分别为:
步骤430,进行地质模型参数定量分析。
步骤440,通过断层处理提取非结构网格法网格模型。
步骤450,构建角点网络构造模型。
步骤460,建立复杂构造模型。
综上所述,在本申请提供的一种实施例中,通过“模型正反演、组合定量分析、非结构化网格建立”为核心的复杂推覆构造三维构造建模技术方法,能够较为精确表征了复杂逆冲推覆构造带不同构造应力造成的不同圈闭类型的平面分布趋势和垂向分布趋势,且可三维可视化,提高了构造认识的精度,为气藏井位部署及后期开发工作奠定了基础。
在另一种逆冲推覆构造带构造的建模方法中,本实施例提供一种方案,该方案能够克服:构造建模方法在稀疏网井条件下,不能充分表征逆冲推覆构造的空间展布特征及圈闭特征的问题。以C气藏为例进行介绍,其中,C气藏构造组合复杂、断层规模较大,断层规模大小不一又相互交织,是一套复杂的断裂***。整体褶皱变形强烈、断裂复杂,形成一系列呈北东南西向展布、***行的褶皱断高构造群。在该构造中,发育8个局部高点。可选地,C气藏的圈闭面积可达41.12平方千米,断层具有47条;构造发育断层类型以逆断层为主,断层样式繁多,断距处于125m至635米的区间。
步骤(1),采集单井的井信息数据、单井测井数据、井点分层数据以及采集微相数据,对所有单井数据进行校正,完成取奇异值、岩心校正等工作,并确定C气藏地层和断层参数。
步骤(2)利用步骤(1)中的井信息数据、单井测井数据和井点分层数据,运用三维地震与井精细标定技术,对C气藏所在的构造区带进行了规律分析及识别划分,确定构造产状为北东向、近东西向、近南北向;构造区带自西向东依次划分为C1逆冲推覆带、C2隐伏前缘带、C3坳陷带。发育多个潜伏构造或封闭于断层的潜伏断块、断鼻,总的趋势仍显示为北西高南东低的特点。
请参见图5,图5是本申请提供的一种C气藏所在的构造区带的沉积微相模型剖面图。在图5中,510是C1逆冲推覆带,520是C2隐伏前缘带,530是C3坳陷带。
步骤(3)在步骤(1)基础上,采用正反演相互引证思路,采用2D声学波动方程建立构造模型开展正演,形成不同类型断块构造样式放大倍数的量版,开展断块构造样式研究,指导复杂逆冲构造带构造解释分析。在相对精细可靠的构造解释基础上,结合区带构造应力模拟分析研究,将构造带构造背景类型划分为五类并进行定量评价(***背斜、滚动背斜、断鼻构造、上升断块、下降断块),指导后续构造模型建立。
请参见图6,图6是本申请提供的一种2D声学波动方程正演模拟Y型断层模板。其中,界面610表示网格10*10,主频40Hz正演模型。界面620表示网格10*10,主频40Hz正演剖面图。界面630表示Y型断层模板。
请参见图7,图7是本申请提供的一种逆冲推覆构造地区三维区内构造类型特征图。在图7中,公开了***背斜、滚动背斜、断鼻、上升断块和下降断块几种构造类型的典型模式、构造圈闭和地震剖面。
步骤(4)在步骤3)构造带构造背景类型划分为五类并进行定量评价的基础上,对断层组合关系进行精细分析,确定复杂构造带内断层的平面组合样式及剖面组合样式。确定复杂构造带发育四种断层组合样式(多组“Y”型、阶梯状、地堑式、地垒式),三级断层模式:一级断裂控制断块发育边界,二级断裂控制区域伸展,三级断裂导致断块复杂化。
步骤(5)在步骤(4)构造定量评价及断层模式分析的基础上,确定不同序级断层在平面分布特征及影响界线。高序级断层(一级、二级)断裂组合关系清晰,断点清楚,构造模型以一级、二级构造为主,控制断层截切及网格连接。低序级断层(三级)为区域应力释放产物,位于高序级断层末端,断层组合关系、断点较复杂,为派生断层。
请参见图8,图8是本申请提供的一种复杂构造带inline1528线不同断层级别、组合关系时深转换地震剖面图。其中,高序级断层810和高序级断层820断点清楚;低序级断层830的断点较复杂。
在构造圈闭类型分析清楚的基础上,开展构造建模研究。通过技术比选,采用非结构化网格算法,确定平面网格划分参数。可选地,平面上网格尺度以最低序级断层集中发育区设置1~2个网格为下限。可选地,复杂推翻构造最低序级断层间距为55米,确定平面网格下限30米。
纵向划分精度利用小波分解技术,建立高精度层序地层格架。综合预测误差滤波分析技术,确定纵向网格划分层序为2个四级层序4个小层。结合储层发育位置,本次建模构造纵向划分为4个zone,储层主要分布在1、2内。利用小波分解的d3曲线将目标组分为8个五级层序。结合求取旋回最小差异厚度,确定纵向网格最小尺寸0.82m。
请参见图9,图9是本申请实施例提供的逆冲复杂构造带非结构化网格建模技术反映复杂构造的示意图。其中,界面910表示角点网格的示意图,界面920表示非机构化网格的示意图。
请参见图10,图10是本申请实施例提供的一种不同低序级断层平面网格分布图。
请参见图11,图11是本申请实施例提供的一种旋回差异算法确定单井最小垂厚对比图。其中,界面1110用于指示旋回最小差异垂厚0.69米,垂深(7452.70-7453.39)的图示。界面1120用于指示旋回最小差异垂厚1.16米,垂深(7216.98~7218.14)的图示。界面1130用于指示旋回最小差异垂厚0.61米,垂深(7160.84~7161.45)的图示。
开展克里金算法、最小曲率法、收敛法多层面差值算法优选,最终优选收敛法建立层面网格模型。克里金算法、最小曲率法、收敛法等方法由于算法不同,精度存在一定差异,在常规三维地质建模中,差异基本可以忽略不计。本申请发现,对于构造复杂、逆冲推覆构造中,部分算法并不适用。克里金算法在应用过程中,出现井间差值变化平缓,但模拟构造面,尤其是逆冲复杂构造部位,层面变形严重,基本无法正确描述构造特征。通过比对优选,克里金算法、最小曲率法与最后选用的收敛法,两者差距达到7.77米。
请参见图12,图12是本申请实施例提供的一种逆冲推覆构造层面中的算法模拟适应性分析图。其中,图示1210表示克里金算法,用于指示井间差值变化相对平缓,模拟构造面边部翘起。图示1220表示最小曲率法,用于指示井间差值变化相对平缓。图示1230表示收敛法,用于指示尽可能降低井间模拟险。
利用三维网格化(英文:Pillar Gridding)构造建模技术,构建更加精确的复杂断层模型。控制交切断层、逆断层、铲状断层断面的闭合及断层尖灭点,优化断层截切关系,建立层面与断层有效接触关系。
构造模型的建立:根据已经建立的层面模型和zone模型,对输入的井点、断层线等数据进行截断变换,去除奇异值,分小层、分层面得到各构造的组合关系,最终建立相对准确的逆冲复杂构造模型。
模型验证:通过建立的逆冲复杂构造模型可以看出,D1构造为上升断块型构造、D2构造为***背斜型构造、D3构造为滚动背斜型。
D1构造、D2构造后经钻井资料证实,两者确实为不连通的构造单元,两个构造由于构造应力、断层组合方式不同,裂缝发育程度出现明显差异。D1构造受到双向挤压,断裂发育,测试气井井均测试产量80万方以上;D2***背斜型构造,断层滑动特征明显,应力释放严重,断裂不是太发育,气井测试产量仅为30万方。因此,所建立的逆冲复杂构造带构造模型具有可信度,可以用于指导气藏的圈闭评价以及后续开发。
请参见图13,图13是本申请实施例提供的一种复杂逆冲推覆构造带三维构造模型图。
下述为本申请装置实施例,可以用于执行本申请方法实施例。对于本申请装置实施例中未披露的细节,请参照本申请方法实施例。
请参考图14,图14示出了本申请一个示例性实施例提供的逆冲推覆构造带构造的建模装置的结构框图。该逆冲推覆构造带构造的建模装置可以通过软件、硬件或者两者的结合实现成为计算设备的全部或一部分。该装置包括:
第一建模模块1410,用于根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;
量版获取模块1420,用于对所述目标构造模型开展正演,获取目标量版,所述目标量版包括至少两种类型且放大n倍的断块构造样式,n大于1;
第二建模模块1430,用于根据所述目标量版,对至少两个不同序级的断层进行展布分析,建立第一层面网格模型;
第三建模模块1440,用于根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
在一个可选的实施例中,所述装置还包括勘探数据获取模块和勘探数据分析模块。勘探数据获取模块,用于获取目标勘探数据,所述目标勘探数据是指示所述逆冲推覆构造带构造在被勘探时收集的地质数据;勘探数据分析模块,用于分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据。
在一个可选的实施例中,所述勘探数据分析模块,用于通过三维地震和井精细标定技术分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据。
在一个可选的实施例中,所述勘探数据获取模块,用于获取目标单井的单井数据,所述目标单井是逆冲推覆构造带中的单井。所述勘探数据分析模块,用于通过小波分解技术处理所述单井数据,获取所述逆冲推覆构造带的所述形态数据,其中,所述形态数据包括气藏地界数据,以及,所述目标单井的断层参数。
在一个可选的实施例中,所述第二建模模块1430,用于根据所述目标量版,确定所述逆冲推覆构造带构造的断层影响界限和断层规模;根据所述断层影响界限和所述断层规模,对至少两个不同序级的断层进行展布分析,建立所述第一层面网格模型。
在一个可选的实施例中,所述第三建模模块1440,用于根据目标多层面差值算法对所述第一层面网格模型进行优化,得到中间层面网格模型;对所述中间层面网格模型进行井分层数据内插,得到所述第二层面网格模型。
在一个可选的实施例中,所述装置还包括筛选模块,用于从候选多层面差值算法中确定出目标多层面差值算法,所述候选多层面差值算法包括克里金算法、最小曲率法和收敛法中至少两种算法。
本申请实施例还提供了一种计算机可读介质,该计算机可读介质存储有至少一条指令,所述至少一条指令由所述处理器加载并执行以实现如上各个实施例所述的逆冲推覆构造带构造的建模方法。
需要说明的是:上述实施例提供的逆冲推覆构造带构造的建模装置在执行逆冲推覆构造带构造的建模方法时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的逆冲推覆构造带构造的建模装置与逆冲推覆构造带构造的建模方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
上述本申请实施例序号仅仅为了描述,不代表实施例的优劣。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本申请的能够实现的示例性的实施例,并不用以限制本申请,凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (9)
1.一种逆冲推覆构造带构造的建模方法,其特征在于,所述方法包括:
根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;
对所述目标构造模型开展正演,获取目标量版,所述目标量版包括指示至少两种类型且放大n倍的断块构造样式,n大于1;
对所述目标量版进行断层级别定量划分,确定所述逆冲推覆构造带的断层组合关系;
对所述断层组合关系进行分析,确定复杂构造带内断层的剖面组合样式以及三级断层模式,所述三级断层模式包括一级断裂、二级断裂以及三级断裂,其中,所述一级断裂控制断块发育边界,所述二级断裂控制区域伸展,所述三级断裂控制断块的复杂程度;
根据所述剖面组合样式以及所述三级断层模式,确定至少两个不同序级的断层的平面分布特征以及影响界线,其中,所述至少两个不同序级中包括以所述一级断裂和所述二级断裂在内的高序级断层以及以所述三级断裂在内的低序级断层,所述低序级断层位于所述高序级断层末端,所述低序级断层为派生断层,由高序级断层控制断层截切以及平面网格的连接;
根据所述平面分布特征以及影响界线,利用非结构网格法建模技术,确定所述平面网格的划分参数,所述平面网格的尺度以最低序级的断层集中发育区设置1-2个网格单元为下限;
利用小波分解技术以及预测误差滤波分析技术,结合储层发育位置对所述至少两个不同序级的断层进行纵向方向的划分,得到划分后的至少两个不同序级的断层,所述划分后的至少两个不同序级的断层的层序精度高于所述至少两个不同序级的断层的层序精度;
基于所述平面网格的划分参数以及所述划分后的至少两个不同序级的断层,建立第一层面网格模型;
根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
2.根据权利要求1所述的方法,其特征在于,在所述根据所述形态数据和二维声学波动方程建立目标构造模型前,所述方法还包括:
获取目标勘探数据,所述目标勘探数据是指示所述逆冲推覆构造带构造在被勘探时收集的地质数据;
分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据。
3.根据权利要求2所述的方法,其特征在于,所述分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据,包括:
通过三维地震和井精细标定技术分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据。
4.根据权利要求3所述的方法,其特征在于,所述获取目标勘探数据,包括:
获取目标单井的单井数据,所述目标单井是逆冲推覆构造带中的单井;
所述通过三维地震和井精细标定技术分析所述目标勘探数据,获取所述逆冲推覆构造带构造的形态数据,包括:
通过小波分解技术处理所述单井数据,获取所述逆冲推覆构造带的所述形态数据,其中,所述形态数据包括气藏地界数据,以及,所述目标单井的断层参数。
5.根据权利要求1至4任一所述的方法,其特征在于,所述根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,包括:
根据目标多层面差值算法对所述第一层面网格模型进行优化,得到中间层面网格模型;
对所述中间层面网格模型进行井分层数据内插,得到所述第二层面网格模型。
6.根据权利要求5所述的方法,其特征在于,所述方法还包括:
从候选多层面差值算法中确定出目标多层面差值算法,所述候选多层面差值算法包括克里金算法、最小曲率法和收敛法中至少两种算法。
7.一种逆冲推覆构造带构造的建模装置,其特征在于,所述装置包括:
第一建模模块,用于根据形态数据和二维声学波动方程建立目标构造模型,所述形态数据用于指示逆冲推覆构造带构造的各地层的形态以及地理位置分布;
量版获取模块,用于对所述目标构造模型开展正演,获取目标量版,所述目标量版包括至少两种类型且放大n倍的断块构造样式,n大于1;
第二建模模块,用于对所述目标量版进行断层级别定量划分,确定所述逆冲推覆构造带的断层组合关系;所述断层组合关系进行分析,确定复杂构造带内断层的剖面组合样式以及三级断层模式,所述三级断层模式包括一级断裂、二级断裂以及三级断裂,其中,所述一级断裂控制断块发育边界,所述二级断裂控制区域伸展,所述三级断裂控制断块的复杂程度;根据所述剖面组合样式以及所述三级断层模式,确定至少两个不同序级的断层的平面分布特征以及影响界线,其中,所述至少两个不同序级中包括以所述一级断裂和所述二级断裂在内的高序级断层以及以所述三级断裂在内的低序级断层,所述低序级断层位于所述高序级断层末端,所述低序级断层为派生断层,由高序级断层控制断层截切以及平面网格的连接;根据所述平面分布特征以及影响界线,利用非结构网格法建模技术,确定所述平面网格的划分参数,所述平面网格的尺度以最低序级的断层集中发育区设置1-2个网格单元为下限;利用小波分解技术以及预测误差滤波分析技术,结合储层发育位置对所述至少两个不同序级的断层进行纵向方向的划分,得到划分后的至少两个不同序级的断层,所述划分后的至少两个不同序级的断层的层序精度高于所述至少两个不同序级的断层的层序精度;基于所述平面网格的划分参数以及所述划分后的至少两个不同序级的断层,建立第一层面网格模型;
第三建模模块,用于根据目标多层面差值算法对所述第一层面网格模型进行优化,建立第二层面网格模型,所述第二层面网格模型的精度高于所述第一层面网格模型,且所述第二层面网格模型是所述逆冲推覆构造带构造的三维模型。
8.一种计算设备,其特征在于,所述计算设备包括处理器和存储器,所述存储器中存储有至少一条指令,所述指令由所述处理器加载并执行以实现如权利要求1至6任一所述的逆冲推覆构造带构造的建模方法。
9.一种计算机可读存储介质,其特征在于,所述存储介质中存储有至少一条指令,所述指令由处理器加载并执行以实现如权利要求1至6任一所述的逆冲推覆构造带构造的建模方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910292415.3A CN111815769B (zh) | 2019-04-12 | 2019-04-12 | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910292415.3A CN111815769B (zh) | 2019-04-12 | 2019-04-12 | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111815769A CN111815769A (zh) | 2020-10-23 |
CN111815769B true CN111815769B (zh) | 2024-04-30 |
Family
ID=72844617
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910292415.3A Active CN111815769B (zh) | 2019-04-12 | 2019-04-12 | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111815769B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117348080B (zh) * | 2023-10-27 | 2024-07-02 | 西南石油大学 | 一种基于构造正演模拟判断逆冲断层发育程度的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236088A (zh) * | 2013-04-16 | 2013-08-07 | 北京金阳普泰石油技术股份有限公司 | 基于复合网格的三维地质体构造建模方法与*** |
CN106501856A (zh) * | 2016-10-20 | 2017-03-15 | 东北石油大学 | 亚地震断层定量预测方法 |
CN106772586A (zh) * | 2017-02-20 | 2017-05-31 | 长江大学 | 一种基于地震信号奇异性的隐蔽性断裂检测方法 |
CN106855636A (zh) * | 2017-03-23 | 2017-06-16 | 西南石油大学 | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 |
CN109035411A (zh) * | 2018-07-30 | 2018-12-18 | 中国石油化工股份有限公司 | 基于正演模拟的三维地层角点网格体元建模方法 |
CN109164487A (zh) * | 2018-09-07 | 2019-01-08 | 中国石油化工股份有限公司 | 一种基于模型建立平均速度场的方法及精细成图的方法 |
-
2019
- 2019-04-12 CN CN201910292415.3A patent/CN111815769B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236088A (zh) * | 2013-04-16 | 2013-08-07 | 北京金阳普泰石油技术股份有限公司 | 基于复合网格的三维地质体构造建模方法与*** |
CN106501856A (zh) * | 2016-10-20 | 2017-03-15 | 东北石油大学 | 亚地震断层定量预测方法 |
CN106772586A (zh) * | 2017-02-20 | 2017-05-31 | 长江大学 | 一种基于地震信号奇异性的隐蔽性断裂检测方法 |
CN106855636A (zh) * | 2017-03-23 | 2017-06-16 | 西南石油大学 | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 |
CN109035411A (zh) * | 2018-07-30 | 2018-12-18 | 中国石油化工股份有限公司 | 基于正演模拟的三维地层角点网格体元建模方法 |
CN109164487A (zh) * | 2018-09-07 | 2019-01-08 | 中国石油化工股份有限公司 | 一种基于模型建立平均速度场的方法及精细成图的方法 |
Non-Patent Citations (1)
Title |
---|
胜利油田高89地区滩坝砂油藏CO_2驱地质建模技术研究;杨明林;;内蒙古石油化工(第06期);第116页-118页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111815769A (zh) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2486336C2 (ru) | Способы имитации разрыва пласта-коллектора и его оценки и считываемый компьютером носитель | |
RU2669948C2 (ru) | Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности | |
US9684084B2 (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
US9026416B2 (en) | System and method of grid generation for discrete fracture modeling | |
AU2009260453B2 (en) | Heterogeneous earth models for a reservoir field | |
US8855986B2 (en) | Iterative method and system to construct robust proxy models for reservoir simulation | |
CN103975341B (zh) | 基于4d饱和度模型和仿真模型的储层建模 | |
CN110308487B (zh) | 一种断溶体型油藏定量表征方法 | |
US10620339B2 (en) | Static earth model calibration methods and systems using tortuosity evaluations | |
CN107402176B (zh) | 裂缝孔隙度的确定方法和装置 | |
CA2869825C (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
KR20160024232A (ko) | 지반정보에 기초한 3차원 공간 모델링 방법 | |
Mahjour et al. | Developing a workflow to represent fractured carbonate reservoirs for simulation models under uncertainties based on flow unit concept | |
CN110705168A (zh) | 构造应力场的模拟方法 | |
CN111815769B (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
CN116449423A (zh) | 一种辫状河单砂体构型分析方法、装置、设备及存储介质 | |
CN115758792A (zh) | 一种基于数字数值一体化的地质灾害评估方法及装置 | |
Regli et al. | GEOSSAV: a simulation tool for subsurface applications | |
CN115469361B (zh) | 一种碎屑岩地层三维地质建模方法 | |
Li et al. | Uncertainty evaluation method and application of reservoir Lithofacies proportion under the condition of few wells | |
CN117852416B (zh) | 基于数字化地质模型的多模态注浆预控分析方法及*** | |
US20230259662A1 (en) | Modeling a karst formation for a wellbore operation | |
CN118260982A (zh) | 三维地应力场特征确定方法及装置 | |
Correia et al. | Using pilot wells to integrate geological modelling and history matching: applied to the Norne Benchmark case | |
WANG et al. | Reserves evaluation for new investment projects in overseas oil and gas field development [J] |
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 |