CN111650638A - 一种地震波旅行时间计算方法 - Google Patents
一种地震波旅行时间计算方法 Download PDFInfo
- Publication number
- CN111650638A CN111650638A CN202010433784.2A CN202010433784A CN111650638A CN 111650638 A CN111650638 A CN 111650638A CN 202010433784 A CN202010433784 A CN 202010433784A CN 111650638 A CN111650638 A CN 111650638A
- Authority
- CN
- China
- Prior art keywords
- travel time
- calculation formula
- basic
- time calculation
- formula
- 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
Images
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. analysis, for interpretation, for correction
Abstract
本发明涉及地震勘探数据处理技术领域,尤其涉及一种地震波旅行时间计算方法;包括以下步骤:S1,确定基本旅行时间计算公式基本形式;S2,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似;S3,对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;本发明实施例通过对基本旅行时间计算公式进行更高阶的近似,然后对更高阶的所述基本旅行时间计算公式降幂处理输出旅行时间;操作简单,运算高效。
Description
技术领域
本发明涉及地震勘探数据处理技术领域,尤其涉及一种地震波旅行时间计算方法。
背景技术
地震波旅行时间在地震数据处理过程中起着重要作用;而部分公式在力求创新的同时加入了参与计算的射线条数、不同形状的反射器等限制,这对后续应用于实际地震数据会产生影响,因此,如何获取更为准确、有效的地震波旅行时间是我们当前需要研究的重点内容之一。
当前应用较多的计算旅行时间的方法有移位双曲法、广义时差近似法,其均能有效计算旅行时间。
现有技术的不足之处在于通过以上方法计算所得的旅行时间操作复杂,效率不高;旅行时间的计算还不准确,而且与初至时间还存在一定差距。
发明内容
为克服现有技术的不足,本发明实施例提供一种地震波旅行时间计算方法,能一定程度上减少计算所得旅行时间与初至时间差距,且操作简单,运算高效。
一方面,本发明实施例提供一种地震波旅行时间计算方法,包括以下步骤:
S1,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
S2,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
S3,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
另一方面,本发明实施例提供一种地震波旅行时间计算系统,包括:
计算优化模块,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
计算输出模块,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
本发明实施例提供一种地震波旅行时间计算方法及系统,通过对基本旅行时间计算公式进行更高阶的近似,然后对更高阶的所述基本旅行时间计算公式降幂处理输出旅行时间;操作简单,运算高效。
附图说明
为了更清楚地说明本发明的技术方案,下面将对本发明技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种地震波旅行时间计算方法流程示意图;
图2是平面水平层状介质模型图;
图3为实验所用模型图;
图4为DEC方法与移位双曲法、广义时差近似法的不同时间算法的旅行时间对比图;
图5为DEC方法与移位双曲法、广义时差近似法的不同旅行时间算法时差分析对比图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例一种地震波旅行时间计算方法流程示意图;如图1所示,S1,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
S2,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
S3,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
具体地,图2是平面水平层状介质模型图;基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;结合图2所示,炮点S和接收点G以偏移距x分开。波从炮点S经过N层界面反射到接收点G的传播时间记为T(x),可用以下参数化的公式表示:
和
其中射线参数p由p=sinβk/vk给出,βk为射线在第k层底面的入射角,vk为第k层的层速度,Δtk为第k层的单程旅行时间;
由公式(2)可知:对于第一层的反射,我们可以消除p并获得熟悉的双曲线旅行时间公式:
其中,T(0)为自激自收时间,x表示偏移距,v1为第一层的层速度。
经多层反射后的旅行时曲线T(x)可以用以下形式表示:
T2(x)=c0+c1x2+c2x4+c3x6+...
其中c0、c1、c2、c3为泰勒级数的系数,c0=T2(0),T(0)为自激自收时间,x表示偏移距。
然后,由公式(1)可知:将式(1)利用泰勒级数展开可扩展为式(4),
或者可表示为
其中,
根据式(5),我们可以得出x2,x4,x6等幂次方的绝对收敛级数;对于x2,我们得到:
其中,
且
..................
对于x2n有:
其中,
通过以上方法,我们可推出T2的幂级数展开式:
其中,
最后,我们将式(9)与式(13)带入式(3),比较等式两边p的幂级数可确定系数ci(i=1,2,...);
用以上理论写出T2(x)的基本形式,如式(15)所示:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+... (15)
其中,
c0=t2
以上,t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数,被定义为:
其次,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似;在2006年已经被提出的连续分式近似旅行时间平方公式为:
我们根据泰勒级数原理,对式(18)进行更高阶的近似,有:
式中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
最后,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;因为Hubral在1980年指出:在实际地震数据测试结果中,由T(x)的高阶项会得到一些几乎没有实用价值的冗余信息。在理论上,高阶项已经被数学上的复杂形式束缚住,通常不建议使用更高阶项做进一步近似处理。
基于以上分析,此时需要对公式中的高阶项进行处理。现利用降幂思维进行处理,得到如式(21)所示:
式中,S2可继续用式(23)-(24)表示:
S2=1+8η (23)
式中,η为各向异性参数,可用ε和δ表示,其中ε表示各向异性强度、δ表示NMO速度和反射振幅的变化,ε和δ均为Thomsen各向异性参数。
本发明实施例提供一种地震波旅行时间计算方法,通过对基本旅行时间计算公式进行更高阶的近似,然后对更高阶的所述基本旅行时间计算公式降幂处理输出旅行时间;操作简单,运算高效。
进一步地,所述步骤S3之后还包括对降幂处理后的所述基本旅行时间计算公式结合近似补偿方法形成降幂补偿公式:
其中,fb为初至时间,TDE(x)为降幂后的所述基本旅行时间计算公式计算出的旅行时间,n为模型层数。
具体地,从泰勒级数高阶近似角度考虑,这种对于公式的“硬”处理方法会改变高阶项中“高”的意义,使计算的旅行时间缺少细微的准确性。因此,对计算后的旅行时间进行补偿是一种弥补细节缺失问题的方式,本文采用一种近似补偿方法如式(25)所示:
最后,将式(21)与补偿方法相结合可形成降幂补偿(DescendingExponentCompensation,DEC)公式,记为DEC方法,如式(26)所示:
其中,fb为初至时间,TDE(x)为式(21)计算出的旅行时间,n为模型层数。
例如,图3为实验所用模型图。如图3所示,实验采用13层水平层状各向异性模型,其中道间距40m,接收器226个,最大深度为6000m,最大偏移距为9000m。初至时间由模型第一层上方的检波器来记录。
图4为DEC方法与移位双曲法、广义时差近似法的不同时间算法的旅行时间对比图。如图4所示,其中(a)深度为800m时旅行时间计算值(b)深度为2600m时旅行时间计算值;(c)深度为3500m时旅行时间计算值(d)深度为4600m时旅行时间计算值;
图5为DEC方法与移位双曲法、广义时差近似法的不同旅行时间算法时差分析对比图。如图5所示,其中(a)深度为800m的时差对比;(b)深度为2600m的时差对比;(c)深度为3500m的时差对比;(d)深度为4600m的时差对比;时差对比实验采用公式(27)来验证差异情况。
其中,以初至时间fb作为参考值,tc表示旅行时间公式计算所得时间。
由图4可知:①DEC计算方法同其他2种方法一样可以有效地计算出旅行时间。②尽管DEC方法的计算值在4个不同深度的近偏移距2000m处所显示的差异不大,但在中偏移距4500m、远偏移距8000m处均能明显看出优势所在;③DEC方法在旅行时间计算对比实验中表现更好。
由图5可知:①DEC方法在近偏移距2000m附近依然显示出与其他2种方法具有较小差异;②从3种方法计算出的时差的绝对值考虑,在中偏移距4500m、远偏移距8000m附近,DEC方法得出的时差值是最小的,这也说明其与初至时间的差距要小于其他方法;③DEC方法提供高了计算的准确性。本发明实施例提供一种地震波旅行时间计算方法,通过对基本旅行时间计算公式进行更高阶的近似,然后对更高阶的所述基本旅行时间计算公式降幂处理后与补偿方法相结合可形成降幂补偿,输出旅行时间;操作简单,运算高效;在计算地层旅行时间问题上DEC方法不但缩短了与初至时间的差距,还提高了计算的准确性。
基于上述实施例,本发明实施例提供了一种地震波旅行时间计算系统,包括:
计算优化模块,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示
计算输出模块,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
本发明实施例提供一种地震波旅行时间计算系统执行上述方法,通过对基本旅行时间计算公式进行更高阶的近似,然后对更高阶的所述基本旅行时间计算公式降幂处理后与补偿方法相结合可形成降幂补偿,输出旅行时间;操作简单,运算高效;在计算地层旅行时间问题上DEC方法不但缩短了与初至时间的差距,还提高了计算的准确性。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (4)
1.一种地震波旅行时间计算方法,其特征在于,包括以下步骤:
S1,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
S2,根据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
S3,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
3.一种地震波旅行时间计算系统,其特征在于,包括:
计算优化模块,确定基本旅行时间计算公式基本形式;所述基本旅行时间计算公式是根据泰勒级数原理展开式并利用各向异性介质模型得出基本旅行时间平方的表达式;所述基本旅行时间计算公式为:
T2(x)=c0+c1x2+c2x4+c3x6+c4x8+...
c0=t2
其中,x为偏移距;c0、c1、c2、c3、c4为泰勒级数的系数;t=T(0)为自激自收时间,v=vNMO表示NMO速度,Sk(k=2,3,4)为异质性参数;
据所述泰勒级数原理和连续分式近似旅行时间平方公式对所述基本旅行时间计算公式进行更高阶的近似,得到下式:
其中,参数B1是基于6阶泰勒展开和4阶泰勒展开的一种系数近似表示,参数B2是基于8阶泰勒展开和6阶泰勒展开的一种系数近似表示;
计算输出模块,以x2为基数对更高阶的所述基本旅行时间计算公式中偏移距的4次方和6次方进行同时降幂处理,输出旅行时间;降幂处理后的所述基本旅行时间计算公式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010433784.2A CN111650638B (zh) | 2020-05-21 | 2020-05-21 | 一种地震波旅行时间计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010433784.2A CN111650638B (zh) | 2020-05-21 | 2020-05-21 | 一种地震波旅行时间计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111650638A true CN111650638A (zh) | 2020-09-11 |
CN111650638B CN111650638B (zh) | 2022-07-05 |
Family
ID=72340643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010433784.2A Active CN111650638B (zh) | 2020-05-21 | 2020-05-21 | 一种地震波旅行时间计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111650638B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040015296A1 (en) * | 2000-07-17 | 2004-01-22 | Emmanuel Causse | Seismic processing with general non-hyperbolic travel-time corrections |
US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
CN102213769A (zh) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 |
CN102721977A (zh) * | 2012-05-31 | 2012-10-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 高精度转换波各向异性叠加速度分析方法 |
CN104076391A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法 |
WO2015155597A2 (en) * | 2014-04-07 | 2015-10-15 | Cgg Services Sa | Attenuating pseudo s-waves in acoustic anisotropic wave propagation |
WO2019071504A1 (zh) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | 一种基于两点射线追踪的地震走时层析反演方法 |
-
2020
- 2020-05-21 CN CN202010433784.2A patent/CN111650638B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040015296A1 (en) * | 2000-07-17 | 2004-01-22 | Emmanuel Causse | Seismic processing with general non-hyperbolic travel-time corrections |
US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
CN102213769A (zh) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 |
CN102721977A (zh) * | 2012-05-31 | 2012-10-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 高精度转换波各向异性叠加速度分析方法 |
WO2015155597A2 (en) * | 2014-04-07 | 2015-10-15 | Cgg Services Sa | Attenuating pseudo s-waves in acoustic anisotropic wave propagation |
CN104076391A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法 |
WO2019071504A1 (zh) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | 一种基于两点射线追踪的地震走时层析反演方法 |
Non-Patent Citations (2)
Title |
---|
孙祥娥 等: "利用Walkaway VSP和全方位地面地震数据求取各向异性参数", 《石油地球物理勘探》 * |
陈可洋: "地震波旅行时计算方法及其模型试验分析", 《石油物探》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111650638B (zh) | 2022-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111239802B (zh) | 基于地震反射波形和速度谱的深度学习速度建模方法 | |
CN101116008B (zh) | 用于针对正常时差拉伸影响进行地震数据真实相对振幅校正的方法和装置 | |
CN103149586B (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
US10670750B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
US9341729B2 (en) | Amplitude contrast seismic attribute | |
EP2689273A2 (en) | System and method for seismic data modeling and migration | |
CN101614826A (zh) | 三维地震数据处理中实现面元均化的方法和装置 | |
CN105353408A (zh) | 一种基于匹配追踪的Wigner高阶谱地震信号谱分解方法 | |
CN110542928A (zh) | 基于vti各向异性传播矩阵的地震响应模拟方法 | |
CN111722283B (zh) | 一种地层速度模型建立方法 | |
WO2016093875A1 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
CN105242313B (zh) | 一种弹性波逆时偏移极性反转的校正方法及系统 | |
CN113552624B (zh) | 孔隙度预测方法及装置 | |
CN111650638B (zh) | 一种地震波旅行时间计算方法 | |
CN107340537A (zh) | 一种p-sv转换波叠前逆时深度偏移的方法 | |
CN103790569B (zh) | 一种利用vsp时深关系对声波测井资料进行校正的方法 | |
Tocheport et al. | A systematic study of source time functions and moment tensors of intermediate and deep earthquakes | |
bin Waheed et al. | A holistic approach to computing first-arrival traveltimes using neural networks | |
CN111273349B (zh) | 一种用于海底浅部沉积层的横波速度提取方法及处理终端 | |
Waheed et al. | A holistic approach to computing first-arrival traveltimes using neural networks | |
Sun et al. | Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm | |
Bashir et al. | Preservation of seismic diffraction to enhance the resolution of seismic data | |
Jiang et al. | Full waveform inversion based on inversion network reparameterized velocity | |
CN116068644A (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 |