CN113219529B - 边界约束的三维构造地质模型构建方法 - Google Patents
边界约束的三维构造地质模型构建方法 Download PDFInfo
- Publication number
- CN113219529B CN113219529B CN202110481452.6A CN202110481452A CN113219529B CN 113219529 B CN113219529 B CN 113219529B CN 202110481452 A CN202110481452 A CN 202110481452A CN 113219529 B CN113219529 B CN 113219529B
- Authority
- CN
- China
- Prior art keywords
- intersection
- fault
- model
- data
- boundary
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000009792 diffusion process Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000005457 optimization Methods 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000007789 sealing Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims 1
- 230000005012 migration Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 238000012360 testing method Methods 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/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/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
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种边界约束的三维构造地质模型构建方法,包括以下步骤:S1、进行交线估计,包括以下子步骤:S11、对原始构造解释数据进行平滑去噪处理;S12、提取层位面上的极值点;S13、根据极值点对层面数据进行分段,找到参与交线估计的有效层位数据;S14、在解释剖面上估计曲面的交点;S15、将交点拟合为交线;S2、用一个交点偏移在二维平面的扩散模型对原始数据的调整量进行建模。本发明根据原始构造解释数据通过计算估计得到构造模型地质曲面的边界约束条件,再在约束下进行地质曲面重构和地质体生成,提高了构造建模方法的有效性和效率。并且有曲面边界约束建模的机制,就能使得建模过程更鲁棒、更准确,减少和代替部分人机交互。
Description
技术领域
本发明涉及一种边界约束的三维构造地质模型构建方法。
背景技术
复杂的地质过程在地下空间的多个尺度上产生明显的结构异质性,对地下空间结构的认识直接影响着人类社会的方方面面:矿产资源评估,油气勘探与生产,地热资源勘探,地质灾害评估,水文地质,岩土工程等。除了通过钻井和挖掘隧道之外不可能直接进入地下,所以大多数对地质构造的认知必须间接来自对地球物理数据的处理。三维构造地质建模作为能够直观展示地质构造形态和分布的关键技术,已经成为了地球信息科学和工业界共同关注的研究热点之一。
构造建模就是以地震解释得到的地质曲面(层位面和断层面)和钻井数据为基础,采用一定的地质曲面重建算法,在三维空间中生成空间曲面,再确定空间位置不同的地质曲面之间的拓扑关系。然后以断层面和层位面为模型的内边界进行地质块体的划分,生成以封闭地质体为基本单位的模型。一方面,随着矿产和油气等资源的地球物理勘探向地质条件复杂、埋藏较深、地层古老的区域发展,三维地震资料获取的成本和难度加大,导致原始地震数据质量偏低;另一方面,钻井资料虽然是最直接和最可靠的样本来源,但由于钻井成本的限制,导致样本数量偏少,分布不均匀。所以作为构造建模输入的构造解释数据通常包含了不可忽视的不确定性(不准确、缺失、不一致性)。构造建模的难点就是基于不确定性数据源,建立合理且准确的模型。传统建模方法适合高数据密度高数据质量的环境,例如在存在大量钻井的区域。但实际应用很少能满足这样的条件,得到的模型通常不能完全满足地质规则或构造解释专家对模型的期望,还需要对模型进行检查并对不合理的局部进行手动调整,每一次调整操作后都需要再完成一遍工作流程来更新模型。其根本原因在于构造解释数据同时记录了地质对象的概念信息和地质模型的可视信息,但传统构造建模方法仅使用了其中的可视信息。单纯可视信息不足以准确描述地质构造,建模过程中又缺乏能减少不确定性的其他约束条件。
在输入的构造解释数据和钻井数据质量不变的情况下(不在本发明讨论范围内),如果计算机能够在模型建立前从现有输入中挖掘出额外的输入信息,并且有这些额外信息约束建模的机制,就能使得建模过程更鲁棒、更准确,减少和代替部分人机交互。
地质构造建模是在将地质,测井,地球物理资料和各种解释结果或者概念模型进行综合分析的基础上,利用计算机图形技术,生成能够反映地质构造形态和关系的三维随机模型。这样建立的模型汇总了各种信息和解释结果。多种三维地质建模方法相继被提出,如:基于钻孔的建模,基于体素的建模,基于垂直剖面的建模,基于曲面的建模,基于多源数据的建模,等等。在本研究中我们提出的是一种基于曲面的建模方法,最终结果是地质构造的三维格架模型,格架是由地质曲面(主要为断层和层位面)及其拓扑关系组成的。地质曲面作为模型的内边界将模型空间分割,形成多个封闭的地质体。
图1给出了通用的基于曲面的建模工作流程。在实际应用中,由于地震数据通常在断层附近的区域信噪比较低,这对构造解释造成了困难,所以原始的构造解释数据带有较大不确定性,例如图2(a)所示的缺失和冲突情况,图2(b)所示的样本少和样本分布不均的情况。
发明内容
本发明的目的在于克服现有技术的不足,提供一种根据原始构造解释数据通过计算估计得到构造模型地质曲面的边界约束条件,再在约束下进行地质曲面重构和地质体生成,提高了构造建模方法的有效性和效率,使得建模过程更鲁棒、更准确的边界约束的三维构造地质模型构建方法。
本发明的目的是通过以下技术方案来实现的:边界约束的三维构造地质模型构建方法,包括以下步骤:
S1、进行交线估计,包括以下子步骤:
S11、对原始构造解释数据进行平滑去噪处理;
S12、提取层位面上的极值点;
S13、根据极值点对层面数据进行分段,找到参与交线估计的有效层位数据;
S14、在解释剖面上估计曲面的交点;
S15、将交点拟合为交线;
S2、用一个交点偏移在二维平面的扩散模型对原始数据进行调整;
S3、在已知全部的交线后,用调整后的原始数据通过克里金插值重构层位面和断层,重构后的曲面自动封闭,形成地质体,并最终得到封闭的完整模型。
进一步地,所述步骤S14具体实现方法为:设原始构造解释数据有n个解释剖面,第i个剖面上的交点为si;有效层位数据为其中表示层位解数据在第i个剖面上的第x个点,1≤x≤α,α表示i个剖面上数据总数量;对Hi进行拟合得到的曲线Pi作为层位H在第i个剖面上的外推表示;断层数据为其中表示断层解释数据在第i个剖面上的第y个点,1≤y≤β,β表示断层数据的总数量;对Fi进行拟合得到的曲线为Qi;模型边界为B,表示模型内部;根据知识图谱得到当前层位H的相邻层位分别为H’和H”,对应F在断层同一侧的交点分别为s′i和s″i;
建立以下优化目标:
(1)交点与层位面趋势的吻合度:
f1=dist(si,Pi)=min(||si-m||2)
其中m表示曲线Pi上任意一点,其坐标为(mx,my,mz);
(2)交点可靠性:
(3)上下盘交点趋势一致性:
优化约束条件为:
进一步地,所述步骤S15具体实现方法为:交线表示为k+1阶B样条曲线的形式,交点si是B样条曲线的控制点:
进一步地,所述步骤S2具体实现方法为:用一个交点偏移在二维平面的扩散模型对原始数据的调整量进行建模,其基本形式为:
设[xmin xmax]×[ymin ymax]为一个交点的偏移影响范围,
(xmin,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax)分别表示交点的偏移影响范的4个角点;其中w表示偏移量,设x方向为交线在这一点处的切线方向,y方向为垂直交线方向,(x,y)为空间二维坐标;w(x,y,t)表示(x,y)位置在t时刻的偏移量,
w(xmin,y,t)、w(xmax,y,t)、w(x,ymin,t)、w(x,ymin,t)分别表示交点偏移影响范围边界处在t时刻的偏移量;vx,vy为扩散系数;设(sx,sy)为交点坐标,f(sx,sy)为偏移模型的初始条件,即交点的偏移情况,对应w(sx,sy,0);gxmin(y,t)、gxmax(y,t)、gymin(x,t)、gymax(x,t)表示偏移模型的边界条件,通常偏移量在边界处应该为0,即gxmin(y,t)=gxmax(y,t)=gymin(x,t)=gymax(x,t)≡0;t在计算偏移量时可以为任意常数。
本发明的有益效果是:本发明设计了一种从原始构造解释数据中估计曲面边界条件(即曲面的交线)并用于约束建模过程的方法。首先根据原始构造解释数据通过计算估计得到构造模型地质曲面的边界约束条件,再在约束下进行地质曲面重构和地质体生成。本发明提出在构造建模方法中加入地质曲面边界约束条件的引导,提高了构造建模方法的有效性和效率。并且有曲面边界约束建模的机制,就能使得建模过程更鲁棒、更准确,减少和代替部分人机交互。
附图说明
图1为通用的基于曲面的构造建模工作流程;
图2为构造解释数据的不确定性例子;
图3为本发明的原始解释数据示意图;
图4为剖面中的构造特征示意图;
图5为本发明的三维构造地质模型构建方法的流程图;
图6为构造解释数据的平滑处理示意图;
图7为层位面的极值点示意图;
图8为图3所示的层位和断层在ε=0时的构造特征离群情况;
图9为可靠性对交线估计的影响图;
图10为本发明具体实施例构造地质模型原始数据;
图11为不同剖面上的曲面交点(T4与F5的第12~15剖面);
图12为曲面间的交线(T4与F5);
图13为重构层位面;
图14为重构断层面;
图15为最终得到封闭的完整模型。
具体实施方式
一个复杂地质模型可以细分成多个相对简单的地质体,当实现了各简单地质体闭合构建时,整个三维地质模型也就形成了。地质体是由多个地质曲面所封闭而成的,而地质曲面又是由模型边界线、地层尖灭线、断层线等边界线组成的,因此可以说地质体也是由边界线是作为轮廓支撑起来的。所以确定了边界线也就确定了地质体之间的接触关系和地质曲面的相交关系。传统构造建模方法首先独立地重构各个曲面,然后基于曲面切割算法计算交线。在这种不考虑曲面相互影响的情况下不能保证得到的交线符合约束。所以,在我们的建模流程中,通过在曲面重构前估计曲面交线来实现对构造建模的约束。
估计交线的本质是一个优化问题,其约束条件是曲面空间关系,其优化目标之一就是交线的“可靠性”。原始解释数据的质量通常参差不齐,如图3所示,图中,(a)为层位面解释数据例子,第13剖面断距发生突变;(b)为断层解释数据例子,第7剖面断层倾角发生突变。不考虑原始数据的质量而直接计算的交线是不可靠的。所以我们首先提出交线的可靠性评估模型作为交线估计的优化目标之一。我们首先在各个解释剖面上估计曲面交点,得到交点后可以进一步计算出一系列构造特征:地层断距(d),垂直断距(v),水平断距(h),断层倾角(ω),断层两侧的层位倾角(α,β),如图4所示。我们知道构造特征在一系列有序剖面中的变化是符合某一规律的,如果存在偏离规律的构造特征,就说明这一剖面附近的交线可靠性较低,可靠性的数值通过构造特征的离群度来度量。记构造特征X的观测值为X={xi}(i=1,2,...,n),有n个解释剖面,xi表示X在第i个剖面的观测值。将X拟合为模型q,纵坐标为所有观测值,横坐标为剖面间的距离。q的目标形式通过观察给出。
定义1:观测点xi的离群因子pi:观测点xi到拟合模型的偏差距离,pi=min(||xi-q||2)。
定义4:离群点:有pi>δ*AP的观测点xi。
定义5:观测点xi的离群度εi:ε的取值以0为中心,以l(l>0)为步长向正负方向改变ε取值(0+m*l),随着ε减小,越来越多的观测点会被判定为离群点,设ε=εmin时所有观测点都被判定为离群点(ε=εmin+l时离群点数<n),ε=εmax时没有离群点(ε=εmax-l时离群点数>0)。观测点xi首次被判定为离群点对应的ε值即为xi的离群度εi(εmin<εi≤εmax)。
离群度εi就是交线在第i剖面附近关于特征X的可靠性度量,εi越大说明第i剖面的交点关于特征X越不可靠。组合多种构造特征的离群度就可以得到交点的可靠性:
模型f表示不同种类构造特征间的联系,可以由专家根据经验给出。
有了对交线的评价标准后,接下来我们提出交线估计的步骤,并建立交线估计的多目标优化模型。下面结合附图进一步说明本发明的技术方案。
如图5所示,本发明的边界约束的三维构造地质模型构建方法,包括以下步骤:
S1、进行交线估计,包括以下子步骤:
S11、对原始构造解释数据进行平滑去噪处理,使用滑动窗口局部回归的方式,如图6所示;
S12、提取层位面上的极值点,如图7所示;
S13、根据极值点对层面数据进行分段,找到参与交线估计的有效层位数据;我们认为离断层最近的极值点到断层的一段层位数据对交点估计的贡献最大,引入其他位置的层位数据反而会对交点的估计产生负面影响;
S14、在解释剖面上估计曲面的交点;
具体实现方法为:设原始构造解释数据有n个解释剖面,第i个剖面上的交点为si;有效层位数据为其中表示层位解数据在第i个剖面上的第x个点,1≤x≤α,α表示i个剖面上数据总数量;对Hi进行拟合得到的曲线Pi作为层位H在第i个剖面上的外推表示;断层数据为其中表示断层解释数据在第i个剖面上的第y个点,1≤y≤β,β表示断层数据的总数量;对Fi进行拟合得到的曲线为Qi;模型边界为B,表示模型内部;根据知识图谱得到当前层位H的相邻层位分别为H’和H”,对应F在断层同一侧的交点分别为s′i和s″i;
建立以下优化目标:
(1)交点与层位面趋势的吻合度:
f1=dist(si,Pi)=min(||si-m||2)
其中m表示曲线Pi上任意一点,其坐标为(mx,my,mz);
(2)交点可靠性:
(3)上下盘交点趋势一致性:
优化约束条件为:
S15、将交点拟合为交线;具体实现方法为:交线表示为k+1阶B样条曲线的形式,交点si是B样条曲线的控制点:
接下来验证加入交线可靠性度量对交线估计的影响。在不考虑可靠性的情况下直接用原始数据对图3中所示的层位面和断层的交点进行估计,构造特征的离群情况如图8所示,基本符合对原始数据的直接观察。图9中可以看到可靠性对交线估计的影响。色图表示直接估计的交点和交线的可靠性,(a)表示考虑可靠性时的交点;(b)表示考虑可靠性时的交线。可以看到可靠性越低的剖面附近的交点和交线相对不考虑可靠性的交点和交线偏离越多,说明可靠性是交线估计的有效约束。
S2、用一个交点偏移在二维平面的扩散模型对原始数据进行调整;
我们证明了在构造特征变化规律的约束下得到的交线相比纯数据驱动得到的交线出现了偏移。如果不对原始层位数据进行调整,那么在重构模型时将会面临原始数据与边界线冲突的问题。并且交线的偏移本身也反应了原始数据质量,交线偏移越多原始数据质量越差。本发明用一个交点偏移在二维平面的扩散模型对原始数据的调整量进行建模,其基本形式为:
设[xmin xmax]×[ymin ymax]为一个交点的偏移影响范围,
(xmin,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax)分别表示交点的偏移影响范的4个角点;其中w表示偏移量,设x方向为交线在这一点处的切线方向,y方向为垂直交线方向,(x,y)为空间二维坐标;w(x,y,t)表示(x,y)位置在t时刻的偏移量,
w(xmin,y,t)、w(xmax,y,t)、w(x,ymin,t)、w(x,ymin,t)分别表示交点偏移影响范围边界处在t时刻的偏移量;vx,vy为扩散系数;设(sx,sy)为交点坐标,f(sx,sy)为偏移模型的初始条件,即交点的偏移情况,对应w(sx,sy,0);gxmin(y,t)、gxmax(y,t)、gymin(x,t)、gymax(x,t)表示偏移模型的边界条件,通常偏移量在边界处应该为0,即gxmin(y,t)=gxmax(y,t)=gymin(x,t)=gymax(x,t)≡0;t在计算偏移量时可以为任意常数。
S3、在已知全部的交线后,用调整后的原始数据通过克里金插值重构层位面和断层,重构后的曲面自动封闭,形成地质体,并最终得到封闭的完整模型。
实验测试环境:本发明所做的实验电脑配置为Inter i5-8500 CPU,8GRAM,64位Windows 10操作系统;采用的软件为MATLAB R2016。实际构造地质模型数据来自四川盆地某油气生产工区。其基础数据包含15个剖面,主要构造为4个地层面(T1,T2,T3,T4),5个逆断层(F1,F2,F3,F4,F5),其中F1断层没有穿过整个模型,如图10所示。
根据本发明所提出的方法依次在剖面上估计曲面的交点,如图11所示,图11为不同剖面上的曲面交点(T4与F5的第12~15剖面);并得到交线如图12所示,图12为曲面间的交线(T4与F5)。在已知全部的曲面边界线后重构层位面如图13所示,图中,(a)为T1面,(b)为T2面,(c)为T3面,(d)为T4面;重构断层面,如图14所示。重构后的曲面自动封闭,形成地质体,并最终得到封闭的完整模型,如图15所示。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (3)
1.边界约束的三维构造地质模型构建方法,其特征在于:包括以下步骤:
S1、进行交线估计,包括以下子步骤:
S11、对原始构造解释数据进行平滑去噪处理;
S12、提取层位面上的极值点;
S13、根据极值点对层面数据进行分段,找到参与交线估计的有效层位数据;
S14、在解释剖面上估计曲面的交点;具体实现方法为:设原始构造解释数据有n个解释剖面,第i个剖面上的交点为si;有效层位数据为其中表示层位解数据在第i个剖面上的第x个点,1≤x≤α,α表示i个剖面上数据总数量;对Hi进行拟合得到的曲线Pi作为层位H在第i个剖面上的外推表示;断层数据为其中表示断层解释数据在第i个剖面上的第y个点,1≤y≤β,β表示断层数据的总数量;对Fi进行拟合得到的曲线为Qi;模型边界为B,表示模型内部;根据知识图谱得到当前层位H的相邻层位分别为H’和H”,对应F在断层同一侧的交点分别为s′i和s″i;
建立以下优化目标:
(1)交点与层位面趋势的吻合度:
f1=dist(si,Pi)=min(||si-m||2)
其中m表示曲线Pi上任意一点,其坐标为(mx,my,mz);
(2)交点可靠性:
其中和分别表示第i个剖面断层上盘和下盘的交点,分别表示构造特征参数地层断距d,垂直断距v,水平断距h,断层倾角ω,断层两侧的层位倾角(α,β)在第i个剖面上的离群度;均表示可靠性,模型f表示不同种类构造特征间的联系;
(3)上下盘交点趋势一致性:
优化约束条件为:
S15、将交点拟合为交线;
S2、用一个交点偏移在二维平面的扩散模型对原始数据进行调整;
S3、在已知全部的交线后,用调整后的原始数据通过克里金插值重构层位面和断层,重构后的曲面自动封闭,形成地质体,并最终得到封闭的完整模型。
3.根据权利要求1所述的边界约束的三维构造地质模型构建方法,其特征在于:所述步骤S2具体实现方法为:用一个交点偏移在二维平面的扩散模型对原始数据的调整量进行建模,其基本形式为:
设[xminxmax]×[yminymax]为一个交点的偏移影响范围,(xmin,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax)分别表示交点的偏移影响范的4个角点;其中w表示偏移量,设x方向为交线在这一点处的切线方向,y方向为垂直交线方向,(x,y)为空间二维坐标;w(x,y,t)表示(x,y)位置在t时刻的偏移量,w(xmin,y,t)、w(xmax,y,t)、w(x,ymin,t)、w(x,ymin,t)分别表示交点偏移影响范围边界处在t时刻的偏移量;vx,vy为扩散系数;设(sx,sy)为交点坐标,f(sx,sy)为偏移模型的初始条件,即交点的偏移情况,对应w(sx,sy,0);gxmin(y,t)、gxmax(y,t)、gymin(x,t)、gymax(x,t)表示偏移模型的边界条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110481452.6A CN113219529B (zh) | 2021-04-30 | 2021-04-30 | 边界约束的三维构造地质模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110481452.6A CN113219529B (zh) | 2021-04-30 | 2021-04-30 | 边界约束的三维构造地质模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113219529A CN113219529A (zh) | 2021-08-06 |
CN113219529B true CN113219529B (zh) | 2022-06-14 |
Family
ID=77090477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110481452.6A Active CN113219529B (zh) | 2021-04-30 | 2021-04-30 | 边界约束的三维构造地质模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113219529B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116819615B (zh) * | 2023-08-30 | 2023-11-21 | 中国石油大学(华东) | 一种地震数据重建方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0681243A1 (en) * | 1994-04-28 | 1995-11-08 | Hewlett-Packard GmbH | Method for performing boolean operations on geometric objects in a computer-aided design system |
CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
CN107644228A (zh) * | 2017-09-21 | 2018-01-30 | 联想(北京)有限公司 | 图像处理方法 |
CN108776993A (zh) * | 2018-05-24 | 2018-11-09 | 北京建筑大学 | 带有孔洞的三维点云的建模方法及地下电缆工井建模方法 |
CN111047700A (zh) * | 2019-12-19 | 2020-04-21 | 山东科技大学 | 一种基于断层产生时序的地质构造三维建模方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2542874C1 (ru) * | 2013-10-15 | 2015-02-27 | Антон Александрович Колесниченко | Способ адаптивной разбивки трехмерной сетки на слои при создании трехмерных цифровых геологических моделей пластов |
CN104597496B (zh) * | 2015-01-30 | 2017-08-25 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种二维地震资料中速度的三维空间归位方法 |
CN109949415B (zh) * | 2019-02-26 | 2023-04-11 | 天衍智(北京)科技有限公司 | 一种三维地表与地质体模型拓扑一致建模的系统及方法 |
CN111796324B (zh) * | 2019-04-09 | 2023-02-10 | 中国石油天然气股份有限公司 | 地震全层位追踪方法及装置 |
-
2021
- 2021-04-30 CN CN202110481452.6A patent/CN113219529B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0681243A1 (en) * | 1994-04-28 | 1995-11-08 | Hewlett-Packard GmbH | Method for performing boolean operations on geometric objects in a computer-aided design system |
CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
CN107644228A (zh) * | 2017-09-21 | 2018-01-30 | 联想(北京)有限公司 | 图像处理方法 |
CN108776993A (zh) * | 2018-05-24 | 2018-11-09 | 北京建筑大学 | 带有孔洞的三维点云的建模方法及地下电缆工井建模方法 |
CN111047700A (zh) * | 2019-12-19 | 2020-04-21 | 山东科技大学 | 一种基于断层产生时序的地质构造三维建模方法及系统 |
Non-Patent Citations (1)
Title |
---|
复杂地质构造智能建模方法研究;鲁才;《中国地球科学联合学术年会 2020》;20200831;第1372页倒数第一段-第1373页第一段 * |
Also Published As
Publication number | Publication date |
---|---|
CN113219529A (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wellmann et al. | 3-D Structural geological models: Concepts, methods, and uncertainties | |
CA2725923C (en) | Heterogeneous earth models for a reservoir field | |
EP0981768B1 (en) | Method for determining barriers to reservoir flow | |
EP3217192B1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
Che et al. | Three-dimensional geological modeling of coal seams using weighted Kriging method and multi-source data | |
US20120271609A1 (en) | Methods and computing systems for hydrocarbon exploration | |
WO2021008630A1 (zh) | 圈闭断层的封闭性定量分析方法、装置及系统 | |
US10914861B2 (en) | Inversion-based workflow for consistent interpretation of nuclear density images in horizontal wells | |
WO2015050530A1 (en) | In-situ wellbore, core and cuttings information system | |
Manchuk et al. | Boundary modeling with moving least squares | |
US20220237891A1 (en) | Method and system for image-based reservoir property estimation using machine learning | |
CN112150582B (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
CN113219529B (zh) | 边界约束的三维构造地质模型构建方法 | |
US11261729B2 (en) | Methods and systems for determining integrity and operational boundaries of subterranean wells | |
WO2022159698A1 (en) | Method and system for image-based reservoir property estimation using machine learning | |
CN106443772B (zh) | 一种去底辟原始地层厚度恢复方法 | |
US11598893B2 (en) | Seismic rock property prediction in forward time based on 4D seismic analysis | |
CN107945271B (zh) | 基于地质块体追踪的三维压力场建模方法 | |
McKean et al. | Geomechanical and Fracture Network Interpretation of a Devonian Outcrop | |
Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
Le et al. | A machine learning approach for calibrating seismic interval velocity in 3D velocity model | |
Mann et al. | Pluto 4D Monitor Two–a second outstanding success | |
Thomason et al. | Model Creation Using Stacked Surfaces | |
CN117130050A (zh) | 一种用于构建三维地质体的方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |