CN100363755C - 绘制含正断层地质体构造等值线图的矩形网网格化方法 - Google Patents
绘制含正断层地质体构造等值线图的矩形网网格化方法 Download PDFInfo
- Publication number
- CN100363755C CN100363755C CNB2005100661096A CN200510066109A CN100363755C CN 100363755 C CN100363755 C CN 100363755C CN B2005100661096 A CNB2005100661096 A CN B2005100661096A CN 200510066109 A CN200510066109 A CN 200510066109A CN 100363755 C CN100363755 C CN 100363755C
- Authority
- CN
- China
- Prior art keywords
- grid
- data
- node
- tomography
- interpolation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000010422 painting Methods 0.000 title claims abstract description 9
- 238000003325 tomography Methods 0.000 claims description 57
- 238000005516 engineering process Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 238000003491 array Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 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 description 2
- 238000009826 distribution Methods 0.000 description 2
- 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 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
- Image Generation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种绘制含正断层地质体构造等值线时的矩形网网格化方法，散点数据预处理，建立索引；断层数据预处理，建立索引；对每个网格点进行数据搜索并插值，不必在全部离散点数组中查找数据点，从而大幅提高数据的搜索速度。
Description
技术领域
本发明涉及地球物理勘探成图技术，尤其是实现自动追踪绘制含正断层地质体构造等值线图的矩形网网格化方法。
技术背景
在石油地质勘探工作时，需要将采集并处理后的地震数据绘制成地质构造图，对于带有复杂断层组合线(多边形)的数据，要先进行网格化才能自动追踪并绘制出等值线。网格化是对于规则的矩形网格上的节点利用周围的离散数据点采用拟合或者插值的方法计算出网格节点的值。
网格化技术是自动绘制等值线图技术的关键，由于地震数据多数情况下含有断层，使得地学曲面非常复杂，由于原始数据常常分布不均，有时呈现出一簇簇的极端情况，这就会造成网格化时某些方向权重过大而另一些方向权重偏小，造成估值不准，影响网格化的合理性，从而使自动追踪绘制的效果不佳，现有的技术不能实现对包含断层组合线的数据进行网格化。
发明内容
本发明所要解决的问题是提供一种插值准确、能够实现自动追踪绘制含正断层地质体构造等值线图的矩形网网格化方法。
本发明提供以下技术方案：绘制含正断层地质体构造等值线图的矩形网网格化方法，首先采集并处理地震数据，还包括以下步骤：
1)散点数据预处理，建立索引；
在用户指定的勘探区域上建立插值网格，在比插值网格稍微大一点的区域上建立索引网格；
将离散数据投影到索引网格的具体单元上并记录散点的索引号；
2)断层组合线数据预处理，建立索引；
对断层组合线在x、y平面投影并且建立索引，记录每个网格单元中穿过单元的断层组合线段数，记录穿过的断层号、断层线的段号；
对网格节点与数据点组成的线段，经过的索引网格单元记录的断层组合线进行相交判断；
3)对每个网格点进行数据搜索并插值；
逐列批量、滚动地判断插值网格的所有节点，落在断层组合线多边形内的节点直接赋予无效值；
对每一个落在断层多边形之外的网格节点，以网格节点为圆心、给定的半径做搜索圆，将搜索圆均分成八个扇区，落在搜索圆内的数据分配到相应的扇区当中，从当前网格节点出发，由内向外、由近及远旋转搜索数据，用落在圆内的与该节点在断层同一侧的层位数据拟合局部曲面或直接插值计算每一个网格节点的值；
4)将网格数据和断裂线数据存入文件，供等值线追踪并绘制含正断层地质体构造等值线图。
本发明还提供以下技术方案：索引网格的密度以每个网格单元包含不超过4到6个数据点为准。
判断插值网格的所有节点时，当在网格的一列中有连续的节点所在的单元有断层组合线穿过，逐个判断节点是否在断层组合线所属的断层多边形内；
当连续的节点所在的单元没有断层段穿过，而在下一个和上一个节点所在的单元有断层段穿过，则：
当下一个和上一个节点所在的单元的断层段所属断层多边形是同一个，则中间的一个节点如果在这个多边形内，那么这些连续的节点就都在此断层多边形内，否则就都不在此断层多边形内；
当下一个和上一个节点所在的单元的断层段所属断层多边形没有相同的，则这些连续的节点就都不在此断层多边形内。
采用沿垂向网格线滚动判断技术对断层多边形数据进行预处理。
对网格节点进行插值的算法包括角度校正的反比距离加权插值法、加权最小二乘曲面拟合、薄板样条曲面插值、克里金插值法。
对于整个绘图区数据点很少的情况，采用滑动加权趋势面误差校正网格化算法。
本发明对离散点数据进行预处理，建立索引网格，不必在全部离散点数组中查找数据点，从而大幅提高数据的搜索速度。本发明矩形网网格化所插值准确、合理，自动追踪绘制含正断层地质体构造等值线图效果好。
附图说明
图1是“簇”状数据集及八象限搜索示意图；
图2是网格节点落在断层多边形内外的快速判断示意图；
图3是数据点(图中井点)较少情况时的网格化效果示意图。
具体实施方式
本发明实施例采用以下步骤：
第一步，对离散点数据进行预处理。根据网格的范围确定一个稍大一点的范围，根据这个范围建立索引网格并将离散数据索引到该网格中。
对离散点数据进行预处理，即在用户指定的区域上建立插值网格，在比插值网格稍微大一点的区域上建立索引网格，索引网格的密度以每个网格单元包含不超过4到6个数据点为准。将离散数据投影到索引网格的具体单元上并记录散点的索引号，根据插值网格节点所在的索引网格单元螺旋状由内向外逐个单元寻找索引网格单元内的离散点。这样不必在全部离散点数组中查找数据点，从而大幅提高数据的搜索速度。
第二步，断层多边形预处理：采用沿垂向网格线滚动判断技术，记录穿过网格单元的断层号、断裂线段号，对网格的每一列采用滚动判断技术快速判断那些网格节点落在断层多边形内，从而直接赋予无效值跳过计算。
本发明所述的绘制含正断层地质体构造等值线时的矩形网网格化方法，对断层组合线数据预处理。在对插值网格节点进行插值时需要对搜索到的数据点进行是否与网格节点在断层同一侧的判断，具体说来就是进行线段相交判断，即网格节点与数据点组成的线段是否与任意一条断层组合线相交。这个过程计算量很大，为了提高运算速度，这里对断层组合线数据进行预处理，即对断层组合线在x、y平面投影并且建立索引：在每个网格单元中记录如下信息：穿过该单元的断层组合线段数；如果穿过该单元的断层组合线段数不为0，还要记录穿过的断层号、断层线的段号。进行相交判断时，只需要对网格节点与数据点组成的线段经过的索引网格单元记录的断层组合线进行少量的相交判断运算，从而提高网格化运算速度，节省计算时间。
在图2中，空心或实心圆表示的连续多个网格节点所在的网格单元中都没有断层多边形穿过，但它们的情况并不一样：实心点所在单元的上下单元穿过的断层多边形是不同的断层多边形，因此这些节点都落在断层多边形外；空心点所在单元的上下单元穿过的断层多边形是同一个断层多边形，此时又分为两种情况：如果这一串空心点其中有一个落在断层多边形内，则它们就都落在该断层多边形内(图2中空心小圆圈所示的网格节点)；否则，就都落在断层多边形之外(图2中同心小圆圈所示的网格节点)。在有断层多边形穿过的网格单元，其上的网格节点是否落在断层多边形内需要逐点判断，但这样的情况已经大为减少，并且知道单元中穿过的断层号，因此断层预处理速度非常快。
本发明在断层组合线数据建立索引的基础上，还对插值网格的所有节点逐列批量、滚动地判断出哪些网格节点落在断层组合线多边形内从而直接赋予无效值而不必计算。方法如下：分成两种情况：第一，在网格的一列中，有一些连续的节点所在的单元都有断层组合线穿过，此时对它们逐个判断是否在断层组合线所属的断层多边形内；第二，有一些连续的节点所在的单元没有断层段穿过，而在它们的下一个和上一个节点所在的单元有断层段穿过，此时又分为两种情况：其一，下一个和上一个节点所在的单元的断层段所属断层多边形是同一个，则中间的一个节点如果在这个多边形内，那么这些连续的节点就都在此断层多边形内，否则就都不在此断层多边形内；其二，下一个和上一个节点所在的单元的断层段所属断层多边形没有相同的，则这些连续的节点就都不在此断层多边形内。该方法可以大大提高运算速度。
第三步，对每一个落在断层多边形之外的网格节点，以该点为圆心、以用户指定的半径作搜索圆，采用八象限分块选取技术搜索散点数据，以待计算网格节点为圆心，以给定的半径做搜索圆，将这个搜索圆均分成八个扇区，落在搜索圆内的数据分配到相应的扇区当中，在每个扇区中只取距离圆心最近的若干个数据点，较好地解决了数据分布不均的问题。
本发明用落在圆内的与该节点在断层同一侧的层位数据拟合局部曲面或直接插值计算每一个网格节点的值。
采用数据螺旋搜索法搜索各象限中的数据，从当前网格节点出发，由内向外由近及远旋转搜索数据，直到所有扇区中的数据点满足要求时自动停止。
本发明由当前网格点出发螺旋向外寻找经过索引的数据点，当所有八个扇区中数据点数满足条件时即自动停止数据搜索。
本发明断层组合线两侧的数据处理原则是：断层两盘的数据特性是不连续的，因此断层两盘的数据不能同时参与运算，即要计算一个网格节点的值，只有与其在同一盘上的数据才能参与运算。
本发明对于整个绘图区数据点很少的情况(如井数据)，采用了加权趋势面加误差校正的网格化方法。先用加权最小二乘一次多项式拟合计算出粗网格(每个节点计算一个方程)。这时的粗网格在有数据的地方是比较精确的，在数据区以外则给出了趋势，对于发生在数据区以外外推的z值超出原始z值范围很多的情况，首先将原始数据z值范围适当外延，得到一个较合理的z值范围，然后将粗网格压缩到这个数值范围，此时粗网格在有数据点的地方不够准确，但它的趋势是合理的。用这个粗网格内插出细网格，然后求出原始数据点处的误差，然后将这些误差分配到细网格上，这时的细网格既保证了在有数据点处的准确性，又有正确的趋势，因此是理想的。
第四步，用搜索到的离散数据，采用插角度校正的反比距离加权插值法或加权最小二乘曲面拟合或薄板样条曲面插值或克里金插值法给当前网格节点插值。
本发明网格化的插值/拟合算法共有四种：角度校正的反比距离加权法、加权最小二乘曲面拟合、简单克里金估计和薄板样条曲面插值法。
本发明关于距离加权插值法，对于如何估计各数据点相对于待估点的权值，基于：距离越小权越大(距离某次方的倒数)和内插比外推权大的原则，在应用权值计算法时同时考虑距离和位置关系，为了求一个数据点的权值，用其余数据点和它本身及待估点形成的夹角余弦被1减对
作影响因子，其越接近内插则夹角余弦越小。
cosang＝1-(((xco-x[m])*(xco-x[n])+(yco-y[m])*
(yco-y[n]))/(dis[m]*dis[n]))
当全部数据点对当前这个数据点的影响值direct计算出来后，每个数据点对待估点的权值就是：
第五步，将网格化后的数据和断裂线数据存入文件中，采用通常的方法对数据进行等值线追踪绘图。
Claims (4)
1.一种绘制含正断层地质体构造等值线图的矩形网网格化方法，首先采集并处理地震数据，其特征还包括以下步骤：
1)散点数据预处理，建立索引；
在用户指定的勘探区域上建立插值网格，在比插值网格稍微大一点的区域上建立索引网格；
将离散数据投影到索引网格的具体单元上并记录散点的索引号；
2)断层组合线数据预处理，建立索引；
对断层组合线在x、y平面投影并且建立索引，记录每个网格单元中穿过单元的断层组合线段数，记录穿过的断层号、断层线的段号；
对网格节点与数据点组成的线段，经过的索引网格单元记录的断层组合线进行相交判断；
3)对每个网格点进行数据搜索并插值；
逐列批量、滚动地判断插值网格的所有节点，落在断层组合线多边形内的节点直接赋予无效值；
对每一个落在断层多边形之外的网格节点，以网格节点为圆心、给定的半径做搜索圆，将搜索圆均分成八个扇区，落在搜索圆内的数据分配到相应的扇区当中，从当前网格节点出发，由内向外、由近及远旋转搜索数据，用落在圆内的与该节点在断层同一侧的层位数据拟合局部曲面或直接插值计算每一个网格节点的值；
4)将网格数据和断裂线数据存入文件，供等值线追踪使用。
索引网格的密度以每个网格单元包含不超过4到6个数据点为准。
判断插值网格的所有节点时，当在网格的一列中有连续的节点所在的单元有断层组合线穿过，逐个判断节点是否在断层组合线所属的断层多边形内；
当连续的节点所在的单元没有断层段穿过，而在下一个和上一个节点所在的单元有断层段穿过，则：
当下一个和上一个节点所在的单元的断层段所属断层多边形是同一个，则中间的一个节点如果在这个多边形内，那么这些连续的节点就都在此断层多边形内，否则就都不在此断层多边形内；
当下一个和上一个节点所在的单元的断层段所属断层多边形没有相同的，则这些连续的节点就都不在此断层多边形内。
2.根据权利要求1所述的绘制含正断层地质体构造等值线时的矩形网网格化方法，其特征在于：采用沿垂向网格线滚动判断技术对断层多边形数据进行预处理。
3.根据权利要求1所述的绘制含正断层地质体构造等值线时的矩形网网格化方法，其特征在于：对网格节点进行插值的算法包括角度校正的反比距离加权插值法、加权最小二乘曲面拟合、薄板样条曲面插值、克里金插值法。
4.根据权利要求1所述的绘制含正断层地质体构造等值线时的矩形网网格化方法，其特征在于：对于整个绘图区数据点很少的情况，采用滑动加权趋势面误差校正网格化算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005100661096A CN100363755C (zh) | 2005-04-21 | 2005-04-21 | 绘制含正断层地质体构造等值线图的矩形网网格化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005100661096A CN100363755C (zh) | 2005-04-21 | 2005-04-21 | 绘制含正断层地质体构造等值线图的矩形网网格化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1854762A CN1854762A (zh) | 2006-11-01 |
| CN100363755C true CN100363755C (zh) | 2008-01-23 |
Family
ID=37195060
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2005100661096A Expired - Fee Related CN100363755C (zh) | 2005-04-21 | 2005-04-21 | 绘制含正断层地质体构造等值线图的矩形网网格化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100363755C (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101930080B (zh) * | 2009-06-26 | 2012-08-15 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101702237B (zh) * | 2009-11-10 | 2011-11-23 | 南京大学 | 一种基于二次误差测度的层位面拟合方法 |
| CN102542136B (zh) * | 2010-12-17 | 2014-09-10 | 中国地质大学（北京） | 计算机离散数据网格化参数提取方法及操作步骤 |
| CN102169025B (zh) * | 2010-12-30 | 2013-03-27 | 中冶连铸技术工程股份有限公司 | 高速温度云图的绘制方法 |
| CN102759755B (zh) * | 2011-04-29 | 2016-08-03 | 黄桂芝 | 采用旋转tin网和非剖面法直接制作平、立面图的地质勘查方法 |
| CN102385067B (zh) * | 2011-08-03 | 2013-09-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 含逆断层等值线绘图方法 |
| CN103425856A (zh) * | 2012-05-24 | 2013-12-04 | 泰瑞数创科技(北京)有限公司 | 一种依据河道测量断面数据自动生成河道地形的方法 |
| CN103424773B (zh) * | 2012-05-25 | 2016-02-03 | 中国石油化工股份有限公司 | 一种基于层位拉平法的古地貌恢复方法 |
| CN102982566B (zh) * | 2012-06-19 | 2016-03-30 | 克拉玛依红有软件有限责任公司 | 一种基于最小曲率法插值的含断层的等值线图自动生成方法 |
| CN104166158B (zh) * | 2013-05-17 | 2016-11-09 | 中国石油天然气集团公司 | 一种地震数据非均匀间距等值线交互追踪方法 |
| CN103473473B (zh) * | 2013-09-26 | 2018-03-02 | 深圳市华傲数据技术有限公司 | 一种基于散点图的数据质量检测方法及系统 |
| CN104635262B (zh) * | 2013-11-13 | 2017-06-06 | 中国石油天然气集团公司 | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 |
| CN107991702B (zh) * | 2016-10-26 | 2019-05-17 | 中国石油化工股份有限公司 | 正断层时序的简单精准判定方法 |
| CN107146267B (zh) * | 2017-05-12 | 2020-08-11 | 中国石油天然气股份有限公司 | 一种绘制等值线图的方法以及系统 |
| CN107688205A (zh) * | 2017-07-20 | 2018-02-13 | 中国石油天然气集团公司 | 相对地质构造图的建立方法和装置 |
| CN109471161B (zh) * | 2018-09-26 | 2020-07-10 | 中国石油天然气集团有限公司 | 地质构造属性剖面的建立方法和装置 |
| CN109751044A (zh) * | 2019-01-29 | 2019-05-14 | 河南理工大学 | 一种油区热储层顶面温度等值线图编制方法 |
| CN114114405A (zh) * | 2020-08-26 | 2022-03-01 | 中国石油天然气股份有限公司 | 基于离散数据网格化的地层真厚度确定方法及装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2165674A (en) * | 1984-10-16 | 1986-04-16 | Conoco Inc | Plotting of hydrocarbon indication using seismic and heat flow data |
-
2005
- 2005-04-21 CN CNB2005100661096A patent/CN100363755C/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2165674A (en) * | 1984-10-16 | 1986-04-16 | Conoco Inc | Plotting of hydrocarbon indication using seismic and heat flow data |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101930080B (zh) * | 2009-06-26 | 2012-08-15 | 中国石油集团东方地球物理勘探有限责任公司 | 曲面拟合海底电缆二次定位方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1854762A (zh) | 2006-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100363755C (zh) | 绘制含正断层地质体构造等值线图的矩形网网格化方法 | |
| CN103425772B (zh) | 一种具有多维信息的海量数据查询方法 | |
| CN102385067B (zh) | 含逆断层等值线绘图方法 | |
| CN103747419B (zh) | 一种基于信号强度差值与动态线性插值的室内定位方法 | |
| MXPA02003097A (es) | Ubicaciones de pozo optimas determinadas a partir del modelo de deposito de 3d. | |
| CN104597494B (zh) | 地震地层体分析方法及装置 | |
| CN107861132A (zh) | 一种gps轨迹优化方法及装置 | |
| CN102222366B (zh) | 复杂空间曲面的拟合方法 | |
| CN104181599A (zh) | 一种基于近地表层的折射波静校正处理方法以及系统 | |
| CN103699678A (zh) | 一种基于多阶段分层采样的层次聚类方法和系统 | |
| CN105473809A (zh) | 用于评估生产策略规划的方法 | |
| CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
| CN107843920A (zh) | 一种聚类分析优化粒子群储层属性预测方法 | |
| CN109829029A (zh) | 一种顾及居住建筑属性的城市人口空间化方法及系统 | |
| CN112037107B (zh) | 一种自然保护地边界重叠特征的量化方法 | |
| CN113139760A (zh) | 基于风雨大数据的台风风险综合评估方法及系统 | |
| CN103971300A (zh) | 基于gis空间信息的变电站站址定位方法 | |
| CN110135652B (zh) | 一种长期汛期径流预测方法 | |
| CN107831516A (zh) | 融合gnss和地面监测网获取大坝实时高精度位移的方法 | |
| CN110543660B (zh) | 一种低冲击开发模拟方法、系统及相关装置 | |
| CN103559389B (zh) | 一种县域耕地等级变化监测样点布设的方法 | |
| CN113435940B (zh) | 矢量地理空间数据在线服务多因子运营计价系统 | |
| CN103886129B (zh) | 将测井数据离散到储层网格模型的方法和装置 | |
| CN103543478A (zh) | 地质形态插值的km方法 | |
| CN108491482B (zh) | 一种顾及邻近度关系的地质图动态综合方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080123 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |