CN105717542A - 基于三维地震的Wheeler域切片砂体雕刻方法 - Google Patents

基于三维地震的Wheeler域切片砂体雕刻方法 Download PDF

Info

Publication number
CN105717542A
CN105717542A CN201610220910.XA CN201610220910A CN105717542A CN 105717542 A CN105717542 A CN 105717542A CN 201610220910 A CN201610220910 A CN 201610220910A CN 105717542 A CN105717542 A CN 105717542A
Authority
CN
China
Prior art keywords
sand body
wheeler
seismics
territory
layer position
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.)
Granted
Application number
CN201610220910.XA
Other languages
English (en)
Other versions
CN105717542B (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.)
Beijing Co-Sail Petroleum Technology CoLtd
Original Assignee
Beijing Co-Sail Petroleum Technology CoLtd
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 Beijing Co-Sail Petroleum Technology CoLtd filed Critical Beijing Co-Sail Petroleum Technology CoLtd
Priority to CN201610220910.XA priority Critical patent/CN105717542B/zh
Publication of CN105717542A publication Critical patent/CN105717542A/zh
Application granted granted Critical
Publication of CN105717542B publication Critical patent/CN105717542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase

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

本发明公开了一种基于三维地震的Wheeler域切片砂体雕刻方法,包括以下步骤:井震结合确定层序顶底界面,计算地震同相轴倾角导向体,根据倾角导向是否突变追踪年代地层,年代地层按沉积顺序排列,代表相对地质时间。追踪年代地层的参数为,数据驱动,所有层都连续,最大迭代为2次;突变点打断,迭代10次,建立起年代地层三维层位体。将层位体的所有小层全部拉平,转换到Wheeler域。Wheeler域切片,能够反映砂体的平面分布范围,对砂体进行空间雕刻,Wheeler域地震属性切片,还能反映砂体内部地震属性的平面分布,研究砂体分布特征。该方法可以提高冲积扇砂体、扇三角洲砂体识别的准确性,提高钻井成功率。

Description

基于三维地震的Wheeler域切片砂体雕刻方法
技术领域
本发明属于砂体雕刻技术领域,具体地涉及一种基于三维地震的Wheeler域切片砂体雕刻方法。
背景技术
冲积扇砂体、扇三角洲砂体识别的原理是基于砂体的密度和速度与周围的泥岩有较明显的差异,从而每一期的砂体反映在地震资料上有一定的特征,以往通过反演、地震属性分析方法对砂体的识别,高度依赖钻井资料,在离井较远的区域,误差很大,离井越远误差越大,钻井成功率仅为60-70%,因此还有30-40%的钻井成功率有待提高。
对于冲积扇和扇三角洲为储层的油藏,高效井一般只分布在主河道相带内,河道间为致密储层分布区,可作侧向遮挡,为了精细刻画冲积扇和扇三角洲的主河道砂体的分布,迫切需要砂体雕刻技术,以进一步提高钻探成功率。而当前国内外对于岩性横向变化大,储层非均质性强,但主河道内产量高的油藏,还缺乏利用地震资料研究砂体的有效手段。
发明内容
针对上述存在的技术问题,本发明旨在提供一种基于三维地震的Wheeler域切片砂体雕刻方法,可以精细刻画砂体的空间分布范围,以及砂体的储层物性特征,并且提高钻井成功率。
为达到上述目的,本发明的技术方案是:
一种基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,包括以下步骤:
S01:井震结合确定层序顶底界面;
S02:计算三维地震倾角导向体,通过自动检测三维地震同相轴,计算每个样点的倾角和方位角;
S03:计算三维地震年代地层层位体,所述层位体为一组可连续追踪的层位,按地层沉积顺序排列;
S04:将层位体的所有小层全部拉平,转到Wheeler域,再把地震数据体按小层拉平,生成Wheeler域的地震数据体;
S05:根据井震标定结果显示Wheeler域切片,确定砂体顶底平面分布范围,做出砂体顶底界面构造图;根据Wheeler域地震切片,分析砂体的地震属性特征,研究砂体的储层物性及含油性分布特征。
优选的,所述步骤S01中井震结合确定层序顶底界面,为等时层序界面,对目的层段,通过钻井合成地震记录,以及测井曲线精细层位标定,以单砂体规模的高精度层序格架进行标定。
优选的,所述步骤S02中计算三维地震倾角导向体的算法选择快速傅立叶变换法FFT,计算步长为3线、3道、2个采样点。
优选的,所述步骤S03的年代地层能够精确反映地震同相轴的形态。
优选的,所述步骤S03包括追踪层位体,所述追踪层位体包括:①追踪连续层位体,确保层位闭合;②根据倾角、方位角是否连续变化把连续的层位体打断,层位重合部分,只保留单层;③层位体补洞消除空白洞。
优选的,所述步骤S04中所述Wheeler域数据体在垂向上代表相对地质时间,表示沉积的先后顺序,最下面的地层先沉积,上面的地层后沉积。
与现有技术相比,本发明的有益效果是:
本发明基于三维地震的Wheeler域切片砂体雕刻方法可以实现砂体平面分布的空间形态以及分析砂体的储层物性分布特征的研究,结合冲积扇、扇三角洲所处的大地构造背景、物源方向、沉积中心分布特征,基准面变化规律,砂体的地震属性特点,开展地质综合研究,寻找有利储层发育区,部署勘探开发井位,可以大大提高钻井成功率,成功率由原来的60-70%提高到90%。本发明中的基于三维地震的Wheeler域切片砂体雕刻方法具有广泛的应用前景,且能够大幅度提高钻井成功率。
附图说明
图1为本发明基于三维地震的Wheeler域切片砂体雕刻方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面结合具体实施方式以及附图,对本发明进一步详细说明。应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
实施例:
图1为本发明的基于三维地震的Wheeler域切片砂体雕刻方法的流程图,如图1所示。
在步骤101,井震结合确定层序顶底界面,为等时层序界面,对目的层段,通过钻井合成地震记录,以及测井曲线精细层位标定,以单砂体规模的高精度层序格架进行标定,为砂体空间雕刻,提供可靠的依据。等时层序界面是计算层位体的终止层,仅仅是限定层位体的计算范围,对层位体的计算结果不产生任何影响。流程进入到步骤102。
在步骤102,所述步骤S02中计算三维地震倾角导向体,算法选择FFT(快速傅立叶变换法),计算地震倾角导向体有多种算法,FFT法更加精准的反映地震同相轴的空间展布形态,计算步长为332(3线、3道、2个采样点),这一组参数能够反映同相轴的变化趋势,又能兼顾砂体横向变化的细节。流程进入到步骤103。
在步骤103,计算三维地震年代地层层位体,所述层位体是一组可连续追踪的层位,按地层沉积顺序排列,层位体能够精确反映地震同相轴的特征。流程进入到步骤104。
在步骤104,追踪层位体的步骤包括:①追踪连续层位体(Continuousevents)104-1,连续层位体的所有层位都连续分布,层位可重合,目的是确保层位闭合,追踪连续层位体的参数为,起始位置为厚度最大处的第一个采样点,年代地层的间距最大为4ms,检测同相轴的最远距离50道,最大迭代次数2;②把连续层位体打断(TruncateHorizonCube)104-2,根据倾角、方位角是否连续变化把连续的层位体打断,层位重合部分,只保留单层;③层位体补洞(AdditerationstoHorizonCube)104-3,连续层位体打断后,层位重合现象被消除,但在剖面上会出现一些没有层位的空白洞,层位体补洞用于消除空白洞,两个年代地层间距最大2ms,超过2ms时增加一层,最大迭代次数10。
在步骤105,将层位体转到Wheeler域,把层位体的小层全部拉平,转到Wheeler域,Wheeler域数据体X、Y方向与常规数据体完全一致,在垂向上不代表深度也不代表时间,而是相对地质时间,表示沉积的先后顺序,最下面的地层先沉积,上面的地层后沉积,再把地震数据体按小层拉平,生成Wheeler域的各种地震属性体,包括常规数据体、反演数据体、相位数据体等。流程进入到步骤106。
在步骤106,根据井震标定结果显示Wheeler域切片,可以确定砂体顶底平面分布范围,做出砂体顶底界面构造图,同时,根据Wheeler域地震属性切片,分析砂体的地震属性特征,研究砂体的储层物性及含油性分布特征。
通过Wheeler域切片对砂体的精细雕刻,研究砂体平面分布形态,分析砂体物性特征,结合冲积扇、扇三角洲所处的大地构造背景,物源方向,沉积中心,基准面变化规律,砂体的地震属性特点,开展地质综合研究,寻找有利储层发育区,部署勘探开发井位,可以大大提高钻井成功率,成功率由原来的60-70%提高到90%。本发明中的基于三维地震的Wheeler域切片砂体雕刻方法具有广泛的应用前景,且能够大幅度提高钻井成功率。
应当理解的是,本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理,而不构成对本发明的限制。因此,在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。此外,本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (8)

1.一种基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,包括以下步骤:
S01:井震结合确定层序顶底界面;
S02:计算三维地震倾角导向体,通过自动检测三维地震同相轴,计算每个样点的倾角和方位角;
S03:计算三维地震年代地层层位体,所述层位体为一组可连续追踪的层位,按地层沉积顺序排列;
S04:将层位体的所有小层全部拉平,转到Wheeler域,再把地震数据体按小层拉平,生成Wheeler域的地震数据体;
S05:根据井震标定结果显示Wheeler域切片,确定砂体顶底平面分布范围,做出砂体顶底界面构造图;根据Wheeler域地震切片,分析砂体的地震属性特征,研究砂体的储层物性及含油性分布特征。
2.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述步骤S01中井震结合确定层序顶底界面,为等时层序界面,对目的层段,通过钻井合成地震记录,以及测井曲线精细层位标定,以单砂体规模的高精度层序格架进行标定。
3.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述步骤S02中计算三维地震倾角导向体的算法选择快速傅立叶变换法FFT,计算步长为3线、3道、2个采样点。
4.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述步骤S03的年代地层能够精确反映地震同相轴的形态。
5.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述步骤S03包括追踪层位体,所述追踪层位体包括:①追踪连续层位体,确保层位闭合;②根据倾角、方位角是否连续变化把连续的层位体打断,层位重合部分,只保留单层;③层位体补洞消除空白洞。
6.根据权利要求5所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述追踪连续层位体的参数为,起始位置为厚度最大处的第一个采样点,年代地层的间距最大为4ms,检测同相轴的最远距离50道,最大迭代次数2。
7.根据权利要求5所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述层位体补洞的参数为两个年代地层间距最大2ms,超过2ms时增加一层,最大迭代次数10。
8.根据权利要求1所述的基于三维地震的Wheeler域切片砂体雕刻方法,其特征在于,所述步骤S04中所述Wheeler域数据体在垂向上代表相对地质时间,表示沉积的先后顺序,最下面的地层先沉积,上面的地层后沉积。
CN201610220910.XA 2016-04-11 2016-04-11 基于三维地震的Wheeler域切片砂体雕刻方法 Active CN105717542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610220910.XA CN105717542B (zh) 2016-04-11 2016-04-11 基于三维地震的Wheeler域切片砂体雕刻方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610220910.XA CN105717542B (zh) 2016-04-11 2016-04-11 基于三维地震的Wheeler域切片砂体雕刻方法

Publications (2)

Publication Number Publication Date
CN105717542A true CN105717542A (zh) 2016-06-29
CN105717542B CN105717542B (zh) 2019-02-22

Family

ID=56159948

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610220910.XA Active CN105717542B (zh) 2016-04-11 2016-04-11 基于三维地震的Wheeler域切片砂体雕刻方法

Country Status (1)

Country Link
CN (1) CN105717542B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107976713A (zh) * 2017-11-08 2018-05-01 中国石油天然气股份有限公司 一种高维地震数据输入下去除沉积背景的方法及装置
CN110244359A (zh) * 2019-07-03 2019-09-17 成都理工大学 一种基于改进地震切片技术的异常体边缘检测计算方法
CN111273346A (zh) * 2018-12-04 2020-06-12 中国石油天然气股份有限公司 去除沉积背景的方法、装置、计算机设备及可读存储介质
US10984590B1 (en) 2019-12-06 2021-04-20 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
CN112711069A (zh) * 2019-10-24 2021-04-27 中国石油化工股份有限公司 基于平原区曲流河河道演化的点沙坝预测方法
US11010969B1 (en) 2019-12-06 2021-05-18 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
CN112904423A (zh) * 2019-12-03 2021-06-04 中国石油天然气股份有限公司 陆相河道的雕刻方法及装置
CN113419274A (zh) * 2021-06-21 2021-09-21 大庆油田有限责任公司 基于高精度层序格架模型的三维地震切片属性体提取方法
US11187826B2 (en) 2019-12-06 2021-11-30 Chevron U.S.A. Inc. Characterization of subsurface regions using moving-window based analysis of unsegmented continuous data
US11249220B2 (en) 2019-08-14 2022-02-15 Chevron U.S.A. Inc. Correlation matrix for simultaneously correlating multiple wells
US11263362B2 (en) 2020-01-16 2022-03-01 Chevron U.S.A. Inc. Correlation of multiple wells using subsurface representation
US11320566B2 (en) 2020-01-16 2022-05-03 Chevron U.S.A. Inc. Multiple well matching within subsurface representation
US11397279B2 (en) 2020-03-27 2022-07-26 Chevron U.S.A. Inc. Comparison of wells using a dissimilarity matrix
US11604909B2 (en) 2019-05-28 2023-03-14 Chevron U.S.A. Inc. System and method for accelerated computation of subsurface representations
CN116413794A (zh) * 2021-12-31 2023-07-11 中国石油天然气股份有限公司 碳酸盐岩的沉积相类型的预测方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009096942A1 (en) * 2008-01-28 2009-08-06 Landmark Graphics Corporation Hybrid stratigraphic layering using pseudo-wheeler space
CN103454678A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 一种地震切片等时性的确定方法及系统
EP2899571A1 (fr) * 2014-01-28 2015-07-29 IFP Energies nouvelles Procédé de construction d'un maillage volumique pour la modelisation de structures geologiques
CN105372705A (zh) * 2015-10-27 2016-03-02 中国石油天然气股份有限公司 一种基于多波资料的地层切片方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009096942A1 (en) * 2008-01-28 2009-08-06 Landmark Graphics Corporation Hybrid stratigraphic layering using pseudo-wheeler space
CN103454678A (zh) * 2013-08-12 2013-12-18 中国石油天然气股份有限公司 一种地震切片等时性的确定方法及系统
EP2899571A1 (fr) * 2014-01-28 2015-07-29 IFP Energies nouvelles Procédé de construction d'un maillage volumique pour la modelisation de structures geologiques
CN105372705A (zh) * 2015-10-27 2016-03-02 中国石油天然气股份有限公司 一种基于多波资料的地层切片方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
单蕊 等: "地震沉积学在川东北YB地区沉积相分析中的应用", 《煤田地质与勘探》 *
白雪峰 等: "倾角导向体约束下的储层预测技术", 《物探化探计算技术》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107976713B (zh) * 2017-11-08 2019-07-09 中国石油天然气股份有限公司 一种高维地震数据输入下去除沉积背景的方法及装置
CN107976713A (zh) * 2017-11-08 2018-05-01 中国石油天然气股份有限公司 一种高维地震数据输入下去除沉积背景的方法及装置
CN111273346A (zh) * 2018-12-04 2020-06-12 中国石油天然气股份有限公司 去除沉积背景的方法、装置、计算机设备及可读存储介质
US11604909B2 (en) 2019-05-28 2023-03-14 Chevron U.S.A. Inc. System and method for accelerated computation of subsurface representations
CN110244359A (zh) * 2019-07-03 2019-09-17 成都理工大学 一种基于改进地震切片技术的异常体边缘检测计算方法
CN110244359B (zh) * 2019-07-03 2022-02-25 成都理工大学 一种基于改进地震切片技术的异常体边缘检测计算方法
US11249220B2 (en) 2019-08-14 2022-02-15 Chevron U.S.A. Inc. Correlation matrix for simultaneously correlating multiple wells
CN112711069A (zh) * 2019-10-24 2021-04-27 中国石油化工股份有限公司 基于平原区曲流河河道演化的点沙坝预测方法
CN112711069B (zh) * 2019-10-24 2024-03-26 中国石油化工股份有限公司 基于平原区曲流河河道演化的点沙坝预测方法
CN112904423A (zh) * 2019-12-03 2021-06-04 中国石油天然气股份有限公司 陆相河道的雕刻方法及装置
CN112904423B (zh) * 2019-12-03 2024-04-30 中国石油天然气股份有限公司 陆相河道的雕刻方法及装置
US10984590B1 (en) 2019-12-06 2021-04-20 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
US11187826B2 (en) 2019-12-06 2021-11-30 Chevron U.S.A. Inc. Characterization of subsurface regions using moving-window based analysis of unsegmented continuous data
US11010969B1 (en) 2019-12-06 2021-05-18 Chevron U.S.A. Inc. Generation of subsurface representations using layer-space
US11320566B2 (en) 2020-01-16 2022-05-03 Chevron U.S.A. Inc. Multiple well matching within subsurface representation
US11263362B2 (en) 2020-01-16 2022-03-01 Chevron U.S.A. Inc. Correlation of multiple wells using subsurface representation
US11397279B2 (en) 2020-03-27 2022-07-26 Chevron U.S.A. Inc. Comparison of wells using a dissimilarity matrix
CN113419274A (zh) * 2021-06-21 2021-09-21 大庆油田有限责任公司 基于高精度层序格架模型的三维地震切片属性体提取方法
CN113419274B (zh) * 2021-06-21 2022-03-01 大庆油田有限责任公司 基于高精度层序格架模型的三维地震切片属性体提取方法
CN116413794A (zh) * 2021-12-31 2023-07-11 中国石油天然气股份有限公司 碳酸盐岩的沉积相类型的预测方法及装置
CN116413794B (zh) * 2021-12-31 2024-07-23 中国石油天然气股份有限公司 碳酸盐岩的沉积相类型的预测方法及装置

Also Published As

Publication number Publication date
CN105717542B (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
CN105717542A (zh) 基于三维地震的Wheeler域切片砂体雕刻方法
CN102053270B (zh) 一种基于沉积地层单元的地震相分析方法
CN102759745B (zh) 一种基于数字地质露头模型正演的碳酸盐岩储层预测方法
Nardin et al. Stratigraphic architecture of a large-scale point-bar complex in the McMurray Formation: Syncrude's Mildred Lake Mine, Alberta, Canada
EP3329307B1 (en) Assignment of systems tracts
Reesink et al. Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina
Pringle et al. Capturing stratigraphic and sedimentological complexity from submarine channel complex outcrops to digital 3D models, Karoo Basin, South Africa
CN103675944B (zh) 一种沉积盆地铀资源勘探的综合地球物理方法
CN104316978A (zh) 地球物理的近地表三维速度场研究方法和装置
Geersen et al. Active tectonics of the North Chilean marine forearc and adjacent oceanic Nazca Plate
CN102200588A (zh) 地震资料的波形相似性体曲率分析方法
CN111899338B (zh) 一种覆盖区地层岩性三维建模的方法、装置及系统
McClymont et al. Characterization of the shallow structures of active fault zones using 3‐D ground‐penetrating radar data
Weismüller et al. Structure of massively dilatant faults in Iceland: lessons learned from high-resolution unmanned aerial vehicle data
Egger et al. Structural controls on geothermal circulation in Surprise Valley, California: A re-evaluation of the Lake City fault zone
Mayall et al. Seismic interpretation workflows for deep-water systems: A practical guide for the subsurface
Zhang et al. Multi-frequencies GPR measurements for delineating the shallow subsurface features of the Yushu strike slip fault
CN110019594A (zh) 一种基于数值模拟的地质知识库建立方法
CN113031066B (zh) 一种隐蔽性断层的识别方法
Galone et al. Integrating near-surface geophysical methods and remote sensing techniques for reconstructing fault-bounded valleys (Mellieha valley, Malta)
Anomneze et al. Application of seismic stratigraphy and structural analysis in the determination of petroleum plays within the Eastern Niger Delta Basin, Nigeria
Sandersen et al. Tectonic impact on Pleistocene and Holocene erosional patterns in a formerly glaciated intra-plate area
CN102841374A (zh) 基于扫描面正演的伪三维快速微地震正演方法
CN106547022A (zh) 年代地层相识别冲积扇微相的方法
Pucci et al. Coseismic ruptures and tectonic landforms along the Düzce segment of the North Anatolian Fault Zone (Ms 7.1, November 1999)

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
GR01 Patent grant
GR01 Patent grant