CN115629420B - 正向地貌构造级次及幅度确定方法、装置、介质及设备 - Google Patents
正向地貌构造级次及幅度确定方法、装置、介质及设备 Download PDFInfo
- Publication number
- CN115629420B CN115629420B CN202211569568.6A CN202211569568A CN115629420B CN 115629420 B CN115629420 B CN 115629420B CN 202211569568 A CN202211569568 A CN 202211569568A CN 115629420 B CN115629420 B CN 115629420B
- Authority
- CN
- China
- Prior art keywords
- landform
- amplitude
- curve
- level
- function
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Abstract
Description
技术领域
本发明涉及石油开采技术领域,具体是关于一种正向地貌构造级次及幅度确定方法、装置、介质及设备。
背景技术
地貌形态对于油气的分布具有控制作用,如对于砂岩油气藏,盆地内的地貌形态是控制沉积砂体的重要因素,决定着储集砂体的发育位置及分布规模。研究表明,地貌形态是制约盆地沉积充填演化的重要因素,宏观地貌决定了盆地的隆坳格局,进而控制着沉积体系的发育演化,微观地貌影响着盆地内局部的水动力条件,控制着局部地区沉积物分散体系的展布规律。对于基岩油气藏,尤其是风化壳型潜山油气藏,裂缝发育带是寻找优质储层的关键因素,研究表明裂缝发育带与风化壳顶面古地貌密切相关,潜山风化壳顶面主要表现为“沟脊”地貌特征,在山脊处,由于位于构造高部位,风化溶蚀作用强烈,易于形成风化溶蚀孔、缝,而山沟谷处,由于地势低,风化溶蚀作用相对较弱,同时山脊处风化剥蚀的产物由于重力作用会堆积于此,进一步减弱了山沟位置的风化作用,因此裂缝发育带主要分布在山脊处,裂缝分布与地貌特征密切相关。
当前,针对地貌表征主要采用三维地震资料进行地貌形态解释,然而受地震资料分辨率制约,地貌形态解释多解性强,而且解释精度与解释人员经验密切相关,另外,对于地貌级次,也未形成统一划分标准。
发明内容
针对上述问题,本发明的目的是提供一种正向地貌构造级次及幅度确定方法、装置、介质及设备,实现了对地貌形态的精细表征。
为实现上述目的,本发明采取以下技术方案:
本发明所述的正向地貌构造级次及幅度确定方法,包括以下步骤:
(1)获取地貌形态曲面函数;
(4)将沟脊幅度曲线作为新的输入曲线重复步骤(1)至(3),直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
所述的正向地貌构造级次及幅度确定方法,优选地,步骤(1)中所述获取地貌形态曲面函数包括如下步骤:
a.根据研究任务,确定目的层层位;
b.在目的层层位的基础上,采用测井声波曲线及密度曲线确定反射系数,结合子波提取制作合成地震记录;
c.在地震记录的基础上,综合利用多种测井曲线,结合VSP测井资料,在过井地震剖面上精细标定目的层;
d.在过井地震剖面的目的层的基础上,综合考虑地震剖面上同相轴的振幅、能量、连续性、相位、频率及波组特征以及地层时间切片,开展地貌形态精细解释工作,最终获取地貌形态曲面;
e.将地貌形态曲面进行不同程度的非线性多项式拟合,形成地貌形态曲面函数。
所述的正向地貌构造级次及幅度确定方法,优选地,步骤(2)中对正向地貌构造级次进行划分,得到第i个平面的级次划分函数曲线,确定级次划分函数曲线的极大值采样点和极小值采样点,具体是通过以下方式进行实现的:
当b>1时,采样点分为极大值采样点和极小值采样点,进行多项式拟合得到极大值拟合函数(x)和极小值拟合函数(x),将(x)的极大值或极小值点与(x)的极大值或极小值点的平方差再开方的中点作为采样点,级次划分函数曲线由进行多项式拟合得到;
其中,b为正向地貌构造形态的划分级次,x为级次划分函数曲线的采样点,w为各采样点的权重系数;n为序号。
所述的正向地貌构造级次及幅度确定方法,优选地,步骤(3)中确定正向地貌构造的沟脊幅度曲线和幅度包括如下步骤:
a.由曲线最大极值点利用多项式拟合生成山脊连线,由曲线最小极值点利用多项式拟合生成山沟连线;山脊连线和山沟连线形成一级古地貌的两条包络线;
b.将两条包络线作为待拟合曲线,利用最小二乘法形成一条含有山脊和山沟幅度特征的幅度曲线;
c.通过幅度曲线的二阶导数求出极大值点和极小值点,求出极大值点到两个相邻极小值点连线的距离为该级别正向地貌构造的幅度H。
本发明提供的正向地貌构造级次及幅度确定装置,包括:
第一处理单元,用于获取地貌形态曲面函数;
第三处理单元,用于确定正向地貌构造的沟脊幅度曲线和幅度:级次划分函数曲线的极大值采样点到相邻两个极小值采样点连线的距离,定义为正向地貌构造幅度H;
第四处理单元,用于将沟脊幅度曲线作为新的输入曲线重复(1)、(2)、(3)操作,直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
本发明提供的一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述的正向地貌构造级次及幅度确定方法步骤。
本发明提供的一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述正向地貌构造级次及幅度确定方法步骤。
本发明由于采取以上技术方案,其具有以下优点:
(1)本发明首次对正向地貌构造形态开展了定量统计,有效解决了地貌形态描述多解性强的难题,方法具有创新性和实用性。
(2)本发明首次实现了对正向地貌构造进行分级分类的定量判别方法。
(3)本发明属于定量判别方法,与传统的定性估算方法相比更具科学性、合理性。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。在整个附图中,用相同的附图标记表示相同的部件。在附图中:
图1是某油田地貌形态解释剖面图;
图2是正向地貌构造级次及幅度确定流程图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
本发明提供一种正向地貌构造级次及幅度确定方法,首次对正向地貌构造形态开展了定量统计,有效解决了地貌形态描述多解性强的难题,且首次实现了对正向地貌构造进行分级分类的定量判别方法,与传统的定性估算方法相比更具科学性和合理性。
本发明提供的正向地貌构造级次及幅度确定方法,包括以下步骤:
(1)获取地貌形态曲面函数;
(4)将沟脊幅度曲线作为新的输入曲线重复(1)、(2)、(3)操作,直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
在上述实施例中,优选地,步骤(1)中所述获取地貌形态曲面函数包括如下步骤:
a.根据研究任务,确定目的层层位;
b.在目的层层位的基础上,采用测井声波曲线及密度曲线确定反射系数,结合子波提取制作合成地震记录;
c.在地震记录的基础上,综合利用多种测井曲线,结合VSP测井资料,在过井地震剖面上精细标定目的层;
d.在过井地震剖面的目的层的基础上,综合考虑地震剖面上同相轴的振幅、能量、连续性、相位、频率及波组特征以及地层时间切片,开展地貌形态精细解释工作,最终获取地貌形态曲面;
e.将地貌形态曲面进行不同程度的非线性多项式拟合,形成地貌形态曲面函数。
当b>1时,采样点分为极大值采样点和极小值采样点,进行多项式拟合得到极大值拟合函数(x)和极小值拟合函数(x),将(x)的极大值或极小值点与(x)的极大值或极小值点的平方差再开方的中点作为采样点,级次划分函数由进行多项式拟合得到;
其中,b为正向地貌构造形态的划分级次,x为级次划分函数曲线的采样点,w为各采样点的权重系数;n为序号。
在上述实施例中,优选地,步骤(3)中确定正向地貌构造的沟脊幅度曲线和幅度包括如下步骤:
a.由曲线最大极值点利用多项式拟合生成山脊连线,由曲线最小极值点利用多项式拟合生成山沟连线;山脊连线和山沟连线形成一级古地貌的两条包络线;
b.将两条包络线作为待拟合曲线,利用最小二乘法形成一条含有山脊和山沟幅度特征的幅度曲线;
c.通过幅度曲线的二阶导数求出极大值点和极小值点,求出极大值点到两个相邻极小值点连线的距离为该级别正向地貌构造幅度h。
本发明还提供一种正向地貌构造级次及幅度确定装置,包括:
第一处理单元,用于获取地貌形态曲面函数;
第三处理单元,用于确定正向地貌构造的沟脊幅度曲线和幅度:级次划分函数曲线的极大值采样点到相邻两个极小值采样点连线的距离,定义为正向地貌构造幅度H;
第四处理单元,用于将沟脊幅度曲线作为新的输入曲线重复(1)、(2)、(3)操作,直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
本发明还提供一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述的正向地貌构造级次及幅度确定方法步骤。
本发明还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述正向地貌构造级次及幅度确定方法步骤。
实施例1:
1)应用LANDMARK商业软件进行潜山油田A地貌形态解释(附图1),得到地貌形态曲面;
6)将G(x)和S(x)的级值平方中点进行多项式拟合,得到2级曲线(如附图2.3);
7)将2级曲线作为采样集合重复步骤4到7,得到2级和3级正向地貌幅度(如附图2.4、2.5),至此第i面的所有级次划分完成;
h={45,70,55}
h'={82}
另外,需要说明的是,正向地貌指的是地壳上相对隆起的构造单元。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种正向地貌构造级次及幅度确定方法,其特征在于,包括以下步骤:
(1)获取地貌形态曲面函数;
(2)基于地貌形态曲面函数,对正向地貌构造级次进行划分,得到第i个平面的级次划分函数曲线,集成每个平面的,形成正向地貌构造级次曲面;确定级次划分函数曲线的极大值采样点和极小值采样点;具体地,步骤(2)是通过以下方式进行实现的:
当b>1时,采样点分为极大值采样点和极小值采样点,进行多项式拟合得到极大值拟合函数(x)和极小值拟合函数(x),将(x)的极大值或极小值点与(x)的极大值或极小值点的平方差再开方的中点作为采样点,级次划分函数曲线由进行多项式拟合得到;
其中,b为正向地貌构造形态的划分级次,x为级次划分函数曲线的采样点,w为各采样点的权重系数;
具体地,步骤(3)中确定正向地貌构造的沟脊幅度曲线和幅度包括如下步骤:
a.由曲线最大极值点利用多项式拟合生成山脊连线,由曲线最小极值点利用多项式拟合生成山沟连线;山脊连线和山沟连线形成一级古地貌的两条包络线;
b.将两条包络线作为待拟合曲线,利用最小二乘法形成一条含有山脊和山沟幅度特征的幅度曲线;
c.通过幅度曲线的二阶导数求出极大值点和极小值点,求出极大值点到两个相邻极小值点连线的距离为正向地貌构造的幅度H;
(4)将沟脊幅度曲线作为新的输入曲线重复步骤(1)至(3),直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
2.根据权利要求1所述的正向地貌构造级次及幅度确定方法,其特征在于,步骤(1)中所述获取地貌形态曲面函数包括如下步骤:
a.根据研究任务,确定目的层层位;
b.在目的层层位的基础上,采用测井声波曲线及密度曲线确定反射系数,结合子波提取制作合成地震记录;
c.在地震记录的基础上,综合利用多种测井曲线,结合VSP测井资料,在过井地震剖面上精细标定目的层,所述VSP是指垂直地震剖面测井;
d.在过井地震剖面的目的层的基础上,综合考虑地震剖面上同相轴的振幅、能量、连续性、相位、频率及波组特征以及地层时间切片,开展地貌形态精细解释工作,最终获取地貌形态曲面;
e.将地貌形态曲面进行不同程度的非线性多项式拟合,形成地貌形态曲面函数。
3.一种基于权利要求1或2的正向地貌构造级次及幅度确定方法的正向地貌构造级次及幅度确定装置,其特征在于,包括:
第一处理单元,用于获取地貌形态曲面函数;
第三处理单元,用于确定正向地貌构造的沟脊幅度曲线和幅度:级次划分函数曲线的极大值采样点到相邻两个极小值采样点连线的距离,定义为正向地貌构造幅度H;
第四处理单元,用于将沟脊幅度曲线作为新的输入曲线重复(1)、(2)、(3)操作,直至n级沟脊,当n级沟脊只有3个拐点时,完成地貌所有级次沟脊识别并确定构造级次。
4.一种计算机存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1或2所述的正向地貌构造级次及幅度确定方法步骤。
5.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1或2所述正向地貌构造级次及幅度确定方法步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211569568.6A CN115629420B (zh) | 2022-12-08 | 2022-12-08 | 正向地貌构造级次及幅度确定方法、装置、介质及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211569568.6A CN115629420B (zh) | 2022-12-08 | 2022-12-08 | 正向地貌构造级次及幅度确定方法、装置、介质及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115629420A CN115629420A (zh) | 2023-01-20 |
CN115629420B true CN115629420B (zh) | 2023-03-17 |
Family
ID=84911206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211569568.6A Active CN115629420B (zh) | 2022-12-08 | 2022-12-08 | 正向地貌构造级次及幅度确定方法、装置、介质及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115629420B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956238A (zh) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | 一种古地貌的恢复方法及装置 |
AU2020102025A4 (en) * | 2020-07-02 | 2020-10-01 | China University Of Geosciences, Beijing | Method for Reconstructing Basin Paleogeomorphology |
CN114265116A (zh) * | 2022-03-02 | 2022-04-01 | 中海油研究总院有限责任公司 | 花岗岩类潜山风化壳沟脊幅度定量统计方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077558B (zh) * | 2012-05-14 | 2015-12-09 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 |
CN109886457A (zh) * | 2019-01-15 | 2019-06-14 | 中国石油大学(华东) | 一种判断复杂岩性古潜山有利风化壳储层发育范围的方法 |
CN111539080A (zh) * | 2020-04-23 | 2020-08-14 | 南京英尼格玛工业自动化技术有限公司 | 一种基于动态网格划分的非标零件载荷分析方法及系统 |
CN115272866A (zh) * | 2022-08-24 | 2022-11-01 | 中国质量认证中心 | 一种基于正负地貌的山区地貌形态分类的方法与应用 |
-
2022
- 2022-12-08 CN CN202211569568.6A patent/CN115629420B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956238A (zh) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | 一种古地貌的恢复方法及装置 |
AU2020102025A4 (en) * | 2020-07-02 | 2020-10-01 | China University Of Geosciences, Beijing | Method for Reconstructing Basin Paleogeomorphology |
CN114265116A (zh) * | 2022-03-02 | 2022-04-01 | 中海油研究总院有限责任公司 | 花岗岩类潜山风化壳沟脊幅度定量统计方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN115629420A (zh) | 2023-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hengl et al. | Mathematical and digital models of the land surface | |
Chaplot et al. | Accuracy of interpolation techniques for the derivation of digital elevation models in relation to landform types and data density | |
Fraisse et al. | Delineation of site-specific management zones by unsupervised classification of topographic attributes and soil electrical conductivity | |
Hamilton et al. | Acoustic seabed segmentation from direct statistical clustering of entire multibeam sonar backscatter curves | |
CN106599129A (zh) | 一种顾及地形特点的多波束点云数据去噪方法 | |
Yang et al. | Filtering of airborne LiDAR bathymetry based on bidirectional cloth simulation | |
Fagbohun et al. | GIS-based estimation of soil erosion rates and identification of critical areas in Anambra sub-basin, Nigeria | |
Ramillien et al. | Global bathymetry derived from altimeter data of the ERS-1 geodetic mission | |
CN110335352B (zh) | 一种机载激光雷达点云的双基元多分辨率层次滤波方法 | |
CN109410265B (zh) | 一种基于往期dem辅助的tin滤波改进算法 | |
CN111340723B (zh) | 一种地形自适应的机载LiDAR点云正则化薄板样条插值滤波方法 | |
CN105868482A (zh) | 一种沉积期微幅度古隆起的推算方法及装置 | |
Dromart et al. | Deposition and erosion of a Light-Toned Yardang-forming unit of Mt Sharp, Gale crater, Mars | |
CN109490978B (zh) | 一种起伏地层的频率域快速高精度正演方法 | |
CN115629420B (zh) | 正向地貌构造级次及幅度确定方法、装置、介质及设备 | |
CN114862715A (zh) | 一种融合地形特征语义信息的tin渐进加密去噪方法 | |
Skentos et al. | Landform analysis using terrain attributes. A GIS application on the Island of Ikaria (Aegean Sea, Greece) | |
Betz-Nutz et al. | Development of the morphodynamics on Little Ice Age lateral moraines in 10 glacier forefields of the Eastern Alps since the 1950s | |
Chakraborty et al. | Acoustic seafloor sediment classification using self-organizing feature maps | |
Kirby et al. | Anomalously large gradients in version 1 of the “GEODATA 9 SECOND” Digital Elevation Model of Australia, and their effects on gravimetric terrain corrections | |
Deng et al. | Terrain analysis | |
Chen et al. | Computational morphology of debris and alluvial fans on irregular terrain using the visibility polygon | |
Schneider et al. | Reconstruction of a palaeosurface and archaeological site location in an anthropogenic drift sand area | |
Yin et al. | A technical note on the application of interpolation method on modeling of 3D surface rock layer | |
Van Nieuwenhuizen | An Analysis of Preprocessing Techniques for the Removal of Transportation Embankments and Surface Roughness in Fine-Resolution LiDAR DEMs |
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 |