CN107179546A - 一种确定速度场的方法及装置 - Google Patents
一种确定速度场的方法及装置 Download PDFInfo
- Publication number
- CN107179546A CN107179546A CN201710328456.4A CN201710328456A CN107179546A CN 107179546 A CN107179546 A CN 107179546A CN 201710328456 A CN201710328456 A CN 201710328456A CN 107179546 A CN107179546 A CN 107179546A
- Authority
- CN
- China
- Prior art keywords
- region
- offset distance
- velocity field
- velocity
- data
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 230000002159 abnormal effect Effects 0.000 claims abstract description 61
- 230000000977 initiatory effect Effects 0.000 claims abstract description 32
- 238000005070 sampling Methods 0.000 claims description 36
- 238000012545 processing Methods 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 14
- 238000012937 correction Methods 0.000 description 13
- 238000003860 storage Methods 0.000 description 10
- 241000209094 Oryza Species 0.000 description 8
- 235000007164 Oryza sativa Nutrition 0.000 description 8
- 235000009566 rice Nutrition 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 230000002547 anomalous effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- IXSZQYVWNJNRAL-UHFFFAOYSA-N etoxazole Chemical compound CCOC1=CC(C(C)(C)C)=CC=C1C1N=C(C=2C(=CC=CC=2F)F)OC1 IXSZQYVWNJNRAL-UHFFFAOYSA-N 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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/303—Analysis for determining velocity profiles or travel times
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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
技术领域
本申请涉及石油地球物理勘探技术领域,特别涉及一种确定速度场的方法及装置。
背景技术
速度分析处理是油气地震勘探资料处理过程中的一个重要环节。通过对目的工区的地震数据中共中心点(Common Middle Point,CMP)道集数据进行速度分析处理,可以得到目的工区的速度场。进而,可以为后续目的工区的地震数据的反演处理奠定基础。
现有技术中通常采用的确定目的工区的速度场的方法的主要过程是:在目的工区的地层为均匀水平层状介质的条件下,也就是在目的工区中不存在异常地质体的条件下,获取目的工区的地震数据中的一个CMP道集数据,基于预设的一组初始均方根速度,对该CMP道集数据进行动校正处理;其中,根据预设的一组初始均方根速度中一个初始均方根速度可以进行一次动校正处理,根据每一次动校正处理后的CMP道集数据,可以得到一组动校正处理后的CMP道集数据;计算每一次动校正处理后的CMP道集数据中各偏移距对应同一采样时间处的振幅之间的相似系数;这样,一个初始均方根速度可以对应得到一个相似系数,一组初始均方根速度可以对应得到多个相似系数;可以将计算得到的多个相似系数中的最大相似系数对应的初始均方根速度作为该CMP道集数据中零偏移距对应该采样时间处的目标均方根速度。重复上述操作,可以得到该CMP道集数据中零偏移距对应所有采样时间处的目标均方根速度;利用迪克斯(Dix)公式可以将目标均方根速度转换为所有采样时间对应的层速度。可以将零偏移距对应所有采样时间处的层速度进行时深转换处理,可以得到零偏移距对应不同地层深度处的层速度。由于目的工区的地层为均匀水平层状介质,该介质中不同地层位置处的层速度均相同,从而,根据零偏移距对应不同深度的层速度,可以确定目的工区的层速度场。
发明人发现现有技术中至少存在如下问题:对于结构复杂的目的工区,例如断层大量发育的盆缘地区、裂缝发育的碳酸盐地区等,其某地层深度处的地层可能为非均匀层状介质,即可能存在异常地质体,该地层中异常地质体对应的地层位置处的层速度与该地层中其他地层位置处的层速度可能不同,但现有技术所得到的层速度场中该地层的所有地层位置处的层速度却均为该CMP道集数据对应的地震射线穿过非均匀层状介质中该地层深度的所有地层位置处的平均层速度,并不能准确反映目的工区中该地层的不同地层位置处的实际层速度。
发明内容
本申请实施例的目的是提供一种确定速度场的方法及装置,以提高所确定的速度场的准确度。
为解决上述技术问题,本申请实施例提供一种确定速度场的方法及装置是这样实现的:
一种确定速度场的方法,包括:
获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,确定所述第一剖面区域的初始层速度场;
基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域;
基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据;所述第一种偏移距地震数据是与所述异常地质区域对应的地震数据;所述第二种偏移距地震数据是与所述第一剖面区域中非异常地质区域对应的地震数据;
根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场;
根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,确定所述第一剖面区域的目标层速度场。
优选方案中,所述基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域,包括:
所述第一剖面区域中的异常地质区域包括:初始层速度大于所述预设层速度阈值第一地层位置;所述第一地层位置是所述第一剖面区域中第一层段处的一个地层位置。
优选方案中,所述预设层速度阈值是所述初始层速度场中第一层段处的平均层速度。
优选方案中,所述根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,包括:
根据预设的一组均方根速度和所述第一种偏移距地震数据中第一共中心点分道集数据,确定所述第一共中心点分道集数据中各个采样时间处的目标均方根速度;
将所述目标均方根速度转换为所述第一共中心点分道集数据中各个采样时间处的目标层速度;对所述第一共中心点分道集数据中各个采样时间处的目标层速度进行时深转换处理,得到所述第一共中心点分道集数据中各个地层深度处的目标层速度;
根据所述第一共中心点分道集数据中各个地层深度处的目标层速度,确定所述第一偏移距区域的目标层速度场。
优选方案中,所述预设的一组均方根速度包括:步长为1米/秒的多个均方根速度;所述步长为1米/秒的多个均方根速度的速度范围为10~20米/秒。
优选方案中,所述第一偏移距区域内同一地层深度的不同地层位置处的目标层速度相同。
优选方案中,所述根据所述第一偏移距区域的目标层速度场和所述第二偏移距区域的目标层速度场,确定所述第一剖面区域的目标层速度场,包括:
根据所述第二偏移距区域的目标层速度场,确定第一剖面剩余区域的目标层速度场;所述第一剖面剩余区域表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域;
将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
优选方案中,所述第一剖面剩余区域与所述第二偏移距区域在同一地层深度的不同地层位置处的目标层速度相同。
一种确定速度场的装置,所述装置包括:初始层速度场确定模块、异常地质区域确定模块、偏移距地震数据划分模块、偏移距区域速度确定模块和目标层速度确定模块;其中,
所述初始层速度场确定模块,用于获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,确定所述第一剖面区域的初始层速度场;
所述异常地质区域确定模块,用于基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域;
所述偏移距地震数据划分模块,用于基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据;所述第一种偏移距地震数据是与所述异常地质区域对应的地震数据;所述第二种偏移距地震数据是与所述第一剖面区域中非异常地质区域对应的地震数据;
所述偏移距区域速度确定模块,用于根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场;
所述目标层速度确定模块,用于根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,确定所述第一剖面区域的目标层速度场。
优选方案中,所述目标层速度确定模块,包括:剩余区域速度确定模块和剖面区域速度确定模块;其中,
所述剩余区域速度确定模块,用于根据所述第二偏移距区域的目标层速度场,确定第一剖面剩余区域的目标层速度场;所述第一剖面剩余区域表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域;
所述剖面区域速度确定模块,用于将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
本申请实施例提供了一种确定速度场的方法及装置,首先,基于预设层速度阈值和所述目的工区中第一剖面区域初始层速度场,可以确定所述第一剖面区域中的异常地质区域;基于所述异常地质区域,可以将所述第一剖面区域的第一共中心点道集数据划分为第一偏移距地震数据和第二偏移距地震数据;所述第一偏移距地震数据可以是与所述异常地质区域对应的地震数据;所述第二偏移距地震数据可以是与所述第一剖面区域中非异常地质区域对应的地震数据;然后,根据所述第一偏移距地震数据可以确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二偏移距地震数据可以确定所述第一剖面区域中第二偏移距区域的目标层速度场;最后,根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,可以确定所述第一剖面区域的目标层速度场。如此,对于复杂结构的目的工区,其异常地质区域的层速度是基于穿过异常地质区域的地震射线对应的地震道数据所确定的,非异常地质区域的层速度是基于穿过非异常地质区域的地震射线对应的地震道数据所确定的,均不是基于穿过异常地质区域和非异常地质区域的地震射线对应的地震道数据所确定的平均层速度,从而可以提高最终确定的目的工区的层速度场的准确度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种确定速度场的方法实施例的流程图;
图2是本申请实施例目的工区中某一剖面区域的初始层速度场的剖面示意图;
图3是本申请实施例第一剖面区域的第一共中心点道集数据对应的偏移距划分示意图;
图4是本申请实施例中基于第一共中心点道集数据确定的第一剖面区域的目标层速度场的剖面示意图;
图5是本申请实施例中基于第一剖面区域的各个共中心点道集数据确定的第一剖面区域的目标层速度场的剖面示意图;
图6是本申请确定速度场的装置实施例的组成结构图;
图7是本申请确定速度场的装置实施例中目标层速度确定模块的组成结构图。
具体实施方式
本申请实施例提供一种确定速度场的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请一种确定速度场的方法实施例的流程图。如图1所示,所述确定速度场的方法包括以下步骤。
步骤S101:获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,确定所述第一剖面区域的初始层速度场。
所述目的工区可以是由主测线维度、联络测线维度和地层深度维度构成的三维区域。所述三个维度两两正交。所述目的工区的某一段地层深度处可以存在异常地质体。该异常地质体位置处的层速度可能明显高于该段地层深度的非异常地质体位置处的层速度。
所述第一剖面区域可以是沿主测线维度或联络测线维度上,且与地表垂直的剖面区域。
可以对所述目的工区进行地震勘探和数据采集,获取所述目的工区的地震数据。所述地震数据可以包括:所述第一剖面区域的多个共中心点(Common Middle Point,CMP)道集数据。所述CMP道集数据可以包括:至少一个地震道数据和所述地震道数据对应的偏移距。所述地震道数据可以包括:至少一个采样时间和所述采样时间处的振幅。在所述CMP道集数据中,所述偏移距、所述采样时间与所述振幅一一对应。
根据所述多个共中心点道集数据,可以确定所述第一剖面区域的初始层速度场。具体地,根据预设的一组均方根速度和所述第一剖面区域的一个CMP道集数据,可以确定该CMP道集数据中零偏移距对应各个采样时间处的初始均方根速度。可以利用迪克斯(Dix)公式将所述初始均方根速度转换为所述各个采样时间处的初始层速度。例如,所述Dix公式可以表示为:
公式(1)中,vn表示零偏移距对应第n个采样时间处的初始层速度,t0,n表示零偏移距对应的第n个采样时间,vrms,n表示零偏移距对应第n个采样时间处的初始均方根速度。可以对所述各个采样时间处的初始层速度进行时深转换处理,可以得到该CMP道集数据中零偏移距对应各个地层深度处的初始层速度,即该CMP道集数据的零偏移距对应地表位置在所述地层深度维度上的一维初始层速度。对于所述第一剖面区域的各个CMP道集数据,根据其零偏移距对应地表位置在所述地层深度维度上的一维初始层速度,可以确定所述第一剖面区域的初始层速度场。
进一步地,所述根据预设的一组均方根速度和所述第一剖面区域的一个CMP道集数据,确定该CMP道集数据的零偏移距对应各个采样时间处的初始均方根速度。具体可以包括,对于所述第一剖面区域的一个CMP道集数据,根据预设的一组均方根速度,可以对该CMP道集数据进行动校正处理。所述预设的一组均方根速度可以包括:步长为1米/秒的多个均方根速度。所述步长为1米/秒的多个均方根速度的速度范围可以为10~20米/秒。其中,根据所述预设的一组均方根速度中一个均方根速度,可以对该CMP道集数据进行一次动校正处理,根据每一次动校正处理后的CMP道集数据,可以得到一组动校正处理后的CMP道集数据。每一次动校正处理的过程可以是类似的。可以计算每一次动校正处理后的CMP道集数据中各偏移距对应第一采样时间处的振幅之间的相似系数。所述第一采样时间可以是该CMP道集数据中零偏移距对应的任一采样时间。这样,一个均方根速度可以对应计算得到一个相似系数,所述预设的一组均方根速度可以得到多个相似系数。可以将所述计算得到的多个相似系数中最大相似系数对应的均方根速度作为该CMP道集数据中零偏移距对应第一采样时间处的初始均方根速度。
进一步地,可以采用下述公式计算每一次动校正处理后的CMP道集数据中各个偏移距对应第一采样时间处的振幅之间的相似系数:
公式(2)中,S表示所述各个偏移距对应第一采样时间处的振幅之间的相似系数,t0表示所述第一采样时间,λ表示预设时窗窗长,xi表示所述动校正处理后的CMP道集数据中第i个地震道数据对应的偏移距,d(t,xi)表示所述动校正处理后的CMP道集数据中偏移距xi对应采样时间t处的振幅,N表示所述CMP道集数据中地震道数据的个数。所述预设时窗窗长可以为5个采样时间间隔。所述采样时间间隔可以是所述地震数据的采样时间间隔。所述采样时间间隔的取值可以为1毫秒。
步骤S102:基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域。
具体地,所述第一剖面区域中的异常地质区域可以包括:初始层速度大于所述预设层速度阈值的第一地层位置。所述第一地层位置可以是所述第一剖面区域中第一层段处的一个地层位置。所述第一层段可以是所述第一剖面区域中具有一定厚度的层段。所述预设层速度阈值可以是所述初始层速度场中第一层段处的平均层速度。
例如,例如,图2是本申请实施例目的工区中某一剖面区域的初始层速度场的剖面示意图。图2中横坐标表示目的工区的主测线维度或联络测线维度的地表位置,单位为米。图2中纵坐标表示目的工区的地层深度,单位为米。图2中的灰度值表示初始层速度,单位为米/秒。图2中黑色框内区域可能存在所述异常地质区域。如图2所示,由于异常地质区域的存在,所得到的初始层速度场中异常地质区域的下伏地层位置的层速度受到较大干扰,甚至出现畸变,导致这些地层位置处的层速度的准确度较低。
步骤S103:基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据。
所述第一种偏移距地震数据可以是与所述异常地质区域对应的地震数据。例如可以包括:所述第一剖面区域的多个共中心点道集数据中穿过所述异常地质区域的地震射线对应的地震道数据。其中,所述第一剖面区域的一个共中心点道集数据中穿过所述异常地质区域的地震射线对应的地震道数据,可以称为所述第一种偏移距地震数据中的一个共中心点分道集数据。
所述第二种偏移距地震数据可以是与所述第一剖面区域中非异常地质区域对应的地震数据。例如可以包括:所述第一剖面区域的多个共中心点道集数据中穿过所述非异常地质区域的地震射线对应的地震道数据。其中,所述第一剖面区域的一个共中心点道集数据中穿过所述非异常地质区域的地震射线对应的地震道数据,可以称为所述第二种偏移距地震数据中的一个共中心点分道集数据。
基于所述异常地质区域,还可以将所述多个共中心点道集数据中的偏移距划分为第一种偏移距和第二种偏移距。所述第一种偏移距可以是与所述第一种偏移距地震数据对应的偏移距。所述第二种偏移距是与所述第二种偏移距地震数据对应的偏移距。
例如,图3是本申请实施例第一剖面区域的某一共中心点道集数据对应的偏移距划分示意图。图3中横坐标表示目的工区的主测线维度或联络测线维度的地表位置,单位为米。图3中纵坐标表示目的工区的地层深度,单位为米。图3中的灰度值表示初始层速度,单位为米/秒。图3中黑色箭头指示的位置为某一共中心点道集数据对应的共中心点位置,区域1和3为均匀水平层状介质,区域2中存在异常地质区域,虚线框内的区域为第一剖面区域中除区域1、2和3以外的剩余区域。图3中共中心点道集数据对应的偏移距被划分为X1、X2和X3。其中,X2可以表示所述第一种偏移距,X1和X3可以表示所述第二种偏移距。
步骤S104:根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场。
所述第一偏移距区域可以包括:穿过所述异常地质区域的地震射线所覆盖的区域。所述第一偏移距区域是与所述第一种偏移距地震数据对应的区域。所述第二偏移距区域可以包括:未穿过所述异常地质区域的地震射线所覆盖的区域。所述第二偏移距区域是与所述第二种偏移距地震数据对应的区域。
根据所述第一种偏移距地震数据可以确定所述第一剖面区域中第一偏移距区域的目标层速度场。具体地,根据预设的一组均方根速度和所述第一种偏移距地震数据中第一共中心点分道集数据,可以确定所述第一共中心点分道集数据中各个采样时间处的目标均方根速度。所述第一共中心点分道集数据可以是所述第一种偏移距地震数据中的一个共中心点分道集数据。可以利用迪克斯公式将所述目标均方根速度转换为所述第一共中心点分道集数据中各个采样时间处的目标层速度。可以对所述第一共中心点分道集数据中各个采样时间处的目标层速度进行时深转换处理,得到所述共中心点分道集数据中各个地层深度处的目标层速度。根据所述第一共中心点分道集数据中各个地层深度处的目标层速度,可以确定所述第一偏移距区域的目标层速度场。所述第一偏移距区域内同一地层深度的不同地层位置处的目标层速度相同。
根据预设的一组均方根速度和所述第一种偏移距地震数据中第一共中心点分道集数据,可以确定所述第一共中心点分道集数据中各个采样时间处的目标均方根速度。具体地,所述确定所述第一共中心点分道集数据中各个采样时间处的目标均方根速度的方法与步骤S101中确定CMP道集数据的零偏移距对应各个采样时间处的初始均方根速度的方法类似,此处不再一一赘述。
根据所述第一共中心点分道集数据中各个地层深度处的目标层速度,可以确定所述第一偏移距区域的目标层速度场。具体地,根据所述第一共中心点分道集数据中各个地层深度处的目标层速度,可以确定所述第一共中心点分道集数据对应的目标层速度场。可以将所述第一共中心点分道集数据对应的目标层速度场作为所述第一偏移距区域的目标层速度场。
在另一种实施方式中,可以重复所述确定第一共中心点分道集数据中各个地层深度处的目标层速度的步骤,可以得到所述第一种偏移距地震数据的每一个共中心点分道集数据对应的目标层速度场,可以将这些目标速度场的平均值作为所述第一偏移距区域的目标层速度场。如此,可以进一步提高所确定的第一偏移距区域的层速度场的准确度。
根据所述第二种偏移距地震数据可以确定所述第一剖面区域中第二偏移距区域的目标层速度场。具体地,所述确定所述第二偏移距区域的目标层速度场的方法与确定所述第一偏移距区域的目标层速度场的方法类似,此处不再一一赘述。
步骤S105:根据所述第一偏移距区域的目标层速度场和所述第二偏移距区域的目标层速度场,确定所述第一剖面区域的目标层速度场。
具体地,根据所述第二偏移距区域的目标层速度场,可以确定第一剖面剩余区域的目标层速度场。所述第一剖面剩余区域可以表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域。例如,所述第一剖面剩余区域可以为如图3所示的剩余区域。所述第一剖面剩余区域与所述第二偏移距区域在同一地层深度的不同地层位置处的目标层速度相同。可以将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
例如,图4是本申请实施例中基于某一个共中心点道集数据确定的第一剖面区域的目标层速度场的剖面示意图。图4中横坐标表示目的工区的主测线维度或联络测线维度的地表位置,单位为米。图4中纵坐标表示目的工区的地层深度,单位为米。图4中的灰度值表示目标层速度,单位为米/秒。图4中黑色箭头指示的位置为共中心点道集数据对应的共中心点位置。图4中虚线表示与图3中相对应的偏移距划分边界线。如图4所示,相比图2中的初始层速度场,本图中的目标层速度场的准确度较高。
例如,图5是本申请实施例中基于第一剖面区域的各个共中心点道集数据确定的第一剖面区域的目标层速度场的剖面示意图。图5中横坐标表示目的工区的主测线维度或联络测线维度的地表位置,单位为米。图5中纵坐标表示目的工区的地层深度,单位为米。图5中的灰度值表示目标层速度,单位为米/秒。如图5所示,基于第一剖面区域的各个共中心点道集数据确定的第一剖面区域的目标层速度场,其准确度进一步得到提高。
所述确定速度场的方法实施例,首先,基于预设层速度阈值和所述目的工区中第一剖面区域初始层速度场,可以确定所述第一剖面区域中的异常地质区域;基于所述异常地质区域,可以将所述第一剖面区域的第一共中心点道集数据划分为第一偏移距地震数据和第二偏移距地震数据;所述第一偏移距地震数据可以是与所述异常地质区域对应的地震数据;所述第二偏移距地震数据可以是与所述第一剖面区域中非异常地质区域对应的地震数据;然后,根据所述第一偏移距地震数据可以确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二偏移距地震数据可以确定所述第一剖面区域中第二偏移距区域的目标层速度场;最后,根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,可以确定所述第一剖面区域的目标层速度场。如此,对于复杂结构的目的工区,其异常地质区域的层速度是基于穿过异常地质区域的地震射线对应的地震道数据所确定的,非异常地质区域的层速度是基于穿过非异常地质区域的地震射线对应的地震道数据所确定的,均不是基于穿过异常地质区域和非异常地质区域的地震射线对应的地震道数据所确定的平均层速度,从而可以提高最终确定的目的工区的层速度场的准确度。
图6是本申请确定速度场的装置实施例的组成结构图。如图6所示,所述确定速度场的装置可以包括:初始层速度场确定模块100、异常地质区域确定模块200、偏移距地震数据划分模块300、偏移距区域速度确定模块400和目标层速度确定模块500。
所述初始层速度场确定模块100,可以用于获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,可以确定所述第一剖面区域的初始层速度场。
所述异常地质区域确定模块200,可以用于基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域。
所述偏移距地震数据划分模块300,可以用于基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据;所述第一种偏移距地震数据可以是与所述异常地质区域对应的地震数据;所述第二种偏移距地震数据可以是与所述第一剖面区域中非异常地质区域对应的地震数据。
所述偏移距区域速度确定模块400,可以用于根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场。
所述目标层速度确定模块500,可以用于根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,确定所述第一剖面区域的目标层速度场。
图7是本申请确定速度场的装置实施例中目标层速度确定模块的组成结构图。如图7所示,图6中目标层速度确定模块500可以包括:剩余区域速度确定模块510和剖面区域速度确定模块520。
所述剩余区域速度确定模块510,可以用于根据所述第二偏移距区域的目标层速度场,确定第一剖面剩余区域的目标层速度场;所述第一剖面剩余区域可以表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域。
所述剖面区域速度确定模块520,可以用于将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
所述确定速度场的装置实施例与确定速度场的方法实施例相对应,可以实现本申请的方法实施例,并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(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。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种确定速度场的方法,其特征在于,包括:
获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,确定所述第一剖面区域的初始层速度场;
基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域;
基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据;所述第一种偏移距地震数据是与所述异常地质区域对应的地震数据;所述第二种偏移距地震数据是与所述第一剖面区域中非异常地质区域对应的地震数据;
根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场;
根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,确定所述第一剖面区域的目标层速度场。
2.根据权利要求1所述的一种确定速度场的方法,其特征在于,所述基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域,包括:
所述第一剖面区域中的异常地质区域包括:初始层速度大于所述预设层速度阈值第一地层位置;所述第一地层位置是所述第一剖面区域中第一层段处的一个地层位置。
3.根据权利要求2所述的一种确定速度场的方法,其特征在于,所述预设层速度阈值是所述初始层速度场中第一层段处的平均层速度。
4.根据权利要求1所述的一种确定速度场的方法,其特征在于,所述根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,包括:
根据预设的一组均方根速度和所述第一种偏移距地震数据中第一共中心点分道集数据,确定所述第一共中心点分道集数据中各个采样时间处的目标均方根速度;
将所述目标均方根速度转换为所述第一共中心点分道集数据中各个采样时间处的目标层速度;对所述第一共中心点分道集数据中各个采样时间处的目标层速度进行时深转换处理,得到所述第一共中心点分道集数据中各个地层深度处的目标层速度;
根据所述第一共中心点分道集数据中各个地层深度处的目标层速度,确定所述第一偏移距区域的目标层速度场。
5.根据权利要求4所述的一种确定速度场的方法,其特征在于,所述预设的一组均方根速度包括:步长为1米/秒的多个均方根速度;所述步长为1米/秒的多个均方根速度的速度范围为10~20米/秒。
6.根据权利要求4所述的一种确定速度场的方法,其特征在于,所述第一偏移距区域内同一地层深度的不同地层位置处的目标层速度相同。
7.根据权利要求1所述的一种确定速度场的方法,其特征在于,所述根据所述第一偏移距区域的目标层速度场和所述第二偏移距区域的目标层速度场,确定所述第一剖面区域的目标层速度场,包括:
根据所述第二偏移距区域的目标层速度场,确定第一剖面剩余区域的目标层速度场;所述第一剖面剩余区域表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域;
将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
8.根据权利要求7所述的一种确定速度场的方法,其特征在于,所述第一剖面剩余区域与所述第二偏移距区域在同一地层深度的不同地层位置处的目标层速度相同。
9.一种确定速度场的装置,其特征在于,所述装置包括:初始层速度场确定模块、异常地质区域确定模块、偏移距地震数据划分模块、偏移距区域速度确定模块和目标层速度确定模块;其中,
所述初始层速度场确定模块,用于获取目的工区中第一剖面区域的多个共中心点道集数据;根据所述多个共中心点道集数据,确定所述第一剖面区域的初始层速度场;
所述异常地质区域确定模块,用于基于预设层速度阈值和所述初始层速度场,确定所述第一剖面区域中的异常地质区域;
所述偏移距地震数据划分模块,用于基于所述异常地质区域,将所述第一剖面区域的多个共中心点道集数据划分为第一种偏移距地震数据和第二种偏移距地震数据;所述第一种偏移距地震数据是与所述异常地质区域对应的地震数据;所述第二种偏移距地震数据是与所述第一剖面区域中非异常地质区域对应的地震数据;
所述偏移距区域速度确定模块,用于根据所述第一种偏移距地震数据确定所述第一剖面区域中第一偏移距区域的目标层速度场,以及根据所述第二种偏移距地震数据确定所述第一剖面区域中第二偏移距区域的目标层速度场;
所述目标层速度确定模块,用于根据所述第一偏移距区域的层速度场和所述第二偏移距区域的层速度场,确定所述第一剖面区域的目标层速度场。
10.根据权利要求9所述的一种确定速度场的装置,其特征在于,所述目标层速度确定模块,包括:剩余区域速度确定模块和剖面区域速度确定模块;其中,
所述剩余区域速度确定模块,用于根据所述第二偏移距区域的目标层速度场,确定第一剖面剩余区域的目标层速度场;所述第一剖面剩余区域表示所述第一剖面区域中除所述第一偏移距区域和所述第二偏移距区域以外的区域;
所述剖面区域速度确定模块,用于将所述第一偏移距区域的目标层速度场、所述第二偏移距区域的目标层速度场和所述第一剖面剩余区域的目标层速度场所构成的速度场作为所述第一剖面区域的目标层速度场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710328456.4A CN107179546B (zh) | 2017-05-11 | 2017-05-11 | 一种确定速度场的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710328456.4A CN107179546B (zh) | 2017-05-11 | 2017-05-11 | 一种确定速度场的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107179546A true CN107179546A (zh) | 2017-09-19 |
CN107179546B CN107179546B (zh) | 2019-01-18 |
Family
ID=59832244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710328456.4A Active CN107179546B (zh) | 2017-05-11 | 2017-05-11 | 一种确定速度场的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107179546B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109188522A (zh) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | 速度场构建方法及装置 |
CN109655884A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 一种速度拾取异常点去除方法及系统 |
CN109884693A (zh) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | 自适应走向速度谱求取方法和系统 |
CN114114417A (zh) * | 2020-08-25 | 2022-03-01 | 中国石油天然气股份有限公司 | 地震成像方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839998A (zh) * | 2009-03-18 | 2010-09-22 | 中国石油天然气集团公司 | 一种高精度的叠前深度偏移方法 |
WO2012005958A3 (en) * | 2010-06-29 | 2012-07-05 | Saudi Arabian Oil Company | Migration velocity analysis of seismic data using common image cube and green's functions |
CN105388520A (zh) * | 2015-10-22 | 2016-03-09 | 中国石油化工股份有限公司 | 一种地震资料叠前逆时偏移成像方法 |
-
2017
- 2017-05-11 CN CN201710328456.4A patent/CN107179546B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101839998A (zh) * | 2009-03-18 | 2010-09-22 | 中国石油天然气集团公司 | 一种高精度的叠前深度偏移方法 |
WO2012005958A3 (en) * | 2010-06-29 | 2012-07-05 | Saudi Arabian Oil Company | Migration velocity analysis of seismic data using common image cube and green's functions |
CN105388520A (zh) * | 2015-10-22 | 2016-03-09 | 中国石油化工股份有限公司 | 一种地震资料叠前逆时偏移成像方法 |
Non-Patent Citations (2)
Title |
---|
易远元 等: "三维地震速度场建立技术——以饶阳凹陷河间南地区为例", 《石油学报》 * |
李丽贤 等: "塔里木盆地孔雀河地区速度场建立及应用探讨", 《中国西部油气地质》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109655884A (zh) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | 一种速度拾取异常点去除方法及系统 |
CN109188522A (zh) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | 速度场构建方法及装置 |
CN109188522B (zh) * | 2018-10-09 | 2020-02-14 | 中国石油天然气股份有限公司 | 速度场构建方法及装置 |
CN109884693A (zh) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | 自适应走向速度谱求取方法和系统 |
CN109884693B (zh) * | 2019-01-18 | 2021-10-15 | 潜能恒信能源技术股份有限公司 | 自适应走向速度谱求取方法和系统 |
CN114114417A (zh) * | 2020-08-25 | 2022-03-01 | 中国石油天然气股份有限公司 | 地震成像方法及装置 |
CN114114417B (zh) * | 2020-08-25 | 2024-03-01 | 中国石油天然气股份有限公司 | 地震成像方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107179546B (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Panning et al. | Importance of crustal corrections in the development of a new global model of radial anisotropy | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN107179546B (zh) | 一种确定速度场的方法及装置 | |
CN107092031B (zh) | 一种声波测井数据校正方法及装置 | |
CN106324676A (zh) | 一种确定断层封闭性的方法及装置 | |
CN106054251B (zh) | 一种初至波拾取方法及装置 | |
CN106842304B (zh) | 一种叠前深度偏移方法及装置 | |
CN108196304B (zh) | 一种多次波压制方法及装置 | |
CN106353818B (zh) | 一种地震数据时差校正的方法及系统 | |
CN106896406A (zh) | 基于阻抗域属性预测优质储层的方法和装置 | |
CN104375178B (zh) | 一种碳酸盐岩缝洞储层预测方法和装置 | |
CN107193045A (zh) | 一种地震数据处理方法及装置 | |
CN106646611A (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN104391319B (zh) | 一种地震资料采集系统的确定方法及装置 | |
CN107144878A (zh) | 一种断层识别方法及装置 | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
CN109143375A (zh) | 一种确定钻井方位的方法及装置 | |
CN106291691A (zh) | 一种地震偏移成像方法及装置 | |
CN109100797A (zh) | 一种确定致密储层含气性的方法及装置 | |
CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
CN106772596B (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN107656309B (zh) | 一种确定断层断距的方法及装置 | |
CN106908842B (zh) | 一种确定单井沉积微相曲线的方法及装置 | |
CN107219562B (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 |