CN106372382B - 一种基于膨胀算法的岩性插值方法 - Google Patents

一种基于膨胀算法的岩性插值方法 Download PDF

Info

Publication number
CN106372382B
CN106372382B CN201610686362.XA CN201610686362A CN106372382B CN 106372382 B CN106372382 B CN 106372382B CN 201610686362 A CN201610686362 A CN 201610686362A CN 106372382 B CN106372382 B CN 106372382B
Authority
CN
China
Prior art keywords
lithology
grid
shot point
expansion
filled
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
Application number
CN201610686362.XA
Other languages
English (en)
Other versions
CN106372382A (zh
Inventor
耿春
李亚林
敬龙江
罗文�
黎书琴
何光明
巫骏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BGP Inc filed Critical BGP Inc
Priority to CN201610686362.XA priority Critical patent/CN106372382B/zh
Publication of CN106372382A publication Critical patent/CN106372382A/zh
Application granted granted Critical
Publication of CN106372382B publication Critical patent/CN106372382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Abstract

本发明提供了一种基于膨胀算法的岩性插值方法。所述岩性插值方法包括以下步骤：对工区进行网格化，确保每个炮点都在网格范围内；用炮点的岩性填充炮点所在的网格；对每个炮点进行膨胀，每个炮点的每一轮膨胀均以与上一轮膨胀被岩性填充的网格为基准点，将与所述基准点相邻且未被填充岩性的网格填充为与炮点相同的岩性，其中，第一轮膨胀的基准点为炮点；每一轮膨胀结束后，判断工区的网格是否均被填充岩性，如果存在未被填充岩性的网格，则继续进行下一轮膨胀，直到整个工区的网格均被填充岩性。本发明利用膨胀的方法来进行岩性的横向插值，解决了岩性插值中的生产难题，为其它类似非数值意义的插值提供了一种新的思路。

Description

一种基于膨胀算法的岩性插值方法
技术领域
本发明涉及地震勘探资料的采集领域，更具体地讲，涉及一种基于膨胀算法的岩性插值方法。
背景技术
插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。传统的插值方法如克里金插值、样条插值等主要适用于数字点的插值。
在地震勘探的近地表结构中，岩性是纵横向变化的，每个激发井的不同深度都有不同的岩性，在平面切片图上，每口井是一些离散的点，岩性用岩性符号(例如，颜色)表示，井与井之间存在的大量空隙，而对于没有具体数字意义的岩性，要想对其插值，目前并没有有效的解决办法。
发明内容
针对现有技术中存在的不足，本发明的目的之一在于解决上述现有技术中存在的一个或多个问题。例如，本发明的目的之一在于针对传统插值算法不能够适用于非数值意义的岩性点，提供一种能够对没有具体数字意义的岩性进行插值的方法。
为了实现上述目的，本发明提供了一种基于膨胀算法的岩性插值方法。所述岩性插值方法包括以下步骤：对工区进行网格化，确保每个炮点都在网格范围内；用炮点的岩性填充炮点所在的网格；对每个炮点进行膨胀，每个炮点的每一轮膨胀均以与上一轮膨胀被岩性填充的网格为基准点，将与所述基准点相邻且未被填充岩性的网格填充为与炮点相同的岩性，其中，第一轮膨胀的基准点为炮点；每一轮膨胀结束后，判断工区的网格是否均被填充岩性，如果存在未被填充岩性的网格，则继续进行下一轮膨胀，直到整个工区的网格均被填充岩性。
根据本发明基于膨胀算法的岩性插值方法，可以按炮点线号和点号的顺序完成所述对每个炮点进行膨胀的步骤。
根据本发明基于膨胀算法的岩性插值方法，在所述用炮点的岩性填充炮点所在的网格的步骤中，如果同一个网格出现多个炮点，则用第一个炮点的岩性填充网格。
根据本发明基于膨胀算法的岩性插值方法，每个炮点的膨胀完成后再进行下一个炮点的膨胀。
根据本发明基于膨胀算法的岩性插值方法，根据工区的大小和计算精度对工区进行网格化。
与现有技术相比，本发明的有益效果包括：采用膨胀算法来对岩性进行插值，本发明采用网格化的方法，能够快速的对岩性进行插值，准确度高，解决了岩性插值中的生产难题。
附图说明
通过下面结合附图进行的描述，本发明的上述和其他目的和特点将会变得更加清楚，其中：
图1是根据本发明示例性实施例的基于膨胀算法的岩性插值方法的流程图。
图2是根据本发明示例性实施例的基于膨胀算法的岩性插值方法对任一岩性点进行膨胀的示意图。
图3是根据本发明示例性实施例的基于膨胀算法的岩性插值方法对岩性点进行网格化的示意图。
图4是根据本发明示例性实施例的基于膨胀算法的岩性插值方法完成第一轮膨胀的示意图。
图5是根据本发明示例性实施例的基于膨胀算法的岩性插值方法完成第二轮膨胀的示意图。
图6是岩性单点显示平面图。
图7是图6所示的岩性点采用根据本发明示例性实施例的基于膨胀算法的岩性插值方法进行岩性插值后的效果图。
附图标记说明：
A、B、C和D-炮点
具体实施方式
在下文中，将结合附图和示例性实施例详细地描述根据本发明的基于膨胀算法的岩性插值方法。
对于没有具体数值意义的岩性进行插值，本申请不是采用传统的插值方法进行数值拟合，而是采用点膨胀的方法。本申请的具体思路为：首先根据计算的精度对工区进行网格化，把所有的岩性点包括在网格范围内，其次把炮点所在的网格用该炮点的岩性进行填充，接下来进行第一轮膨胀，把每个炮点所在网格周边的8个网格填充与中间网格相同的岩性，膨胀完后，检查是否填充满整个网格，如果没有填充完，再继续第下一轮的膨胀，直到填充完整个工区网格为止。
图1是根据本发明示例性实施例的基于膨胀算法的岩性插值方法的流程图。图2是根据本发明示例性实施例的基于膨胀算法的岩性插值方法对任一岩性点进行膨胀的示意图。
如图1所示，根据本发明示例性实施例的基于膨胀算法的岩性插值方法首先根据计算的精度对工区进行网格化，确保所有的炮点(岩性点)都在网格范围内；然后用炮点的岩性填充炮点所在的网格；接着按顺序对每个炮点进行膨胀，并在每一轮膨胀结束后判断工区的网格是否均被填充岩性，如果存在未被填充岩性的网格，则继续进行下一轮膨胀，直到整个工区的网格均被填充岩性。其中，膨胀是将与物体接触的所有背景点合并到该物体中，使边界向外部扩张的过程。在膨胀的过程中，如果出现待扩张(或待填充)的网格已填充岩性，则忽略(或跳过)该网格。
如图2所示，对于任一炮点A其膨胀过程为：第一轮膨胀是将与炮点A所在网格相邻的8个网格(图2中标1的网格)填充为与炮点A所在网格相同的岩性；第二轮膨胀是将在第一轮膨胀的基础上，以第一轮膨胀被岩性填充的网格(即图2中标1的网格)为基准点，将与基准点相邻的网格(即图2中标2的网格)填充为与炮点A相同的岩性；以此类推，进行第三轮膨胀，……，第N轮膨胀。
按炮点线号与点号的顺序完成对每个炮点进行膨胀，其中，炮点的线号与点号是炮点的编号，用于区别炮点身份的参数，例如，炮点线号代表横向上的编号，点号代表纵向上的编号(需要说明的是线号不一定是横向、点号不一定是纵向，这里只是举例说明而已)。每个炮点的膨胀完成后再进行下一个炮点的膨胀。如果同一个网格出现多个炮点，则用第一个炮点的岩性填充网格。
为了更好地理解本发明的上述示例性实施例，下面结合具体示例和附图对其进行进一步说明。
示例
图3是根据本发明示例性实施例的基于膨胀算法的岩性插值方法对岩性点进行网格化的示意图。图4是根据本发明示例性实施例的基于膨胀算法的岩性插值方法完成第一轮膨胀的示意图。图5是根据本发明示例性实施例的基于膨胀算法的岩性插值方法完成第二轮膨胀的示意图。
某工区有四个炮点A、B、C和D，以该工区为例说明本申请基于膨胀算法的岩性插值方法。
第一步，网格化(如图3所示)：计算观测系统的方位角，按照方位角把整个观测系统进行旋转，切换到水平模式。再根据最大最小的坐标范围与网格的大小(计算精度)对整个观测系统进行网格化，确保每个炮点都在网格范围内。
第二步：膨胀扩张：按照炮点线号、点号的顺序(例如图4和5中的A→B→C→D)，把炮点所在的网格用该炮点的岩性进行填充(同一个网格出现多个炮点，则填充第一个)，接下来按照顺序进行第一轮膨胀，把每个炮点所在网格周边的8个网格填充与中间网格相同的岩性(图4所示)，如果出现待扩张的网格已填充岩性，则忽略该网格。
第三步：第一轮膨胀完后，检查整个工区网格是否被岩性填充满，如果没有填充满，再继续第二轮的膨胀(图5所示)，以此类推，直到填充完整个工区网格为止。
第四步：获得工区岩性插值结果图。图6是岩性单点显示平面图。图7是图6所示的岩性点进行岩性插值后的效果图。如图6和图7所示，填充完整个网格后，获得岩性插值的结果图，可以输出显示结果。
综上所述，针对传统插值算法不能够适用于非数值意义的岩性点，本发明利用膨胀的方法来进行岩性的横向插值，解决了岩性插值中的生产难题，为其它类似非数值意义的插值提供了一种新的思路。
尽管上面已经通过结合示例性实施例描述了本发明，但是本领域技术人员应该清楚，在不脱离权利要求所限定的精神和范围的情况下，可对本发明的示例性实施例进行各种修改和改变。

Claims (5)

1.一种基于膨胀算法的岩性插值方法，其特征在于，所述岩性插值方法包括以下步骤：
对工区进行网格化，确保每个炮点都在网格范围内；
用炮点的岩性填充炮点所在的网格；
对每个炮点进行膨胀，每个炮点的每一轮膨胀均以与上一轮膨胀被岩性填充的网格为基准点，将与所述基准点相邻且未被填充岩性的网格填充为与炮点相同的岩性，其中，第一轮膨胀的基准点为炮点；
每一轮膨胀结束后，判断工区的网格是否均被填充岩性，如果存在未被填充岩性的网格，则继续进行下一轮膨胀，直到整个工区的网格均被填充岩性。
2.根据权利要求1所述的基于膨胀算法的岩性插值方法，其特征在于，按炮点线号和点号的顺序完成所述对每个炮点进行膨胀的步骤。
3.根据权利要求2所述的基于膨胀算法的岩性插值方法，其特征在于，在所述用炮点的岩性填充炮点所在的网格的步骤中，如果同一个网格出现多个炮点，则用第一个炮点的岩性填充网格。
4.根据权利要求1所述的基于膨胀算法的岩性插值方法，其特征在于，每个炮点的膨胀完成后再进行下一个炮点的膨胀。
5.根据权利要求1所述的基于膨胀算法的岩性插值方法，其特征在于，根据工区的大小和计算精度对工区进行网格化。
CN201610686362.XA 2016-08-19 2016-08-19 一种基于膨胀算法的岩性插值方法 Active CN106372382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610686362.XA CN106372382B (zh) 2016-08-19 2016-08-19 一种基于膨胀算法的岩性插值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610686362.XA CN106372382B (zh) 2016-08-19 2016-08-19 一种基于膨胀算法的岩性插值方法

Publications (2)

Publication Number Publication Date
CN106372382A CN106372382A (zh) 2017-02-01
CN106372382B true CN106372382B (zh) 2019-03-29

Family

ID=57878103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610686362.XA Active CN106372382B (zh) 2016-08-19 2016-08-19 一种基于膨胀算法的岩性插值方法

Country Status (1)

Country Link
CN (1) CN106372382B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101388967A (zh) * 2008-10-20 2009-03-18 四川虹微技术有限公司 一种用于视图合成的空洞填充方法
CN101482970A (zh) * 2009-01-08 2009-07-15 上海交通大学 图像处理中的改进形态学滤波方法
CN102073997A (zh) * 2011-01-12 2011-05-25 东华理工大学 小型文档图像透视复原方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2918418C (en) * 2013-08-16 2020-02-18 Landmark Graphics Corporation Converting reserve estimates in a reservoir model to a standard format for dynamic comparison

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101388967A (zh) * 2008-10-20 2009-03-18 四川虹微技术有限公司 一种用于视图合成的空洞填充方法
CN101482970A (zh) * 2009-01-08 2009-07-15 上海交通大学 图像处理中的改进形态学滤波方法
CN102073997A (zh) * 2011-01-12 2011-05-25 东华理工大学 小型文档图像透视复原方法

Also Published As

Publication number Publication date
CN106372382A (zh) 2017-02-01

Similar Documents

Publication Publication Date Title
US10422923B2 (en) Systems and methods for modeling fracture networks in reservoir volumes from microseismic events
CN104504754B (zh) 一种油气储层多点统计建模的方法及装置
CN110346654B (zh) 基于普通克里金插值的电磁频谱地图构建方法
CN104850682A (zh) 基于位置的多点地质统计学建模方法
CN103217710B (zh) 获取共偏移距共方位角域道集的方法及系统
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN104252549A (zh) 一种基于克里金插值的分析布井方法
CN103529486B (zh) 一种地球化学异常圈定方法
CN105626023A (zh) 低渗透油藏垂直压裂裂缝方位试井确定方法
CN105093319A (zh) 基于三维地震数据的地面微地震静校正方法
CN104181593B (zh) 一种三维无射线追踪回折波层析成像方法及装置
CN105549077A (zh) 基于多级多尺度网格相似性系数计算的微震震源定位方法
CN107730582B (zh) 基于海洋遥感数据的海浪三维显示方法
CN105005075B (zh) 基于地震频率信息的多波匹配方法
US20200249375A1 (en) Process for defining the locations of a plurality of wells in a field, related system and computer program product
CN106908836B (zh) 采集脚印压制方法及系统
CN105301638B (zh) 一种提取风化层底界面的方法和装置
US20150032377A1 (en) System and method for remaining resource mapping
CN104573339B (zh) 页岩气储层的地质参数确定方法和装置
CN106372382B (zh) 一种基于膨胀算法的岩性插值方法
CN102830430B (zh) 一种层位速度建模方法
CN110287601B (zh) 一种毛竹胸径年龄二元联合分布精确估算方法
CN103886129A (zh) 将测井数据离散到储层网格模型的方法和装置
CN105223630B (zh) 基于地质模型的全方位观测系统参数论证方法
CN114236624B (zh) 基于电磁法估算压裂改造空间体积的方法和系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant after: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20180402

Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Applicant after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201109

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.