CN105629304B - 一种基于多属性的砂体叠合模式识别方法 - Google Patents

一种基于多属性的砂体叠合模式识别方法 Download PDF

Info

Publication number
CN105629304B
CN105629304B CN201511005038.9A CN201511005038A CN105629304B CN 105629304 B CN105629304 B CN 105629304B CN 201511005038 A CN201511005038 A CN 201511005038A CN 105629304 B CN105629304 B CN 105629304B
Authority
CN
China
Prior art keywords
sand body
work area
kinds
overlapping
characteristic sequence
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
CN201511005038.9A
Other languages
English (en)
Other versions
CN105629304A (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 Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
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 China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201511005038.9A priority Critical patent/CN105629304B/zh
Publication of CN105629304A publication Critical patent/CN105629304A/zh
Application granted granted Critical
Publication of CN105629304B publication Critical patent/CN105629304B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6169—Data from specific type of measurement using well-logging

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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于多属性的砂体叠合模式识别方法，包括以下步骤：1)结合目标工区的地质信息参数，分别建立五种叠合模式的地球物理理论模型，并分别进行模型正演，得到对应的合成地震数据；2)对五种叠合模式所得到的地震数据分别提取地震属性；3)选取一部分地震属性作为特征属性，并分别建立五种叠合模式的特征序列模板；4)将实际工区中的不连续砂体分离出来，并提取与建立特征序列模板所用的相同的地震属性；5)将五种不同的特征序列模板分别与实际工区分离出并提取的地震属性进行相关性计算，求解出相关系数；6)取相关系数最大的一个特征序列模板，该特征序列模板所对应的叠合模式即为实际工区的砂体叠合模式。

Description

一种基于多属性的砂体叠合模式识别方法
技术领域
本发明涉及一种基于多属性的砂体叠合模式识别方法，属于石油勘探与开发技术领域。
背景技术
随着油田工区勘探开发的成熟度越来越高，建立三维储层定量地质模型成为储层研究的重点目标，这也说明储层研究发展的阶段更高了。对油藏进行科学的评价、开发管理、以及对油藏进行三维模拟均要求有一个比较准确的三维储层地质模型，即表征储层的地质属性在三维空间上的变化与分布的数字化模型，这一模型具有常规储层地质的二维图件无可比拟的优点。而随着开发的成熟度越来越高，新发现的油藏大多具有以下特点：构造油气藏越来越难找，油气藏的规模越来越小，储层岩性的横向变化越来越剧烈，这就大大增加了人们对地下地质构造的认识难度，也不利于后续的开发生产。面对这些状况，利用现有的各类数据来建立起高分辨率的三维储层模型，进而实现对储层分布规律及其非均质性的精细描述，是一个较为理想的应对方法。建立三维储层地质模型的其中一个重要的步骤是对砂体的分布进行识别。
当前，在地震储层研究中通过地震属性、波阻抗等信息寻找砂体分布，已经成为油气分析的重要技术。但在沙泥薄互层的情况下，受限于地震记录的分辨率，单纯基于地震数据得到的砂体反演结果只能够显示砂体的整体特征，难以了解砂体内幕。实际上，在沙泥岩互层，尤其是岩层较薄的情况下，基于地震记录识别的砂体可能是多个砂体、泥质夹层的复合体，如何基于现有数据分析砂体内幕成为当前急待解决的问题。本发明就是通过已有的地质模式作为约束条件，构建相应的砂体叠合模式；利用这些叠合模式作为先验条件，然后从地震数据中提取相应属性，通过对比属性间的相似关系，找到各砂体最佳的叠合模式，从能够了解各个砂体内部的沙泥岩分布特征。
发明内容
针对上述问题，本发明的目的是提供一种能够获知砂体内部沙泥岩分布特征的基于多属性的砂体叠合模式识别方法。
为实现上述目的，本发明采用以下技术方案：一种基于多属性的砂体叠合模式识别方法，包括以下步骤：1)结合目标工区的地质信息参数，分别建立五种叠合模式的地球物理理论模型，并分别进行模型正演，得到对应的合成地震数据；2)对五种叠合模式所得到的地震数据分别提取地震属性；3)选取一部分地震属性作为特征属性，并分别建立五种叠合模式的特征序列模板；4)将实际工区中的不连续砂体分离出来，并提取与建立特征序列模板所用的相同的地震属性；5)将五种不同的特征序列模板分别与实际工区分离出并提取的地震属性进行相关性计算，求解出相关系数；6)取相关系数最大的一个特征序列模板，该特征序列模板所对应的叠合模式即为实际工区的砂体叠合模式。
所述步骤3)中，所述特征序列为A＝(A1,A2,A3,A4,A5)，A1＝(a11,a12,a13,a14)，a11＝(p1,p2,p3)，其中，Ai代表第i(i＝1,2,…,5)个模板的特征序列，aij代表第i(i＝1,2,…,5)模板的第j(j＝1,2,…,4)个特征属性，p1、p2、p3代表某个特征属性的具体数据。
所述步骤2)中所提取的地震属性包括振幅属性、频率属性、相位属性。
所述步骤1)中建立地球物理理论模型所用到的目标工区的地质信息参数包括不同岩性的速度、密度、孔隙度、厚度。
本发明由于采取以上技术方案，其具有以下优点：1、本发明方法通过识别砂体的叠合模式，从而能够判断出砂体内部的砂泥岩分布特征。2、本发明方法可以由计算机来实现，从而实现自动识别过程，识别效率高。3本发明能够对用于建模的数据提供更为严格的约束，有利于提高模型的精准度。
附图说明
图1为本发明的方法流程示意图；
图2为本发明所用到的五种砂体叠合模式的地球物理模型图；
图3为本发明中五种叠合模式模型及其合成地震记录；其中，图(a)、(b)、(c)、(d)、(e)分别代表第一、二、三、四、五种模式的模型图及合成地震数据；
图4为五种叠合模式合成地震记录中提取的最大最小振幅信息；
图5为具体实施例秦皇岛32-6油田过井沉积相分布；
图6为具体实施例秦皇岛32-6油田明砂体叠合模式反演结果；
图7为具体实施例秦皇岛32-6油田沉积相平面分布；
图8为具体实施例秦皇岛32-6油田砂体叠合模式反演结果。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图2所示，目前在地质学上，一项受到普遍认可的成果是将河道砂体简化总结为五种地质概念叠合模式(由井涌泉博士根据岳大力老师的研究成果首先提出的)：
第一种模式为楔形砂体模型，为最简单的砂体叠合模式，模型中只有单个砂体且砂体的厚度不断变化；
第二种模式为两个砂体组合，砂体厚度恒定但砂体间的泥岩厚度不断变化；
第三种模式为两个楔形体组合，砂体间泥岩厚度恒定；
第四种为三个砂体组合，砂体间为泥岩，中间砂体的位置不断变化；
第五种为四个砂体组合，砂体间为泥岩，中间两个砂体的位置不断变化。
本发明的是以上述理论成果为基础而实施的，本发明的步骤如下(如图1所示)：
1)结合目标工区的地质信息参数(包括工区地下底层中砂岩、泥岩、火山岩等不同岩性的速度、密度、孔隙度、厚度等)，建立五种叠合模式的地球物理理论模型，并分别进行模型正演，得到对应的合成地震数据。
模型正演时首先对模型进行赋值参数：砂岩的密度为2.1g/cm3，砂岩的速度为2000m/s；泥岩的密度为2.3g/cm3，泥岩的速度为2580m/s。正演时用到的地震子波为雷克子波，频率为40Hz。
2)对五种叠合模式所得到的地震数据分别提取地震属性(如图3所示)，包括振幅属性、频率属性、相位属性等。
3)选取一部分地震属性作为特征属性，并分别建立五种叠合模式的特征序列模板。
将特征属性转化为特征序列，即指将选取的特征属性按顺序放到同一个一维序列中。为了减少用于叠合模式识别时的计算量，可以对选取的特征属性进行抽取，如每个属性选择三个数作为该特征属性的特征值，依次放到同一个一维序列中。如假设有四个特征属性，每个特征属性我们用三个数据来表征，那么最终的一维特征序列为A＝(A1,A2,A3,A4,A5)，A1＝(a11,a12,a13,a14)，a11＝(p1,p2,p3)，其中，Ai代表第i(i＝1,2,…,5)个模板的特征序列，aij代表第i(i＝1,2,…,5)模板的第j(j＝1,2,…,4)个特征属性，p1、p2、p3代表某个特征属性的具体数据。
4)将实际工区中的不连续砂体分离出来，并提取与建立特征序列模板所用的相同的地震属性。
5)将五种不同的特征序列模板分别与实际工区分离出并提取的地震属性进行相关性计算，并求解出相关系数，其中，相关系数的计算公式如下：
式中，xi表示目标砂体数据，n表示采样点数，yi表示模板数据。
6)取相关系数最大的一个特征序列模板，该特征序列模板所对应的叠合模式即为实际工区的砂体叠合模式。
至此，完成了对实际工区砂体的叠合模式的识别，一旦识别出了实际工区砂体的叠合模式，那么砂体内部的特征也就清楚了。
下面以秦皇岛32-6油田为具体实施例，来具体说明本发明的基于多属性的砂体叠合模式识别方法流程。
该油田构造位于渤中坳陷石臼坨凸起中西部，凸起周边被渤中、秦南和南堡三大富油凹陷所环绕，是渤海海域油气富集最有利地区之一，构造面积约110km2，含油面积39.7km2，东南距427油田约20km，东距428西油田42km，西北距京塘港约20km。根据资料显示，在凸起石臼坨上，发育着两套含油层系，自下而上分别为：前第三系潜山地层与上覆新生界地层构成储盖组合；由新生界储盖组合构成的含油层系，其中包括四个层组，分别为：明化镇、馆陶、东营及沙河街组。该油田工区钻遇了第二套含油层系中的前三套层组，其中明化镇组下段为该油田的主力含油层段。为了方便研究，我们在该油田北区划定一定范围作为研究的目标工区。划定的目标工区面积约13.5km2，工区内有A、B两个平台，共46口井，其中包括两口直井。工区范围为线号范围196至366，道号范围481至737，采样时间范围为900至1300ms，采样间隔为2ms。目的层段为新近系的明化镇组，包括明化镇1、明化镇2、明化镇3和明化镇4共4个油组，地震解释层位8层。
图5为秦皇岛32-6油田目标工区过井的沉积相分布图，图6为砂体叠合模式的反演结果。对比图5和图6二者之间整体相态分布较为相似。由于沉积相分析基于测井数据，而砂体叠合模式依据地震数据进行计算，因此二者在局部位置存在差异。
此外，沉积相分析虽然能够划分不同相带，但同一沉积相内往往视为均质；与之相比，砂体叠合模式反演方法能够得到砂体的叠合类型，且充分考虑所在位置的地震属性信息，因此能够更加精确的刻画砂体的非均质性。
如图7所示为秦皇岛32-6油田目标工区沉积相的平面分布图，图8为砂体叠合模式的平面分布图，二者之间存在较大的差异。这种差异来自于所用信息的不同，基于测井的沉积相分析与基于地震的砂体叠合模式分析由于所有数据的差异，必然会导致结果的不对应。沉积相分析从井信息出发，在井附近能够很好的吻合，但在远离井的位置，其可靠性将受到挑战。反观砂体叠合模式反演，虽然其在井附近难以与井完全匹配，由于其依靠地震信息进行分析研究，因此在远井位置依然具有较高的可靠性。
本发明仅以上述实施例进行说明，各部件的结构、设置位置及其连接都是可以有所变化的，在本发明技术方案的基础上，凡根据本发明原理对个别部件进行的改进和等同变换，均不应排除在本发明的保护范围之外。

Claims (4)

1.一种基于多属性的砂体叠合模式识别方法，包括以下步骤：
1)结合目标工区的地质信息参数，分别建立五种叠合模式的地球物理理论模型，并分别进行模型正演，得到对应的合成地震数据；
2)对五种叠合模式所得到的地震数据分别提取地震属性；
3)选取一部分地震属性作为特征属性，并分别建立五种叠合模式的特征序列模板；
4)将实际工区中的不连续砂体分离出来，并提取与建立特征序列模板所用的相同的地震属性；
5)将五种不同的特征序列模板分别与实际工区分离出并提取的地震属性进行相关性计算，求解出相关系数；
6)取相关系数最大的一个特征序列模板，该特征序列模板所对应的叠合模式即为实际工区的砂体叠合模式。
2.如权利要求1所述的一种基于多属性的砂体叠合模式识别方法，其特征在于：所述步骤3)中，所述特征序列为A＝(A1,A2,A3,A4,A5)，A1＝(a11,a12,a13,a14)，a11＝(p1,p2,p3)，其中，Ai代表第i，i＝1,2,…,5个模板的特征序列，aij代表第i，i＝1,2,…,5模板的第j，j＝1,2,…,4个特征属性，p1、p2、p3代表所选取的特征属性的具体数据。
3.如权利要求1所述的一种基于多属性的砂体叠合模式识别方法，其特征在于：所述步骤2)中所提取的地震属性包括振幅属性、频率属性、相位属性。
4.如权利要求1所述的一种基于多属性的砂体叠合模式识别方法，其特征在于：所述步骤1)中建立地球物理理论模型所用到的目标工区的地质信息参数包括不同岩性的速度、密度、孔隙度、厚度。
CN201511005038.9A 2015-12-29 2015-12-29 一种基于多属性的砂体叠合模式识别方法 Active CN105629304B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511005038.9A CN105629304B (zh) 2015-12-29 2015-12-29 一种基于多属性的砂体叠合模式识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511005038.9A CN105629304B (zh) 2015-12-29 2015-12-29 一种基于多属性的砂体叠合模式识别方法

Publications (2)

Publication Number Publication Date
CN105629304A CN105629304A (zh) 2016-06-01
CN105629304B true CN105629304B (zh) 2017-10-03

Family

ID=56044415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511005038.9A Active CN105629304B (zh) 2015-12-29 2015-12-29 一种基于多属性的砂体叠合模式识别方法

Country Status (1)

Country Link
CN (1) CN105629304B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106707340A (zh) * 2016-12-13 2017-05-24 中国石油天然气股份有限公司大港油田分公司 一种火山岩岩相的预测方法
CN106646606B (zh) * 2016-12-14 2018-08-10 中国石油大学(华东) 基于地震特征参数模式识别的薄砂体表征方法
CN106842317B (zh) * 2016-12-22 2019-02-15 中国石油天然气集团公司 一种预测油砂体分布的方法及装置
CN110727027B (zh) * 2018-07-17 2021-11-02 中国石油天然气股份有限公司 多期次河道砂的精细刻画方法
CN109358323B (zh) * 2018-11-05 2022-01-21 中国矿业大学(北京) 一种地质雷达数据解释方法、装置，电子设备及存储介质
CN111399040B (zh) * 2020-03-06 2021-11-16 西南石油大学 一种基于地震属性负差异特征的叠置砂体识别模型及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0531515B1 (en) * 1991-03-27 1995-08-23 Exxon Production Research Company Displaying n dimensional seismic data in an n-1 dimensional format
CN104991275B (zh) * 2015-07-09 2018-03-09 北京诺克斯达石油科技有限公司 一种特征切片薄互层分析法

Also Published As

Publication number Publication date
CN105629304A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN105629304B (zh) 一种基于多属性的砂体叠合模式识别方法
CN107817535B (zh) 薄夹层的确定方法和装置
CN101158724B (zh) 基于偶极小波的储层厚度预测方法
Barwis et al. Sandstone petroleum reservoirs
CN105044770A (zh) 致密砂砾岩气藏储层定量预测方法
CN105093306B (zh) 一种地球物理勘探中储层自动解释与厚度求取方法
CN109870719B (zh) 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统
CN103513278B (zh) 利用地震波组厚度进行储层预测的方法
CN107515957A (zh) 泥页岩层序地层划分方法
CN103217714A (zh) 基于频率-尺度匹配的地震沉积学解释方法
CN109164487A (zh) 一种基于模型建立平均速度场的方法及精细成图的方法
CN109541685A (zh) 一种河道砂体识别方法
CN106570262A (zh) 一种储层构型结构的描述方法
Schneider et al. Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA
CN105572726A (zh) 一种层序地层格架下多期次叠置扇体的精细刻画方法
CN106353806B (zh) 一种检测生物礁储层连通性的方法
CN106526697B (zh) 一种三角洲建设与破坏主要转换面识别方法
CN106338764A (zh) 生物启发计算地层圈闭油气藏超剥线识别方法
CN106033127A (zh) 基于横波速度变化率的地应力方位地震预测方法
CN103278852B (zh) 利用地震数据体波形结构特征模型预测油气的方法
Singleton Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas
CN104360386A (zh) 一种针对花岗岩地层划分对比的面元法
Curtis et al. Gulf Coast Cenozoic: a model for the application of stratigraphic concepts to exploration on passive margins
CN105093330B (zh) 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法
CN109765621B (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
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC Research Institute

Patentee before: China National Offshore Oil Corporation