CN108897040B - 一种确定陆检与海底之间的耦合性的方法及装置 - Google Patents

一种确定陆检与海底之间的耦合性的方法及装置 Download PDF

Info

Publication number
CN108897040B
CN108897040B CN201810440808.XA CN201810440808A CN108897040B CN 108897040 B CN108897040 B CN 108897040B CN 201810440808 A CN201810440808 A CN 201810440808A CN 108897040 B CN108897040 B CN 108897040B
Authority
CN
China
Prior art keywords
seismic data
detection
seismic
land
root
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
CN201810440808.XA
Other languages
English (en)
Other versions
CN108897040A (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
China National Petroleum Corp
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 China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201810440808.XA priority Critical patent/CN108897040B/zh
Publication of CN108897040A publication Critical patent/CN108897040A/zh
Application granted granted Critical
Publication of CN108897040B publication Critical patent/CN108897040B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例公开了一种确定陆检与海底之间的耦合性的方法及装置。所述方法包括:获取水检共检波点道集数据/陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;分别确定第一水检地震数据、第二水检地震数据、第一陆检地震数据和第二陆检地震数据的地震振幅属性值;基于第一水检地震数据、第二水检地震数据、第一陆检地震数据和第二陆检地震数据的地震振幅属性值,确定目标海域的陆检与海底之间的耦合等级。本申请实施例提供的技术方案,可以简化评价陆检与海底之间的耦合性的过程,满足野外采集现场质量监控需求。

Description

一种确定陆检与海底之间的耦合性的方法及装置
技术领域
本申请涉及海底地震勘探技术领域,特别涉及一种确定陆检与海底之间的耦合性的方法及装置。
背景技术
海底地震勘探是指将压电检波器(水检)和速度检波器(陆检)一起放置在海底进行地震数据采集的技术,具体包括海底电缆地震采集和海底节点采集两种方式。水检由地层反射引起周围海水压力变化来产生地震响应,与海底之间不存在耦合问题;陆检是通过感知海底表面质点的振动来记录地震信号,其与大地组成共振耦合系统。由于受洋流、海底地形及施工方法等方面的影响,陆检与海底构成的共振耦合系统的耦合性可能较差,这样,该共振耦合系统可能会严重改造陆检信号,即改变陆检接收地震信号的振幅、频率及相位特征,从而降低地震信号的信噪比。因此,陆检与海底之间的耦合性是评价地震资料品质的重要手段之一,也是野外地震施工的重点工作。
目前常规确定陆检与海底之间耦合性的方法主要是利用“单自由度有阻尼的自由振动系统”来描述陆检与海底之间的耦合效应,并通过求取陆检与海底构成的耦合系统的传输函数来评价该耦合效应。由于该方法需要求取耦合系统的传输函数来评价耦合效应,导致评价过程过于繁琐,不满足野外采集现场质量监控需求。
发明内容
本申请实施例的目的是提供一种确定陆检与海底之间的耦合性的方法及装置,以简化评价陆检与海底之间的耦合性的过程,满足野外采集现场质量监控需求。
为解决上述技术问题,本申请实施例提供一种确定陆检与海底之间的耦合性的方法及装置是这样实现的:
一种确定陆检与海底之间的耦合性的方法,提供有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据;所述方法包括:
获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
优选方案中,分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,包括:
分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值;
将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
优选方案中,采用下述公式计算所述第一水检地震数据的均方根振幅值:
Figure BDA0001655821050000021
其中,κhy表示所述第一水检地震数据的均方根振幅值,HyPi,j表示所述第一指定时窗内的第i个地震道中第j个采样点的振幅值,i=1,2,3,…,m1,m1表示所述第一指定时窗内的地震道的数量,j=1,2,3,…,n1,n1表示所述第一指定时窗内的一个地震道中的采样点的数量。
优选方案中,基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级,包括:
将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值;
根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
优选方案中,根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级,包括:
将所述第一均方根比值减去所述第二均方根比值,得到差异值;
根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级。
优选方案中,根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级,包括:
当所述差异值大于或等于所述预设差异阈值时,确定所述目标海域的陆检与海底之间的耦合等级为一级;
当所述差异值小于所述预设差异阈值时,确定所述目标海域的陆检与海底之间的耦合等级为二级。
一种确定陆检与海底之间的耦合性的装置,所述装置提供目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据;所述装置包括:地震数据获取模块、振幅属性值确定模块和耦合等级确定模块;其中,
所述地震数据获取模块,用于获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
所述振幅属性值确定模块,用于分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
所述耦合等级确定模块,用于基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
优选方案中,所述振幅属性值确定模块用于分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值,并将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
优选方案中,所述耦合等级确定模块用于将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值,并根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
一种确定陆检与海底之间的耦合性的装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述存储器中存储有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据,所述计算机程序被所述处理器运行时执行以下步骤:
获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
由以上本申请实施例提供的技术方案可见,本申请实施例提供的确定陆检与海底之间的耦合性的方法及装置,可以获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;可以分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;可以基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。如此,采用本申请方法无需计算陆检与海底构成的共振耦合系统的传输函数,便可以确定陆检与海底之间的耦合性,从而可以简化评价陆检与海底之间的耦合性的过程,满足野外采集现场质量监控需求。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种确定陆检与海底之间的耦合性的方法实施例的流程图;
图2是本申请实施例中水检共检波点道集数据、第一水检地震数据和第二水检地震数据的剖面示意图;
图3是本申请实施例中陆检共检波点道集数据、第一陆检地震数据和第二陆检地震数据的剖面示意图;
图4是本申请实施例中不同检波点位置处对应第一均方根比值和第二均方根比值的对比示意图;
图5是本申请确定陆检与海底之间的耦合性的装置的一种实施例的组成结构示意图;
图6是本申请确定陆检与海底之间的耦合性的装置的另一种实施例的组成结构示意图。
具体实施方式
本申请实施例提供一种确定陆检与海底之间的耦合性的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
本申请实施例提供了一种确定陆检与海底之间的耦合性的方法。所述确定陆检与海底之间的耦合性的方法可以提供有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据。
在本实施方式中,所述目标海域可以是指已进行海底地震勘探的海域。
在本实施方式中,可以通过海底电缆地震采集或海底节点采集的方式,获取所述目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据。其中,所述指定检波点位置可以用于表征海底排布的检波线上的任一检波点位置。所述水检共检波点道集数据可以通过在所述指定检波点位置处的水检,例如压电检波器,采集所述目标海域中各个炮点发射并经过地层反射的地震波后得到的。所述陆检共检波点道集数据可以通过在所述指定检波点位置处的陆检,例如速度检波器,采集所述目标海域中各个炮点发射并经过地层反射的地震波后得到的。
图1是本申请一种确定陆检与海底之间的耦合性的方法实施例的流程图。如图1所示,所述确定陆检与海底之间的耦合性的方法,包括以下步骤。
步骤S101:获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波。
在本实施方式中,可以获取所述水检共检波点道集数据中第一指定时窗内的第一水检地震数据和第二指定时窗内的第二水检地震数据。还可以获取所述陆检共检波点道集数据中所述第一指定时窗内的第一陆检地震数据和所述第二指定时窗内的第二陆检地震数据。其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间。所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间。所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波。
由于检波器采集的地震数据中采样时间小于初始时间的部分地震波信号通常是背景噪音,而从炮点发射并经过地层仅一次反射后被检波器采集的地震波信号通常是有效信号或强反射信号,这样,所述水检共检波点道集数据中第一指定时窗内的第一水检地震数据和所述陆检共检波点道集数据中第一指定时窗内的第一陆检地震数据分别用于表征所述水检共检波点道集数据中背景噪音时窗内的水检背景噪音信号和所述陆检共检波点道集数据中背景噪音时窗内的陆检背景噪音信号。所述水检共检波点道集数据中第二指定时窗内的第二水检地震数据和所述陆检共检波点道集数据中第二指定时窗内的第二陆检地震数据分别用于表征所述水检共检波点道集数据中有效信号时窗内的水检有效信信号和所述陆检共检波点道集数据中有效信号时窗内的陆检有效信号。
例如,图2是本申请实施例中水检共检波点道集数据、第一水检地震数据和第二水检地震数据的剖面示意图。图3是本申请实施例中陆检共检波点道集数据、第一陆检地震数据和第二陆检地震数据的剖面示意图。图2和图3中的横坐标为炮点桩号和检波点桩号,纵坐标为采样时间,单位为毫秒(ms)。其中,图2和图3中的time表示采样时间。在本实施例中,共检波点桩号为5776。所述第一指定时窗内的地震道为道号39至62的地震道,共23道,采样时间的范围为100ms至500ms,采样时间间隔为2ms,共201个采样点。所述第一指定时窗内的地震道为道号66至89的地震道,共23道,采样时间的范围为1100ms至1500ms,采样时间间隔为2ms,共201个采样点。
步骤S102:分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征。
在本实施方式中,分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,具体可以包括,可以分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值。可以将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
在本实施方式中,可以采用下述公式计算所述第一水检地震数据的均方根振幅值:
Figure BDA0001655821050000071
其中,κhy表示所述第一水检地震数据的均方根振幅值,HyPi,j表示所述第一指定时窗内的第i个地震道中第j个采样点的振幅值,i=1,2,3,…,m1,m1表示所述第一指定时窗内的地震道的数量,j=1,2,3,…,n1,n1表示所述第一指定时窗内的一个地震道中的采样点的数量。
可以采用下述公式计算所述第一陆检地震数据的均方根振幅值:
Figure BDA0001655821050000072
其中,κge表示所述第一陆检地震数据的均方根振幅值,GePi,j表示所述第一指定时窗内的第i个地震道中第j个采样点的振幅值,i=1,2,3,…,m1,m1表示所述第一指定时窗内的地震道的数量,j=1,2,3,…,n1,n1表示所述第一指定时窗内的一个地震道中的采样点的数量。
可以采用下述公式计算所述第二水检地震数据的均方根振幅值:
Figure BDA0001655821050000081
其中,χhy表示所述第二水检地震数据的均方根振幅值,HyNk,v表示所述第二指定时窗内的第k个地震道中第v个采样点的振幅值,k=1,2,3,…,m2,m2表示所述第二指定时窗内的地震道的数量,v=1,2,3,…,n2,n2表示所述第二指定时窗内的一个地震道中的采样点的数量。
可以采用下述公式计算所述第二陆检地震数据的均方根振幅值:
Figure BDA0001655821050000082
其中,χhy表示所述第二陆检地震数据的均方根振幅值,GeNk,v表示所述第二指定时窗内的第k个地震道中第v个采样点的振幅值,k=1,2,3,…,m2,m2表示所述第二指定时窗内的地震道的数量,v=1,2,3,…,n2,n2表示所述第二指定时窗内的一个地震道中的采样点的数量。
步骤S103:基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
在本实施方式中,基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级,具体可以包括,可以将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,例如,采用α表示所述第一均方根比值,则α=κhyge,其中,κhy表示所述第一水检地震数据的均方根振幅值,κge表示所述第一陆检地震数据的均方根振幅值。还可以将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值,例如,采用β表示所述第二均方根比值,则β=χhyge,其中,χhy表示所述第二水检地震数据的均方根振幅值,χge表示所述第二陆检地震数据的均方根振幅值。可以根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
在本实施方式中,根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级,具体可以包括,可以将所述第一均方根比值减去所述第二均方根比值,得到差异值。可以根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级。
在本实施方式中,根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级,具体可以包括,当所述差异值大于或等于所述预设差异阈值时,可以确定所述目标海域的陆检与海底之间的耦合等级为一级。或者,当所述差异值小于所述预设差异阈值时,可以确定所述目标海域的陆检与海底之间的耦合等级为二级。其中,一级对应的耦合性比二级对应的耦合性好。通常,当所述目标海域的陆检与海底之间的耦合等级为一级时,表明所述目标海域的陆检与海底之间的耦合效果较好,当所述目标海域的陆检与海底之间的耦合等级为二级时,表明所述目标海域的陆检与海底之间的耦合效果较差。
在本实施方式中,所述预设差异阈值的取值范围可以包括50~150。在实际应用过程中,由于所述预设差异阈值的取值通常与上述有效信号和背景噪音信号的时窗窗长、地震勘探施工环境、检波器类型、海底介质等多方面因素相关联,因此,所述预设差异阈值的具体取值可以根据实际情况而定。
例如,图4是本申请实施例中不同检波点位置处对应第一均方根比值和第二均方根比值的对比示意图。其中,图4中的横坐标和纵坐标分别为检波点桩号和均方根比值,图4中的虚线和实线分别表示第一均方根比值和第二均方根比值。图4中的检波点桩号5776对应的均方根比值是基于图2和图3中的选取的指定时窗内的地震数据计算得到的。针对共检波点桩号5776,所述第一指定时窗内的地震道为道号39至62的地震道,共23道,即m1=23;采样时间的范围为100ms至500ms,采样时间间隔为2ms,共201个采样点,即n1=201。所述第一指定时窗内的地震道为道号66至89的地震道,共23道,即m2=23;采样时间的范围为1100ms至1500ms,采样时间间隔为2ms,共201个采样点,即n2=201。这样,可以计算得到共检波点桩号5776位置对应的第一水检地震数据的均方根振幅值
Figure BDA0001655821050000101
可以计算得到共检波点桩号5776位置对应的第一陆检地震数据的均方根振幅值
Figure BDA0001655821050000102
可以计算得到共检波点桩号5776位置对应的第二水检地震数据的均方根振幅值
Figure BDA0001655821050000103
可以计算得到共检波点桩号5776位置对应的第二陆检地震数据的均方根振幅值
Figure BDA0001655821050000104
那么,可以计算得到共检波点桩号5776位置对应的第一均方根比值α=κhyge=188.4235,以及可以计算得到共检波点桩号5776位置对应的第二均方根比值β=χhyge=1.229。
所述确定陆检与海底之间的耦合性的方法实施例,可以获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;可以分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;可以基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。如此,采用本申请方法无需计算陆检与海底构成的共振耦合系统的传输函数,便可以确定陆检与海底之间的耦合性,从而可以简化评价陆检与海底之间的耦合性的过程,满足野外采集现场质量监控需求。
图5是本申请确定陆检与海底之间的耦合性的装置的一种实施例的组成结构示意图。所述确定陆检与海底之间的耦合性的装置可以提供目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据。如图5所示,所述确定陆检与海底之间的耦合性的装置可以包括:地震数据获取模块100、振幅属性值确定模块200和耦合等级确定模块300。
所述地震数据获取模块100,可以用于获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波。
所述振幅属性值确定模块200,可以用于分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征。
所述耦合等级确定模块300,可以用于基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
在本实施方式中,所述振幅属性值确定模块200可以用于分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值,并将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
在本实施方式中,所述耦合等级确定模块300可以用于将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值,并根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
图6是本申请确定陆检与海底之间的耦合性的装置的另一种实施例的组成结构示意图。如图6所示,所述确定陆检与海底之间的耦合性的装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述存储器中存储有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据,所述计算机程序被所述处理器运行时执行以下步骤:
步骤S101:获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
步骤S102:分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
步骤S103:基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级。
所述确定陆检与海底之间的耦合性的装置实施例与所述确定陆检与海底之间的耦合性的方法实施例相对应,可以实现确定陆检与海底之间的耦合性的方法实施例的技术方案,并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的装置、模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (8)

1.一种确定陆检与海底之间的耦合性的方法,其特征在于,提供有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据;所述方法包括:
获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级;其中,包括:将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值;根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
2.根据权利要求1所述的方法,其特征在于,分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,包括:
分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值;
将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
3.根据权利要求2所述的方法,其特征在于,采用下述公式计算所述第一水检地震数据的均方根振幅值:
Figure FDA0002358061410000021
其中,κhy表示所述第一水检地震数据的均方根振幅值,HyPi,j表示所述第一指定时窗内的第i个地震道中第j个采样点的振幅值,i=1,2,3,…,m1,m1表示所述第一指定时窗内的地震道的数量,j=1,2,3,…,n1,n1表示所述第一指定时窗内的一个地震道中的采样点的数量。
4.根据权利要求1所述的方法,其特征在于,根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级,包括:
将所述第一均方根比值减去所述第二均方根比值,得到差异值;
根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级。
5.根据权利要求4所述的方法,其特征在于,根据预设差异阈值和所述差异值,确定所述目标海域的陆检与海底之间的耦合等级,包括:
当所述差异值大于或等于所述预设差异阈值时,确定所述目标海域的陆检与海底之间的耦合等级为一级;
当所述差异值小于所述预设差异阈值时,确定所述目标海域的陆检与海底之间的耦合等级为二级。
6.一种确定陆检与海底之间的耦合性的装置,其特征在于,所述装置提供目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据;所述装置包括:地震数据获取模块、振幅属性值确定模块和耦合等级确定模块;其中,
所述地震数据获取模块,用于获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
所述振幅属性值确定模块,用于分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
所述耦合等级确定模块,用于基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级;其中,包括:将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值;根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
7.根据权利要求6所述的装置,其特征在于,所述振幅属性值确定模块用于分别计算所述第一水检地震数据的均方根振幅值、所述第二水检地震数据的均方根振幅值、所述第一陆检地震数据的均方根振幅值和所述第二陆检地震数据的均方根振幅值,并将所述第一水检地震数据的均方根振幅值作为所述第一水检地震数据的地震振幅属性值、将所述第二水检地震数据的均方根振幅值作为所述第二水检地震数据的地震振幅属性值、将所述第一陆检地震数据的均方根振幅值作为所述第一陆检地震数据的地震振幅属性值,以及将所述第二陆检地震数据的均方根振幅值作为所述第二陆检地震数据的地震振幅属性值。
8.一种确定陆检与海底之间的耦合性的装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,其特征在于,所述存储器中存储有目标海域中指定检波点位置处对应的水检共检波点道集数据和陆检共检波点道集数据,所述计算机程序被所述处理器运行时执行以下步骤:
获取所述水检共检波点道集数据/所述陆检共检波点道集数据中第一指定时窗内的第一水检地震数据/第一陆检地震数据,以及第二指定时窗内的第二水检地震数据/第二陆检地震数据;其中,所述第一指定时窗内的采样时间小于所述第一指定时窗内的地震道对应的初至时间;所述第二指定时窗内的采样时间表示检波器采集到指定地震波的时间;所述指定地震波表示从所述目标海域的炮点发射并经过地层一次反射后的地震波;
分别确定所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值;其中,所述地震振幅属性值用于表征地震数据的振幅特征;
基于所述第一水检地震数据、所述第二水检地震数据、所述第一陆检地震数据和所述第二陆检地震数据的地震振幅属性值,确定所述目标海域的陆检与海底之间的耦合等级;其中,包括:将所述第一水检地震数据的均方根振幅值除以所述第一陆检地震数据的均方根振幅值,得到第一均方根比值,以及将所述第二水检地震数据的均方根振幅值除以所述第二陆检地震数据的均方根振幅值,得到第二均方根比值;根据所述第一均方根比值和所述第二均方根比值,确定所述目标海域的陆检与海底之间的耦合等级。
CN201810440808.XA 2018-05-10 2018-05-10 一种确定陆检与海底之间的耦合性的方法及装置 Active CN108897040B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810440808.XA CN108897040B (zh) 2018-05-10 2018-05-10 一种确定陆检与海底之间的耦合性的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810440808.XA CN108897040B (zh) 2018-05-10 2018-05-10 一种确定陆检与海底之间的耦合性的方法及装置

Publications (2)

Publication Number Publication Date
CN108897040A CN108897040A (zh) 2018-11-27
CN108897040B true CN108897040B (zh) 2020-05-08

Family

ID=64342709

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810440808.XA Active CN108897040B (zh) 2018-05-10 2018-05-10 一种确定陆检与海底之间的耦合性的方法及装置

Country Status (1)

Country Link
CN (1) CN108897040B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147716A (zh) * 2020-09-15 2020-12-29 中海油田服务股份有限公司 一种水陆检信号的合并方法和装置
CN112394415B (zh) * 2020-10-21 2024-06-25 中国石油天然气集团有限公司 海底节点水平方位角的确定方法及装置
CN112327360B (zh) * 2020-10-29 2022-12-09 中海油田服务股份有限公司 一种实现噪声处理的方法、装置、计算机存储介质及终端
CN112415595B (zh) * 2020-10-30 2024-05-28 中国石油天然气集团有限公司 海底节点二次定位的方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104865601A (zh) * 2015-05-18 2015-08-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 确定地震勘探中采集的地震记录的类型的方法
CN105223610A (zh) * 2015-10-22 2016-01-06 中国船舶重工集团公司第七六○研究所 判断海底地震仪与航行舰船地震波信号耦合效果的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107576986B (zh) * 2017-08-30 2019-03-01 中国石油集团东方地球物理勘探有限责任公司 一种确定海底地震反射波旅行时的方法和装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104865601A (zh) * 2015-05-18 2015-08-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 确定地震勘探中采集的地震记录的类型的方法
CN105223610A (zh) * 2015-10-22 2016-01-06 中国船舶重工集团公司第七六○研究所 判断海底地震仪与航行舰船地震波信号耦合效果的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Geophone-seabed coupling effect and its corretion;Zhang Bao-Qing et.al;《APPLIED GEOPHYSICS》;20160331;第13卷(第1期);145-155 *
提高双检检波器陆检信噪比的研究;韩学义 等;《物探装备》;20130831;第23卷(第4期);226-231 *

Also Published As

Publication number Publication date
CN108897040A (zh) 2018-11-27

Similar Documents

Publication Publication Date Title
CN108897040B (zh) 一种确定陆检与海底之间的耦合性的方法及装置
CN107193045B (zh) 一种地震数据处理方法及装置
CN108279436B (zh) 一种拾取地震波初至的方法及装置
CN108252710B (zh) 一种确定油藏含水饱和度的方法及装置
CN107679358B (zh) 一种确定储层渗透率的方法及装置
CN108196304B (zh) 一种多次波压制方法及装置
CN108363096B (zh) 一种水陆检波器地震数据上下行波场分离方法及装置
CN106054251B (zh) 一种初至波拾取方法及装置
CN106646611B (zh) 一种纵波和转换波地震数据时间匹配的方法及装置
CN106842304B (zh) 一种叠前深度偏移方法及装置
CN107092031B (zh) 一种声波测井数据校正方法及装置
CN107367762B (zh) 一种确定储层参数的方法及装置
CN107589455B (zh) 一种水陆检波器地震数据合并处理的方法及装置
US10921477B2 (en) Method, apparatus, and system for eliminating frequency dispersion effect
US6985404B2 (en) Method and apparatus of performing amplitude variation with offset analysis
CN108121018B (zh) 检波点定位准确度评价方法和装置
CN107144878B (zh) 一种断层识别方法及装置
CN108386185B (zh) 一种确定砂体连通性的方法及装置
CN108035708B (zh) 一种去除地层界面反射波的方法及装置
US10585198B2 (en) Noise attenuation of multicomponent microseismic data
CN109143375B (zh) 一种确定钻井方位的方法及装置
CN106772599B (zh) 一种计算地层横波速度的方法及装置
CN111239832A (zh) 3d3c-vsp成像处理方法、装置及设备
CN106772596B (zh) 一种确定叠前时间偏移速度场的方法及装置
CN111665558B (zh) 一种油气分布检测方法和系统

Legal Events

Date Code Title Description
PB01 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