CN106443791B - 求取倾斜地层或各向异性地层横波剩余静校正量的方法 - Google Patents

求取倾斜地层或各向异性地层横波剩余静校正量的方法 Download PDF

Info

Publication number
CN106443791B
CN106443791B CN201611038311.2A CN201611038311A CN106443791B CN 106443791 B CN106443791 B CN 106443791B CN 201611038311 A CN201611038311 A CN 201611038311A CN 106443791 B CN106443791 B CN 106443791B
Authority
CN
China
Prior art keywords
wave
data
static correction
converted
value
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
CN201611038311.2A
Other languages
English (en)
Other versions
CN106443791A (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 CN201611038311.2A priority Critical patent/CN106443791B/zh
Publication of CN106443791A publication Critical patent/CN106443791A/zh
Application granted granted Critical
Publication of CN106443791B publication Critical patent/CN106443791B/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/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/53Statics correction, e.g. weathering layer or transformation to a datum

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

本发明提供了一种求取倾斜地层或各向异性地层横波剩余静校正量的方法。所述方法包括:收集目标倾斜地层或各向异性地层的转换波地震道集数据、转换波叠加速度和切除参数;保留转换波地震道集中的正偏移距数据或负偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据或负偏移距数据进行叠加处理,以获得半偏移距叠加剖面;采用共接收波点叠加反射层位匹配法处理所述半偏移距叠加剖面,从而求得目标倾斜地层或各向异性地层的横波剩余静校正量。本发明尤其适合于求取倾斜地层或各向异性明显地层的横波剩余静校正量;并能够提高倾斜地层或各向异性地层的横波剩余静校正量求取的准确度。

Description

求取倾斜地层或各向异性地层横波剩余静校正量的方法
本发明涉及地球物理勘探调查技术领域,主要应用于石油、页岩气、煤炭等地震勘探的转换波地震资料处理,更具体地讲,涉及一种半偏移距叠加求取倾斜地层横波剩余静校正量的方法。
背景技术
通常,多波地震勘探利用了地震勘探过程中的多种矢量波场信息,在岩性分析、流体识别、裂缝检测、气云区成像等方面显示出其独特优势,在实际工业生产中越来越受到重视。例如,近年来,四川盆地已采集超过5000km2三维多波地震资料,勘探目标也从中浅层平缓构造拓展到深层具有陡倾地层的复杂构造。目前,陆上多波地震主要是通过炸药激发纵波,三分量检波器同时接收纵波(可简称为PP波)和转换横波(可简称为转换波,或PS波),其中,近地表横波静校正处理是影响陆上多波地震勘探数据成像的最关键的因素之一。
一般来讲,在生产实践中,横波低降速带横向变化异常剧烈,与纵波低降速带在厚度和速度两方面都具有很大差异,因此,建立近地表横波速度模型难度较大,接收点处的近地表横波静校正量求取误差也较大。总体来讲,近地表横波静校正量求取大致有三类方法。第一类是通过表层结构调查建立横波速度模型而获得,其主要问题是,野外地表条件复杂多样,施工难度大,获得的横波资料差,调查点密度低,导致横波速度模型精度低;第二类是将纵波静校正量乘以一固定的比例系数作为横波静校正量,该方法是预设一系列比例系数,分别与已知的纵波静校正量相乘作为横波静校正量,并分别应用到转换波叠加处理中,根据叠加剖面成像的好坏来选择最佳的比例系数,该方法不适用于横向变化的纵、横波速度模型结构;第三类延迟时静校正方法需要在单炮上拾取纵波、横波初至,适用于水平层状地层,拾取工作量大且横波初至不易拾取。
上述三种近静校正量已达到纵波的2-8倍,严重影响PS转换波地震资料的聚焦成像。目前,针对横波剩余静校正量大的问题,常规方法是采用共接收波点叠加反射层位匹配法求取接收点处大的横波剩余静校正量。该方法是通过分别拾取纵波、转换波共接收点全偏移距叠加剖面上同一反射层位并匹配的方法求取接收点处纵横波速度比,再根据已知的纵波静校正量和纵横波速度比求取大的横波剩余静校正量。但是,该方法在遇到陡倾地层或明显的地层各向异性时,由于转换波的共转换点在倾斜地层界面上非常分散,导致全偏移距道集中数据间的反射时差较大,使倾斜地层在叠加剖面上不能成像,因而在叠加剖面上无法拾取相应的反射层,导致该方法无法应用。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。
本发明的目的之一在于提供一种适合于求取倾斜地层或各向异性地层的横波剩余静校正量的方法。
本发明的一方面提供了一种求取倾斜地层横波剩余静校正量的方法。所述方法采用共接收波点叠加反射层位匹配法求取横波剩余静校正量,并包括以下步骤:收集目标倾斜地层的转换波地震道集数据、转换波叠加速度和切除参数;保留转换波地震道集中的正偏移距数据或负偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据或负偏移距数据进行叠加处理,以获得半偏移距叠加剖面;采用所述共接收波点叠加反射层位匹配法处理所述半偏移距叠加剖面,从而求得目标倾斜地层横波剩余静校正量。
本发明的另一方面提供了一种求取各向异性地层横波剩余静校正量的方法。所述方法采用共接收波点叠加反射层位匹配法求取横波剩余静校正量,并包括以下步骤:收集目标各向异性地层的转换波地震道集数据、转换波叠加速度和切除参数;保留转换波地震道集中的正偏移距数据或负偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据或负偏移距数据进行叠加处理,以获得半偏移距叠加剖面;采用所述共接收波点叠加反射层位匹配法处理所述半偏移距叠加剖面,从而求得目标各向异性地层横波剩余静校正量。
与现有技术相比,本发明的有益效果包括:能够适合于求取倾斜地层或各向异性地层的横波剩余静校正量;能够提高倾斜地层或各向异性地层的横波剩余静校正量求取的准确度。
具体实施方式
在下文中,将结合示例性实施例来详细说明本发明的求取倾斜地层或各向异性地层横波剩余静校正量的方法。
总地来讲,为了能够在陡倾地层转换波资料静校正处理中有效地求取横波剩余静校正量,发明人采用半偏移距叠加方法,以削弱共转换点在倾斜地层界面上分散的效应,从而保证陡倾地层(或称倾斜地层)或各向异性明显地层能够在转换波叠加剖面上成像,进而可通过常规的共接收波点叠加反射层位匹配法求得倾斜地层或各向异性明显地层的横波剩余静校正量。
常规的共接收波点叠加反射层位匹配法求取接收点处大的横波剩余静校正量方法要求在共接收波点叠加剖面上拾取对应的地震反射层位,但是地下陡倾地层在转换波全偏移距叠加剖面上无法成像,导致该方法无法应用。本发明有效地解决了该问题。
在常用的中间激发、两边接收、多次覆盖的多波地震观测系统中,转换波共转换点分布在共接收点两侧,在地震资料处理中可分为正偏移距道集和负偏移距道集,两者的集合统称为全偏移距道集。常规的共检收点叠加剖面就是将所有共检收点处全偏移距道集叠加,在横向上可形成叠加剖面。对于各向同性、水平层状地层,动校正处理后的共检收点道集中各转换点具有相同的反射时差,具有同相叠加的特点,在叠加剖面上可以看到较为清晰的地层反射,并能够方便地拾取相应的地震反射层位。
然而,发明人发现:对于地下地层,普遍存在各向异性且常常发育有复杂的倾斜地层,因此,在倾斜地层界面上,共接收点两侧正、负偏移距道集的转换波共转换点呈分散状分布,地层倾斜度越大,各向异性越明显,共转换点的分散度越高,正、负偏移距道集数据之间的反射时差也越大。如果正偏移距道集的共转换点分布在倾斜地层的上倾方向,则正偏移距道集的共转换点分布在倾斜地层的下倾方向;如果正偏移距道集的共转换点分布在倾斜地层的下倾方向,则正偏移距道集的共转换点分布在倾斜地层的上倾方向。因此,动校正处理后的共接收点道集中各转换点具有不同的反射时差,尤其是正、负偏移距道集数据之间的反射时差相对最大,导致包含正、负偏移距道集的全偏移距道集不能同相叠加或叠加效果差,在叠加剖面上无法看到较为清晰的倾斜地层反射,因而无法拾取相应的地震反射层位。
针对地层倾角和各向异性因素造成的转换波共接收点全偏移距道集中正、负偏移距数据之间存在反射时差,全偏移距叠加不能成像的问题,发明人采取仅选择全偏移距数据中正或负偏移距数据进行叠加,以消除或削弱数据间的反射时差,使倾斜地层在转换波叠加剖面上可以成像,为常规的共接收波点叠加反射层位匹配法求取接收点处大的横波剩余静校正量方法能够在陡倾地层和各向异性明显的探区应用奠定了基础。
在本发明的一个示例性实施例中,求取倾斜地层横波剩余静校正量的方法采用共接收波点叠加反射层位匹配法求取横波剩余静校正量。具体来讲,求取倾斜地层横波剩余静校正量的方法通过以下步骤来实现:
收集目标倾斜地层的转换波地震道集数据、转换波叠加速度和切除参数。这里,转换波地震道集数据可以为三维数据或二维数据。另外,为了方便计算,可利用常规计算软件平台来实现本步骤。例如,输入转换波共接收点地震道集、转换波叠加速度、切除参数文件到常规的地震数据处理软件中。常规的地震数据处理软件中都具有切除功能,数据处理人员根据业内专业知识和经验从地震道集处理过程中获得切除文件。
可仅保留转换波地震道集中的正偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据进行叠加处理,以获得半偏移距叠加剖面。另外,也可以仅保留转换波地震道集中的负偏移距数据,并利用转换波叠加速度和切除参数对所保留的负偏移距数据进行叠加处理,以获得半偏移距叠加剖面。常规的叠加处理是同时保留地震道集中的正偏移距数据和负偏移距数据进行叠加。
采用常规的共接收波点叠加反射层位匹配法处理前述半偏移距叠加剖面,从而求得目标倾斜地层横波剩余静校正量。例如,(1)拾取纵波半偏移距叠加剖面上的横向连续性好、地震反射特征明显的某一地震地质反射层位,获得反射时间Tp;(2)在转换横波共接收点半偏移距叠加剖面上拾取同一地震地质反射层位,获得反射时间Tps(由于横波速度小于纵波速度,导致同一地质层位在纵波和转换横波共接收点叠加剖面上具有不同的地震反射时间,因此,反射层位匹配是指必须对比解释纵波和转换横波共接收点半偏移距叠加剖面,才能准确拾取转换横波共接收点叠加剖面上对应的层位);(3)根据以下公式计算纵横波速度比r0
r0=(2Tps-1)/Tp
其中,r0、Tp、Tps分别是纵横波速度比、拾取层位的纵波反射时间、拾取层位的转换横波反射时间;(4)根据已知的纵波静校正量和纵横波速度比求取目标倾斜地层横波剩余静校正量,例如,求取接收点处大的横波剩余静校正量。
在本发明的另一个示例性实施例中,求取各向异性地层横波剩余静校正量的方法采用共接收波点叠加反射层位匹配法求取横波剩余静校正量。
具体来讲,求取各向异性地层横波剩余静校正量的方法通过以下步骤来实现:
收集目标各向异性地层的转换波地震道集数据、转换波叠加速度和切除参数。这里,转换波地震道集数据可以为三维数据或二维数据。另外,为了方便计算,可利用常规计算软件平台来实现本步骤。例如,输入转换波共接收点地震道集、转换波叠加速度、切除参数文件到常规的地震数据处理软件中。常规的地震数据处理软件中都具有切除功能,数据处理人员根据业内专业知识和经验从地震道集处理过程中获得切除文件。
仅保留转换波地震道集中的正偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据进行叠加处理,以获得半偏移距叠加剖面。常规的叠加处理是同时保留地震道集中的正偏移距数据和负偏移距数据进行叠加。另外,也可以仅保留转换波地震道集中的负偏移距数据,并利用转换波叠加速度和切除参数对所保留的负偏移距数据进行叠加处理,以获得半偏移距叠加剖面。
采用常规的共接收波点叠加反射层位匹配法处理前述半偏移距叠加剖面,从而求得目标各向异性地层横波剩余静校正量。例如,(1)拾取纵波半偏移距叠加剖面上的横向连续性好、地震反射特征明显的某一地震地质反射层位,获得反射时间Tp;(2)在转换横波共接收点半偏移距叠加剖面上拾取同一地震地质反射层位,获得反射时间Tps(由于横波速度小于纵波速度,导致同一地质层位在纵波和转换横波共接收点叠加剖面上具有不同的地震反射时间,因此,反射层位匹配是指必须对比解释纵波和转换横波共接收点半偏移距叠加剖面,才能准确拾取转换横波共接收点叠加剖面上对应的层位);(3)根据以下公式计算纵横波速度比r0
r0=(2Tps-1)/Tp
其中,r0、Tp、Tps分别是纵横波速度比、拾取层位的纵波反射时间、拾取层位的转换横波反射时间;(4)根据已知的纵波静校正量和纵横波速度比求取目标各向异性地层横波剩余静校正量,例如,求取接收点处大的横波剩余静校正量。
综上所述,本发明的方法针对因地层倾角或各向异性因素造成的转换波共接收点全偏移距道集中正、负偏移距数据之间存在反射时差,全偏移距叠加不能成像的问题,通过采取仅选择全偏移距数据中的正偏移距数据或负偏移距数据进行叠加,从而消除或削弱了数据间的反射时差,使倾斜地层在转换波叠加剖面上可以成像,进而通过常规的共接收波点叠加反射层位匹配法求取接收点处大的横波剩余静校正量。本发明的方法尤其适合于在陡倾地层和各向异性明显的探区应用。本发明的方法能应用于二维和三维转换波多波地震资料处理,对油气、页岩气、矿产资源勘探有重要的应用价值。
尽管上面已经结合示例性实施例描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。

Claims (4)

1.一种求取倾斜地层横波剩余静校正量的方法,所述方法采用共接收点叠加反射层位匹配法求取横波剩余静校正量,其特征在于,所述方法包括以下步骤:
收集目标倾斜地层的转换波地震道集数据、转换波叠加速度和切除参数;
保留转换波地震道集中的正偏移距数据或负偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据或负偏移距数据进行叠加处理,以获得半偏移距叠加剖面;
采用所述共接收点叠加反射层位匹配法处理所述半偏移距叠加剖面,从而求得目标倾斜地层横波剩余静校正量,
其中,所述采用共接收点叠加反射层位匹配法处理半偏移距叠加剖面的步骤包括:通过分别拾取纵波、转换横波共接收点半偏移距叠加剖面上同一反射层位并匹配的方法求取接收点处纵横波速度比,再根据已知的纵波静校正量和纵横波速度比求取目标倾斜地层横波剩余静校正量。
2.根据权利要求1所述的求取倾斜地层横波剩余静校正量的方法,其特征在于,所述转换波地震道集数据为三维数据或二维数据。
3.一种求取各向异性地层横波剩余静校正量的方法,所述方法采用共接收点叠加反射层位匹配法求取横波剩余静校正量,其特征在于,所述方法包括以下步骤:
收集目标各向异性地层的转换波地震道集数据、转换波叠加速度和切除参数;
保留转换波地震道集中的正偏移距数据或负偏移距数据,并利用转换波叠加速度和切除参数对所保留的正偏移距数据或负偏移距数据进行叠加处理,以获得半偏移距叠加剖面;
采用所述共接收点叠加反射层位匹配法处理所述半偏移距叠加剖面,从而求得目标各向异性地层横波剩余静校正量,
其中,所述采用共接收点叠加反射层位匹配法处理半偏移距叠加剖面的步骤包括:通过分别拾取纵波、转换横波共接收点半偏移距叠加剖面上同一反射层位并匹配的方法求取接收点处纵横波速度比,再根据已知的纵波静校正量和纵横波速度比求取目标各向异性地层横波剩余静校正量。
4.根据权利要求3所述的求取各向异性地层横波剩余静校正量的方法,其特征在于,所述转换波地震道集数据为三维数据或二维数据。
CN201611038311.2A 2016-11-11 2016-11-11 求取倾斜地层或各向异性地层横波剩余静校正量的方法 Active CN106443791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611038311.2A CN106443791B (zh) 2016-11-11 2016-11-11 求取倾斜地层或各向异性地层横波剩余静校正量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611038311.2A CN106443791B (zh) 2016-11-11 2016-11-11 求取倾斜地层或各向异性地层横波剩余静校正量的方法

Publications (2)

Publication Number Publication Date
CN106443791A CN106443791A (zh) 2017-02-22
CN106443791B true CN106443791B (zh) 2018-07-10

Family

ID=58220364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611038311.2A Active CN106443791B (zh) 2016-11-11 2016-11-11 求取倾斜地层或各向异性地层横波剩余静校正量的方法

Country Status (1)

Country Link
CN (1) CN106443791B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112305613B (zh) * 2019-07-25 2024-03-01 中国石油天然气集团有限公司 转换横波检波点静校正方法及装置
CN110471112A (zh) * 2019-09-20 2019-11-19 西南石油大学 基于叠加速度变化的倾斜层层速度异常反演方法
CN112946745B (zh) * 2019-12-11 2024-06-25 中国石油天然气集团有限公司 一种基于引力拟合的转换波剩余静校正方法及系统
CN113514888B (zh) * 2021-06-18 2024-04-16 中国石油化工股份有限公司 一种基于非剩余时差的多次波压制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5663928A (en) * 1994-10-19 1997-09-02 Elf Aquitaine Production Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN103487834A (zh) * 2013-09-09 2014-01-01 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波共检波点叠加静校正方法
CN103630934A (zh) * 2012-08-23 2014-03-12 中国石油天然气集团公司 一种确定转换波检波点大的横波静校正量的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5663928A (en) * 1994-10-19 1997-09-02 Elf Aquitaine Production Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法
CN103630934A (zh) * 2012-08-23 2014-03-12 中国石油天然气集团公司 一种确定转换波检波点大的横波静校正量的方法
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN103487834A (zh) * 2013-09-09 2014-01-01 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波共检波点叠加静校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于偏移成像道集的剩余静校正方法;徐锦承 等;《地球物理学报》;20160229;第59卷(第2期);第655-663页 *

Also Published As

Publication number Publication date
CN106443791A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN107526101B (zh) 一种获取地震反射波的采集和处理方法
CN104570125B (zh) 一种利用井数据提高成像速度模型精度的方法
CN106094032B (zh) 一种构建地层速度模型的方法
Urosevic et al. Targeting nickel sulfide deposits from 3D seismicreflection data at Kambalda, Australia
AU2014227491B2 (en) Construction and application of angle gathers from three-dimensional imaging of multiples wavefields
CN104678434B (zh) 一种预测储层裂缝发育参数的方法
CN106443791B (zh) 求取倾斜地层或各向异性地层横波剩余静校正量的方法
US20080137480A1 (en) Method of Building a Subsurface Velocity Model
CN102053261A (zh) 一种地震数据处理方法
CN111123359B (zh) 随钻测井与地层格架约束的井周地震成像探测方法及装置
CN105093320A (zh) 针对高速结晶盐壳覆盖区层析静校正初至拾取方法
CN105425299B (zh) 确定地层裂缝分布的方法和装置
CN104483705A (zh) 一种三维剩余静校正方法
Ahmadi et al. High-resolution 2D seismic imaging and forward modeling of a polymetallic sulfide deposit at Garpenberg, central Sweden
CN110687602A (zh) 浅层地震多波联合勘探方法
Brodic et al. Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging
Takougang et al. Extraction and characterization of faults and fractures from 3D VSP data in a carbonate reservoir: A workflow
Zhang et al. Simultaneous inversion for microseismic event location and velocity model in Vaca Muerta Formation
CN107656308B (zh) 一种基于时间深度扫描的共散射点叠前时间偏移成像方法
Landa et al. Point and edge diffractions in three dimensions
CN104422955B (zh) 一种利用旅行时变化量进行各向异性参数提取的方法
CN102053260A (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN102053262A (zh) 获得地震转换波的方位速度的方法及处理地震数据的方法
US20140257707A1 (en) Progressive 3D Vertical Seismic Profiling Method
CN109143398B (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
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: CNPC CHUANQING DRILLING ENGINEERING Co.,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: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

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

Effective date of registration: 20180330

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

Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

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

Applicant before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

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

Effective date of registration: 20201111

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.