CN104360386B - 一种针对花岗岩地层划分对比的面元法 - Google Patents
一种针对花岗岩地层划分对比的面元法 Download PDFInfo
- Publication number
- CN104360386B CN104360386B CN201410737737.1A CN201410737737A CN104360386B CN 104360386 B CN104360386 B CN 104360386B CN 201410737737 A CN201410737737 A CN 201410737737A CN 104360386 B CN104360386 B CN 104360386B
- Authority
- CN
- China
- Prior art keywords
- attribute
- division
- sensitive
- well
- substratum
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种针对花岗岩地层划分对比的面元法，包括以下步骤：1)敏感地震属性选取；2)垂向等比例划分地震属性体；3)提取各小层敏感属性平面图；4)地层对比面元划分；5)以面元为范围，进行垂向地层格架划分；6)按顺序连接各井关键分隔点，完成地层划分对比。本发明适用于花岗岩储层划分与对比，对使用者几乎无经验要求，并且具有较强的可操作性。
Description
技术领域
本发明涉及一种地层划分对比方法，特别是一种针对花岗岩地层划分对比的面元法。
背景技术
地层作为地质体，具有多方面的特征，如地质特征(矿物成分、化学组分、岩石的结构和构造，层里层面特征、古生物化石)、地球物理特征(地磁学特征、地震波反射特征、测井曲线特征、放射性同位素等)、地球化学特征(化学元素、分子化石等)和年龄特征等。在不同的地质及资料条件下，往往根据地层各自的特征，来开展地层划分和对比。地层对比主要应用多井资料进行，即多井对比。多井对比的主要问题是：虽然地层本身具有侧向连续性，但由于井间具有较大的距离(往往数百米以上)，难于直接追踪其连续性，可依据的资料主要为岩心和测井资料。目前国内基于常规碎屑岩和碳酸盐地层已经探索出了诸如地震地层学，层序地层学，年代地层学等多套划分对比方法，并形成了相应的理论体系及技术手段。近年来，与花岗岩地层相关的石油勘探开发行为日益频繁，但是花岗岩地层发育规律与现有沉积岩地层差异较大，相关地层对比理论和技术很不成熟，致使相关地层划分对比工作对地质经验的依赖性非常强，且具有较强的主观性和不确定性。因此明确提出花岗岩地层划分对比方法非常必要。
发明内容
针对上述问题，本发明的目的是提供一种基于地震属性信息的针对花岗岩地层划分对比的面元法。
为实现上述目的，本发明采用以下技术方案：一种针对花岗岩地层划分对比的面元法，它包括以下步骤：1)敏感地震属性选取：通过对目标油田的不同地震属性的比较，选出敏感的地震属性，并计算出敏感地震属性的时间厚度。敏感的地震属性是指能够反映地下地质特征的属性，其选取标准为井点测井电性参数值与地震属性值相符合；2)垂向等比例划分地震属性体：以目的层顶面为基准面，以三维地震数据采样率作为间隔，向下将敏感地震属性体分为若干小层，分别记为n1，n2，…，ni(1≤i≤+∞)；3)提取各小层敏感属性平面图：以所划分的小层为计算单元，计算各小层敏感地震属性值，形成各自小层对应的平面图；4)地层对比面元划分：在地层平面上，以各钻井为中心、以小于四分之一井间距的长度为半径，划分出数量与井相同的面元；5)以面元为范围，进行垂向地层格架划分：在面元范围内，根据各小层自上而下属性变化规律，在垂向上将敏感地震属性特征相似的相邻小层进行合并，形成规模较大的若干大地层F1，F2，…，并记录各井垂向地层分层点，也称作关键分隔点；6)按顺序连接各井关键分隔点，完成地层划分对比：按序连接各井关键分隔点，即完成花岗岩地层划分与对比。
所述步骤5)中，敏感地震属性的特征包括界线深浅、连续性、长度、数量以及分布，对于相邻两层地震属性平面图，若相邻两层属性中的界线深浅、连续性、长度、数量和分布之一的特征比值介于0.7～1.4，即判定两层的敏感地震属性特征相似。
本发明由于采取以上技术方案，其具有以下优点：1、本发明适用于花岗岩储层划分与对比，对使用者几乎无经验要求，并且具有较强的可操作性。2、本发明方法实施效率高，应该本方法进行地层划分对比可将现有方法的效率提高一倍以上。
附图说明
图1是实施例中P油田的n1小层对应的敏感属性时间切片示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明一种针对花岗岩地层划分对比的面元法，其包括以下步骤(以P油田为例进行说明)：
1)敏感地震属性选取
通过对目标油田的不同地震属性的比较，选出敏感的地震属性，并计算出敏感地震属性的时间厚度。敏感的地震属性是指能够反映地下地质特征的属性，其选取标准为井点测井电性参数值与地震属性值相符合。
在本实施例中，以P油田为例，选取蚂蚁体属性作为基本分析手段，经过斯伦贝谢公司的商业软件Geoframe软件中LPM模块计算，得到该地震属性体垂向上时间厚度为46ms。
2)垂向等比例划分地震属性体
以目的层顶面为基准面，以三维地震数据采样率作为间隔，向下将敏感地震属性体分为若干小层，分别记为n1，n2，…，ni(1≤i≤+∞)。
在本实施例中，P油田三维地震敏感属性(蚂蚁体属性)数据体采样率为2ms。从花岗岩顶面开始，以2ms为时间间隔，自上而下共划分了n1、n2、……n23等23个平行于顶面的小层。其中，n1为T8～T8+2ms、n2为T8+2～T8+4ms、……、n23为T8+44～T8+46ms范围内地层(T代表地震解释的标准同相轴，T后面的数字代表底层年代，T8代表中生界顶面；对于本实施例P油田而言，其顶面为T8面)。
3)提取各小层敏感属性平面图
以所划分的小层为计算单元，计算各小层敏感地震属性值，形成各自小层对应的平面图。
在本实施例中，P油田共划分为了23个小层，以该23个小层为计算单元，计算各个小层敏感属性值，将敏感属性(蚂蚁体属性)值作平面插值形成各自小层对应的地震属性平面图。图1给出了P油田n1小层对应的敏感地震属性(蚂蚁体属性)平面图，从图中，可以看到工区范围内尤其是井点周围属性特征，以井1为例，可以明确属性强度、界线连续性、界线延伸长度、界线数量以及界线分布特征等，n2～n23小层的敏感地震属性时间切片可同样得到上述信息。
4)地层对比面元划分
在地层平面上，以各钻井为中心、以小于四分之一井间距的长度为半径，划分出数量与井相同的面元。
在本实施例中，P油田据此在平面上划分为P1、P2、……、P19等共19个面元，每个面元上有1口钻井。
5)以面元为范围，进行垂向地层格架划分
在面元范围内，根据各小层自上而下属性变化规律，在垂向上将敏感地震属性特征相似的相邻小层进行合并，形成规模较大的若干大地层F1，F2，…，并记录各井垂向地层分层点，也称作关键分隔点。
敏感地震属性的特征包括界线深浅、连续性、长度、数量以及分布；对于相邻两层地震属性平面图，若相邻两层属性中的任一特征的比值介于0.7～1.4，即可判定属性特征相似，则可以合并为一层。
在本实施例中，P油田面元P1垂向上小层n1、n2、n3、……、n9属性相似，可以合并为地层F1；类似，面元P1垂向上小层n10～n15可合并为大地层F2；面元P1垂向上小层n16～n23可合并为大地层F3。面元P1地层F1、F2、F3顶、底面在1井井点上自上而下有4个关键分隔点，记为P1T0、P1F1、P1F2、P1F3。
类似地，P油田各面元均可划分为F1、F2、F3等3套大地层，其关键分隔点分别记为P2T0、P2F1、P2F2、P2F3；……；P19T0、P19F1、P19F2、P19F3。
6)按顺序连接各井关键分隔点，完成地层划分对比
按序连接各井关键分隔点，即完成花岗岩地层划分与对比。
在本实施例中，将P1T0、P2T0、……P23T0相连，作为第一个地层界面；将P1F1、P2F1、……P23F1相连，作为第二个地层界面；将P1F2、P2F2、……P23F2相连，作为第三个地层界面；将P1F3、P2F3、……P23F3相连，作为第四个地层界面。至此，将P油田划分为4套地层。
本发明仅以上述实施例进行说明，各部件的结构、设置位置及其连接都是可以有所变化的，在本发明技术方案的基础上，凡根据本发明原理对个别部件进行的改进和等同变换，均不应排除在本发明的保护范围之外。
Claims (2)
1.一种针对花岗岩地层划分对比的面元法，它包括以下步骤：
1)敏感地震属性选取：通过对目标油田的不同地震属性的比较，选出敏感的地震属性，并计算出敏感地震属性的时间厚度，敏感的地震属性是指能够反映地下地质特征的属性，其选取标准为井点测井电性参数值与地震属性值相符合；
2)垂向等比例划分地震属性体：以目的层顶面为基准面，以三维地震数据采样率作为间隔，向下将敏感地震属性体分为若干小层，分别记为n1，n2，…，ni(1≤i≤+∞)；
3)提取各小层敏感属性平面图：以所划分的小层为计算单元，计算各小层敏感地震属性值，形成各自小层对应的平面图；
4)地层对比面元划分：在地层平面上，以各钻井为中心、以小于四分之一井间距的长度为半径，划分出数量与井相同的面元；
5)以面元为范围，进行垂向地层格架划分：在面元范围内，根据各小层自上而下属性变化规律，在垂向上将敏感地震属性特征相似的相邻小层进行合并，形成规模较大的若干大地层F1，F2，…，并记录各井垂向地层分层点，也称作关键分隔点；
6)按顺序连接各井关键分隔点，完成地层划分对比：按序连接各井关键分隔点，即完成花岗岩地层划分与对比。
2.如权利要求1所述的一种针对花岗岩地层划分对比的面元法，其特征在于：所述步骤5)中，敏感地震属性的特征包括界线深浅、连续性、长度、数量以及分布，对于相邻两层地震属性平面图，若相邻两层属性中的界线深浅、连续性、长度、数量和分布之一的特征比值介于0.7～1.4，即判定两层的敏感地震属性特征相似。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410737737.1A CN104360386B (zh) | 2014-12-04 | 2014-12-04 | 一种针对花岗岩地层划分对比的面元法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410737737.1A CN104360386B (zh) | 2014-12-04 | 2014-12-04 | 一种针对花岗岩地层划分对比的面元法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104360386A CN104360386A (zh) | 2015-02-18 |
| CN104360386B true CN104360386B (zh) | 2016-09-21 |
Family
ID=52527666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410737737.1A Active CN104360386B (zh) | 2014-12-04 | 2014-12-04 | 一种针对花岗岩地层划分对比的面元法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104360386B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105089658B (zh) * | 2015-07-01 | 2018-04-06 | 中国石油天然气股份有限公司 | 基于不确定度的地层对比方法及装置 |
| CN105807318B (zh) * | 2016-03-11 | 2019-01-04 | 中国地质大学（北京） | 以高分辨率层序地层学为基础的油层划分与对比方法 |
| CN106094051B (zh) * | 2016-06-08 | 2017-12-15 | 中国海洋石油总公司 | 一种花岗岩地层划分对比方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5892732A (en) * | 1996-04-12 | 1999-04-06 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| US6092026A (en) * | 1998-01-22 | 2000-07-18 | Bp Amoco Corporation | Seismic signal processing and exploration |
| CN102096098A (zh) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | 地震属性提取方法 |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
| CN103529474A (zh) * | 2012-07-02 | 2014-01-22 | 中国石油化工股份有限公司 | 采用岩性细分实现岩相精细描述的方法 |
| CN104181591A (zh) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | 多维地震属性融合方法及装置 |
-
2014
- 2014-12-04 CN CN201410737737.1A patent/CN104360386B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5892732A (en) * | 1996-04-12 | 1999-04-06 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| US6092026A (en) * | 1998-01-22 | 2000-07-18 | Bp Amoco Corporation | Seismic signal processing and exploration |
| CN102096098A (zh) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | 地震属性提取方法 |
| CN103529474A (zh) * | 2012-07-02 | 2014-01-22 | 中国石油化工股份有限公司 | 采用岩性细分实现岩相精细描述的方法 |
| CN103454686A (zh) * | 2013-08-12 | 2013-12-18 | 中国石油天然气股份有限公司 | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 |
| CN104181591A (zh) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | 多维地震属性融合方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| High-resolution sequence stratigraphic division and distribution of tidal deposits in the Zhuhai Formation, Huizhou Sag, Pearl River Mouth Basin, South China Sea;WENBO ZHENG等;《Geological Magazine》;20121231;第149卷(第4期);第722-728页 * |
| 复杂三维地层模型地层划分研究;夏艳华 等;《安徽理工大学学报（自然科学版）》;20110930;第31卷(第3期);第35-38页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104360386A (zh) | 2015-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jenkins et al. | Quantifying and predicting naturally fractured reservoir behavior with continuous fracture models | |
| Malinowski et al. | Deep seismic reflection profile in Central Europe reveals complex pattern of Paleozoic and Alpine accretion at the East European Craton margin | |
| Dolberg et al. | Porosity prediction from seismic inversion, Lavrans field, Halten Terrace, Norway | |
| Treadgold et al. | Eagle Ford shale prospecting with 3D seismic data within a tectonic and depositional system framework | |
| Pranter et al. | Sandstone-body dimensions in a lower coastal-plain depositional setting: lower Williams Fork Formation, Coal Canyon, Piceance Basin, Colorado | |
| Bergmann et al. | Combination of seismic reflection and constrained resistivity inversion with an application to 4D imaging of the CO 2 storage site, Ketzin, Germany | |
| Gernigon et al. | Magnetic expression of salt diapir–related structures in the Nordkapp Basin, western Barents Sea | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| Tyiasning et al. | A comparison of diffraction imaging to incoherence and curvature | |
| Bawazer et al. | Characterization of a fractured basement reservoir using high-resolution 3D seismic and logging datasets: A case study of the Sab’atayn Basin, Yemen | |
| Lyu et al. | Fault damage zone characterization in tight-oil sandstones of the Upper Triassic Yanchang Formation in the southwest Ordos Basin, China: Integrating cores, image logs, and conventional logs | |
| CN104360386B (zh) | 一种针对花岗岩地层划分对比的面元法 | |
| Tetyukhina et al. | Acoustic nonlinear full-waveform inversion on an outcrop-based detailed geological and petrophysical model (Book Cliffs, Utah) | |
| Gouveia et al. | Jotun 4D: Characterization of Fluid Contact Movement from Time-lapse Seismic and Production-logging-tool Data | |
| Cadena et al. | Seismic and sequence stratigraphic interpretation of the area of influence of the Magdalena submarine fan, offshore northern Colombia | |
| Chopra et al. | Introduction to this special section: Reservoir characterization | |
| Tari et al. | Seismic expression of salt tectonics in the Sab’atayn Basin, onshore Yemen | |
| Shahin et al. | Optimal dynamic rock-fluid physics template validated by petroelastic reservoir modeling | |
| Toosy et al. | Depositional architecture of Eocene and Cretaceous reservoirs (blending seismic and well data) in Mubarik Area, Central Indus Basin, Pakistan | |
| Christensen et al. | Seismically driven reservoir characterization using an innovative integrated approach: Syd Arne Field | |
| Weston Bellman | Incorporating 3C seismic data quantitatively for enhanced geologic detail in an oil sands reservoir | |
| Meng et al. | A case study of complex carbonate reservoir connectivity analysis, Tarim Basin, China | |
| Ge et al. | Facies delineation and sandstone prediction using seismic sedimentology and seismic inversion in the Eocene Huizhou Depression, Pearl River Mouth Basin, China | |
| Zhang et al. | Detailed characterization of a multilayered coalbed methane field using high-resolution sequence stratigraphy: examples from the Surat Basin in Australia | |
| Wei et al. | Interpretation of multicomponent seismic data across Wister geothermal field, Imperial Valley, California |
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 | ||
| 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 |