平码五不中公式规律
  • / 17
  • 下载费用:30 金币  

一种地下速度模型的构建方法.pdf

关 键 ?#21097;?/dt>
一种 地下 速度 模型 构建 方法
  专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
摘要
申请专利号:

CN201510524461.3

申请日:

2015.08.24

公开号:

CN106483559A

公开日:

2017.03.08

当前法律状态:

授权

有效性:

有权

法?#19978;?#24773;: 授权|||实质审查的生效IPC(主分类):G01V 1/28申请日:20150824|||公开
IPC分类号: G01V1/28 主分类号: G01V1/28
申请人: 中国石油化工股份有限公司; 中国石油化工股份有限公司石油物探技术研究院
发明人: 倪瑶; ?#25506;?#38596;
地址: 100728 ?#26412;?#24066;朝阳区朝阳门北大街22号
优?#28909;ǎ?/td>
专利代理机构: ?#26412;?#32895;宏知识产权代理有限公司 11372 代理人: 张文娟;朱绘
PDF完整版下载: PDF下载
法律状态
申请(专利)号:

CN201510524461.3

授权公告号:

||||||

法律状态公告日:

2018.11.23|||2017.04.05|||2017.03.08

法律状态类型:

授权|||实质审查的生效|||公开

摘要

本发明公开了一种地下速度模型的构建方法,属于油气勘探开发技术领域,避免了叠前地震同相轴的人工/半自动拾取的过程,加快了地下速度模型的构建过程。该地下速度模型的构建方法包括:获取共反射面元叠加算子;根据所述共反射面元叠加算子,获取地震波的运动学波场属性参数;基于建立的正则化层析方程,结合运动学波场属性参数和所建立的初始地下速度模型,构建地下速度模型。

权利要求书

1.一种地下速度模型的构建方法,其特征在于,包括:
获取共反射面元叠加算子;
根据所述共反射面元叠加算子,获取地震波的运动学波场属性参数;
基于建立的正则化层析方程,结合运动学波场属性参数和所建立的初始地下
速度模型,构建地下速度模型。
2.根据权利要求1所述的方法,其特征在于,基于建立的正则化层析方程,
结合运动学波场属性参数和所建立的初始地下速度模型,构建地下速度模型包
括:
?#34903;鑑、基于建立的正则化层析方程,结合运动学波场属性参数,计算模型
更新量;
?#34903;鑒、将当前的地下速度模型叠加模型更新量,获得更新的地下速度模型;
?#34903;鑓、正演计算,获得更新的地下速度模型的运动学波场属性参数;
?#34903;鑔、判断更新的地下速度模型的运动学波场属性参数与地震波观测运动
学波场属性参数之间的数据残差是否小于或等于预设残差,若是,执行?#34903;鑕,
若否,执行?#34903;鑖;
?#34903;鑕、输出该更新的地下速度模型为所构建的地下速度模型;
?#34903;鑖、根据更新的地下速度模型,更新正则化层析方程,将更新的地下速
度模型作为当前的地下速度模型,返回执行?#34903;鑑。
3.根据权利要求2所述的方法,其特征在于,根据更新的地下速度模型,
更新正则化层析方程包括:
获取数据与条件矩阵和更新的地下速度模型对应的弗雷歇导数矩阵;
获取与地下速度模型空间维度相同的单位矩阵、对B样条基函数系数的更新
量在横向和纵向上的一阶差分算子;
获取所设定的单位矩阵的乘数,获取两个一阶差分算子各自的乘数;
根据数据残差、弗雷歇导数矩阵、单位矩阵、对B样条基函数系数的更新量
在横向和纵向上的一阶差分算子和所获取的乘数,更新正则化层析方程。
4.根据权利要求3所述的方法,其特征在于,根据所述共反射面元叠加算
子,获取地震波的运动学波场属性参数包括:
对零偏地震数据上的一样点,设定不同的运动学属性参数对,所述运动学属
性参数对包括零偏射线由地下反射点传播至地表的角度与局部波前曲率半径;
利用所述共反射面元叠加算子,得到各运动学属性参数对在该样点的叠加能
量谱;
获取叠加能量谱值最大的运动学属性参数对;
根据所获取的运动学属性参数对,基于射线理论,获取适合于后续建立正则
化层析方程的地震波运动学波场属性参数。
5.根据权利要求4所述的方法,其特征在于,
共反射面元叠加算子涉及地表?#26800;?#21644;半偏移距下地震波传播的走时关系式,
该表达式不?#35272;?#20110;地下宏观速度模型。
6.根据权利要求5所述的方法,其特征在于,
运动学波场属性参数包括正入射点波到达地表的传播时间、波前曲率、横向
坐标以及慢度矢量的水平分量。

说明书

一种地下速度模型的构建方法

技术领域

本发明涉及油气勘探开发技术领域,具体地说,涉及一种地下速度模型的构
建方法。

背景技术

估计地下介质的速度(尤其是其低波数部分,也称宏观速度)分布是地震勘
探的核心问题之一,对地下构造偏移?#19978;瘛?#39044;测油气储层参数都?#35272;?#20110;一个准确
的速度模型。地震层析?#19978;?#25216;术具有较高的反演精度,是估计地下速度模型的重
要方法,其中应用最广泛的是反射层析?#19978;瘛?br />

反射层析具有效率高、容?#36164;?#29616;等优点,然而其需要在叠前地震道集上拾取
沿整个剖面连续分布的同相轴。这一操作十分耗时(特别是在三维情况下),且
对于低信噪比叠前地震数据而言,几乎是不可能完成的任务。

另外一类在工业?#26800;?#21040;广泛应用的宏观速度估计技术是偏移速度分析
(MVA)。该类技术将层析反演的原理与偏移?#19978;?#29983;成的共?#19978;?#28857;道集(Common
Image Gather,简称CIG)结合起来,?#26723;?#20102;速度估计技术对叠前地震数据信噪比
的要求,具有较高的实用价值。然而该类技术需要多次实施叠前深度偏移并人机
交互拾取共?#19978;?#28857;道集,这使得该类技术在三维情况下十分低效,一定程度上降
低了其实用价值。

发明内容

本发明的目的在于提供一种地下速度模型的构建方法,避免了叠前地震同相
轴的人工/半自动拾取的过程,加快了地下速度模型的构建过程。

本发明提供了一种地下速度模型的构建方法,包括:

获取共反射面元叠加算子;

根据所述共反射面元叠加算子,获取地震波的运动学波场属性参数;

基于建立的正则化层析方程,结合运动学波场属性参数和所建立的初始地下
速度模型,构建地下速度模型。

可选?#27169;?#22522;于建立的正则化层析方程,结合运动学波场属性参数和所建立的
初始地下速度模型,构建地下速度模型包括:

?#34903;鑑、基于建立的正则化层析方程,结合运动学波场属性参数,计算模型
更新量;

?#34903;鑒、将当前的地下速度模型叠加模型更新量,获得更新的地下速度模型;

?#34903;鑓、正演计算,获得更新的地下速度模型的运动学波场属性参数;

?#34903;鑔、判断更新的地下速度模型的运动学波场属性参数与地震波观测运动
学波场属性参数之间的数据残差是否小于或等于预设残差,若是,执行?#34903;鑕,
若否,执行?#34903;鑖;

?#34903;鑕、输出该更新的地下速度模型为所构建的地下速度模型;

?#34903;鑖、根据更新的地下速度模型,更新正则化层析方程,将更新的地下速
度模型作为当前的地下速度模型,返回执行?#34903;鑑。

可选?#27169;?#26681;据更新的地下速度模型,更新正则化层析方程包括:

获取数据与条件矩阵和更新的地下速度模型对应的弗雷歇(Fréchet)导数矩
阵;

获取与地下速度模型空间维度相同的单位矩阵、对B样条基函数系数的更新
量在横向和纵向上的一阶差分算子;

获取所设定的单位矩阵的乘数,获取两个一阶差分算子各自的乘数;

根据数据残差、弗雷歇导数矩阵、单位矩阵、对B样条基函数系数的更新量
在横向和纵向上的一阶差分算子和所获取的乘数,更新正则化层析方程。

可选?#27169;?#26681;据所述共反射面元叠加算子,获取地震波的运动学波场属性参数
包括:

对零偏地震数据上的一样点,设定不同的运动学属性参数对,所述运动学属
性参数对包括零偏射线由地下反射点传播至地表的角度与局部波前曲率半径;

利用所述共反射面元叠加算子,得到各运动学属性参数对在该样点的叠加能
量谱;

获取叠加能量谱值最大的运动学属性参数对;

根据所获取的运动学属性参数对,基于射线理论,获取适合于后续建立正则
化层析方程的地震波运动学波场属性参数。

可选?#27169;?#20849;反射面元叠加算子涉及地表?#26800;?#21644;半偏移距下地震波传播的走时
关系式,该表达式不?#35272;?#20110;地下宏观速度模型。

可选?#27169;?#36816;动学波场属性参数包括正入射点波到达地表的传播时间、波前曲
率、横向坐标以及慢度矢量的水平分量。

本发明带来了以下有益效果:本发明实施例提供了一种地下速度模型的构建
方法,该构件方法中利用共反射面元叠加算子来提取、观测地震波场的运动学波
场属性参数,并用于层析速度反演,避免了叠前地震同相轴的人工/半自动拾取,
理论上较反射层析更加适应于低信噪比叠前地震数据。

本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书
中变得显而?#20934;?#25110;者通过实施本发明而了解。本发明的目的和其他优点可通过
在说明书、权利要求书以及附图中所特别指出的结构?#35789;?#29616;和获得。

附图说明

为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需
要的附图做简单的介绍:

图1是本发明实施例中的构建速度模型的流程示意图;

图2是本发明实施例中的由s0到s1的透射过程的示意图;

图3是本发明实施例中的层析迭代的过程示意图;

图4是本发明实施例中的真实理论模型的示意图;

图5是本发明实施例中的初始速度模型的示意图;

图6是本发明实施例中的层析迭代1次后的速度模型的示意图;

图7是本发明实施例中的层析迭代10次后的速度模型的示意图;

图8是本发明实施例中的层析迭代20次后的速度模型的示意图;

图9是本发明实施例中的误差泛函下降曲线的示意图;

图10是本发明实施例中的部分炮道集的示意图;

图11是本发明实施例中的Kirchhoff叠前深度偏移结果的示意图。

具体实施方式

以下将结合附?#25216;?#23454;施例来详细说明本发明的实施方式,借此对本发明如何
应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实
施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的
各个特征可以相互结合,所形成的技术方案均在本发明的保护?#27573;?#20043;内。

本发明提供了一种地下速度模型的构建方法,如图1所示,该方法包括:

?#34903;鑃101、获取共反射面元叠加算子。

二维情况下,对地震数据中的某一零偏样点(ξm,T),其中T表示正入射点波
的传播时间,ξm是该零偏样点的地表?#26800;?#22352;标。共反射面元叠加算子是该零偏样
点邻域?#27573;?#20869;的?#26800;?#22352;标ξ与半偏移距h的函数,可由运动学属性参数对
(α,RNIP)表示如下:


其中,v0为近地表速度,RN是正入射点波的波前曲率半径。对于任一零偏样
点(ξm,T),其运动学波场属性d可通过对叠前数据进行多参数相干分析得到。式
(1)反映了?#26800;?半偏移距坐标系下地震波场的走时关系,可以看出共反射面元
叠加方法只需要得知近地表速度,是一种完全数据驱动的方法。?#26723;?#27880;意的是运
动学属性参数对有明确的物理意义,其中α,RNIP分别对应于零偏射线由地下反射
点传播至地表的角度与局部波前曲率半径,对应于法向入射点(Normal Incident
Point,简称NIP)波层析的数据分量,这里法向入射点波即为零偏射线。

?#34903;鑃102、根据所述共反射面元叠加算子,获取地震波的运动学波场属性参
数。

接下来,对零偏地震数据上任一样点(ξ,T),通过给定最大、最小值?#27573;В?br />等间隔地选取从而设定不同的运动学属性参数对(α,RNIP)i(i为?#25105;?#27491;整数)。对
每组运动学属性参数对,利用式(1)描述的共反射面元叠加算子进行叠加,得
到各组运动学属性参数对的叠加能量谱。

对样点(ξ,T),比较不同运动学属性参数对的叠加能量谱值,选取最大的叠
加能量值对应的运动学属性参数对作为该点的运动学属性参数对(α,RNIP)。

根据射线理论,可知有如下两个定量关系式:



其中,式(2)中的pξ为射线理论中慢度矢量的水平分量,代表了零偏数据
的同相轴斜?#21097;琈为地震波场的局部波前曲率。

结合式(2)和(3),对零偏地震数据上的?#25105;?#19968;个样点可以得到一组运动
学波场属性参数d=(T,M,ξ,pξ),即该运动学波场属性参数包括对NIP波到达地
表的传播时间T、波前曲率M、地表?#26800;?#27178;坐标?#25105;?#21450;慢度矢量的水平分量pξ。
本发明中,运动学波场属性参数的提取过程完全自动化,并作为运动学波场属性
参数作为正则化层析方程的建立与求解的输入。

由于引入了同相轴斜率与波前曲率参数用来约束速度模型,本发明实施例
中,反演得到的速度模型具有较高的精度,满足后续偏移?#19978;?#30340;要求。

?#34903;鑃103、基于建立的正则化层析方程,结合运动学波场属性参数和所建立
的初始地下速度模型,构建地下速度模型。

本发明的层析反演属于数据拟合类的速度估计方法,即?#19994;?#19968;个如式(5)定义
的速度模型m,使得正演计算得到的数据dmod=f(m)与式(4)定义的拾取数据d的
误差最小化。f表示层析的非线性正算子,即反射点至地表方向的初值射线追踪。
显然,对于一个零偏射线的反射,假如反射点的坐标、出射角以及速度模型正确
时,由反射点以出射角出射射线,当到达观测面时,零偏射线刚好到以正确的出
射方向、正确的走时到达正确的地表位置。

本发明实施例中的层析数据分量和模型分量如下:

d=(T,M,ξ,pξ)i,i=1,...,ndata (4)

m=[(x,z,θ)i,i=1,...,ndata;vjk,j=1,...,nx,k=1,...,nz] (5)

式(4)中数据分别为拾取到的ndata个运动学波场属性参数。同样,每个数据点
d对应于一个零偏反射过程,在模型中对应于式(5)的一组模型参数(x,z,θ)i,即反
射点横坐标、纵坐标、出射角。

以二范数衡量数据残差,本发明实施例中的构建地下速度模型的过程可归结
为如下泛函的求极值问题:


其中,Δd(m)=d-f(m)。对角阵也称为数据协方差矩阵,起加权不同数
据分量的作用。

由于反问题固有的病态性与多解性,本发明实施例中,引入额外的先验模型
进?#24615;?#26463;,以保证反演结果的可靠性。改造式(7)的数据残差泛函如下:


其中为较小的正实数(这三个参数分布在0到1之间,
这三个参数也需要用户通过测试自定义给出,当反演的初始模型(形式如式(5)
定义)越好、数据?#20998;试?#39640;则这几个参数给定越小,反之给定越大),用于加权
数据拟合项与规则化项之间的权重;mr为某给定的先验模型,先验模型的形式与
式(5)定义的速度模型一致,实施反演时需要用户自定义给出先验模型。给定
先验模型的方法很多,本发明实施例取先验模型mr=0。

为阻尼项,这项的引入使得每次迭代的模型更新量较小。
和分别约束速度模型在横向和纵向上的光滑程
度,D1、D2分别是对B样条基函数系数更新量在横向和纵向上的一阶差分算子。

对一个有7个数的数列而言,其一阶差分算子定义为如下6×7的矩阵,


同理,对于有N个点的数列,一阶差分算子是形如上式的(N-1)×N矩阵。

根据反演的局部最优化理论,式(7)刻画的最优化问题可等效于如下线性
方程组的迭代求解,通过多次迭代达到反演解收敛的目?#27169;?#20854;中方程组(8)称
为本发明实施例层析反演的正则化层析方程。


式(8)中,Δd是反演的数据残差,Δd(m)=d-f(m);I是一个与模型空间
维度相同的单位矩阵其中Nmod el是本发明层析反
演所用到的模型分量的维度;F是Fréchet(弗雷歇)导数矩阵,是一个Ndata×Nmod el
的矩阵,其中Ndata是本发明实施例中层析反演所用到的数据分量的维度;D是数
据预条件矩阵,用来加权本发明层析反演的不同量纲的数据空间分量,对每个运
动学波场属性参数(T,M,ξ,pξ)i,D由如下形式定义:


本发明实施例中取(σT,σM,σξ,σp)=(103,109,10,106)。

其中f表示透射射线追踪算子(在油气地球物理勘探领域较为成熟),m是
式(5)定义的模型分量;F为在mn处取值的Fréchet导数矩阵,其中
Fij表示Fréchet导数矩阵的第i?#26800;趈列的元素,具体如下:

本发明实施例中层析的Fréchet导数的计算基于射线扰动理论
与射线传播矩阵的性?#21097;?#32771;虑到射线扰动理论在射线中心坐标系下的简洁形式,
本发明实施例在该坐标系下经过推导,得到了层析反演所需的Fréchet导数。

如图2所示,考虑如下一个由s0到s1的透射过程,有两边全微分
得:


由式(1)即可计算出

的计算可以由射线扰动理论求得。射线扰动理论本质上是建立
了观测射线场的扰动量、初始射线场以及速度的扰动量之间的线性关系,由该线
性关系可以换算出(ξ,pξ,M)的扰动量与(x0,z0,θ,v)的扰动量的线性关系式。

射线扰动理论在射线中心坐标系下的形式最为简单,此处推导时先在射线中
心坐标系下应用射线扰动理论,然后根据直角坐标与射线中心坐标的几何关?#26723;?br />出(ξ,pξ,M)与(x0,z0,θ,v)的线性关系式,即可求得

根据射线扰动理论,扰动坐标q与相应的慢度矢量分量p可由下式所示:


其中:Π(s1,s0)称为从s0点传播
到s1点的射线传播矩阵,其为一2×2的矩阵,记为:


(Q1(s1,s0),P1(s1,s0))T与(Q2(s1,s0),P2(s1,s0))T分别为傍轴射线追踪系统初值为
(1,0)T与(0,1)T的解,物理上分别对应于归一化的平面波震源与点震源。射线中心
坐标系下的傍轴射线追踪系统为:


射线传播矩阵有如下性?#21097;?br />

Π(s1,s)=Π(s1,s0)Π(s0,s)=Π(s1,s0)Π-1(s,s0)

可知:


由A(2)、A(3)得:


其中,

考虑到射线中心坐标系与直角坐标系的几何关系有:



q0=cosθΔx-sinθΔz A(5)


s0=sinθΔx+cosθΔz

上述的几何关系A(5)代入A(4)得到:



另外,由得,


而根据射线理论,有:


S=ΔS1(Δv)+ΔS2(Δq,Δp) A(10)



根据A(6)、A(7)即可求得导数联立A(5)、A(8)至A(12)可求得
本发明实施例中用B样条插?#24403;?#36798;速度模型,反演时在迭代中更新
的速度参数为B样条基函数的系数,因此需要求取数据分量对B样条基函数系数
的导数,这只需要将代入式A(6)、A(7)、A(11)即得到
相应于B样条基函数系数的全微分表达式即可。其中,式A(12)中的vq、vqq、
vqqq分别代表速度模型在垂直于射线传播方向上的一阶方向导数、二阶方向导数
与三阶方向导数。

本发明利用最小二乘QR方法(LSQR)求解矩阵方程组(8)(方程组(8)中Δm
是未知量,其他均是已知量),该方法是一种迭代的方法,可以在最小二乘意义
下高效地求解大规模稀疏矩阵。

具体?#27169;?#21453;演过程中,每?#20540;?#20195;用运动学波动属性参数d来计算方程组(8)
的右端项,其中f(m)是利用上一?#20540;?#20195;的模型正演模拟得到的。

具体?#27169;?#22914;图3所示,设当前的速度模型分量为mk,则第k+1次迭代的过
程为:

?#34903;鑑、基于建立的正则化层析方程,结合运动学波场属性参数,计算模型
更新量Δmk+1。

?#34903;鑒、将当前的地下速度模型mk叠加模型更新量Δmk+1,获得更新的地下
速度模型mk+1。

?#34903;鑓、正演更新的地下速度模型,获得更新的地下速度模型的运动学波场
属性参数f(mk+1)。

?#34903;鑔、判断更新的地下速度模型的运动学波场属性参数与地震波的运动学
波场属性参数之间的数据残差是否小于或等于预设残差,若是,执行?#34903;鑕,若
否,执行?#34903;鑖。

?#34903;鑕、输出该更新的地下速度模型为所构建的地下速度模型。

?#34903;鑖、根据更新的地下速度模型,更新正则化层析方程,将更新的地下速
度模型作为当前的地下速度模型,返回执行?#34903;鑑。

获取数据与条件矩阵和更新的地下速度模型对应的弗雷歇导数矩阵;获取与
地下速度模型空间维度相同的单位矩阵、对B样条基函数系数的更新量在横向和
纵向上的一阶差分算子;获取所设定的单位矩阵的乘数,获取两个一阶差分算子
各自的乘数;根据数据残差、以及上述获取的各参数,更新正则化层析方程。

显然,整个层析反演迭代过程全自动化,不需要进行人工?#31245;ぃ?#36825;提高了整
个速度估计过程的效?#21097;?#20943;少了人工工作量,缩短了处理周期。

基于上述构建速度模型的方法,在本发明实施例的具体实施场景中,图4为
原始新疆岩性模型,模型大小为横向30公里、纵向6公里。为了验证层析对强
非均质性模型的适应能力,对图4进行光滑并运用运动学与动力学射线追踪为其
提供理想的数据输入。二维B样条基函数之间的间隔剖分为横向200m、纵向
100m。图5为初始给定的常梯度模型,所前文所述,只需将B样条基函数系数在
纵向给成相应的常梯度即可完成。由地表的数据分量在初始模型中由地表位置向
下出射射线,当单程旅行时用完时停止,可以完成对反射点位置坐标和出射角的
初始化。需要说明的是,该层析反演方法在理想无噪音数据中进行。

输入数据分量与初始模型分量,便可开始层析反演的模型迭代更新过程。图
6为第一次迭代更新后的速度场,可以看出,表层的速度得到有效修改。随着迭
代的进行,速度模型逐渐收敛,迭代10轮后的速度模型如图7所示,地表?#20102;?br />公里?#27573;?#20869;的速度模型得到了较好的?#25351;矗?#28145;部的模型由于采集孔径问题覆盖不
足,未能完全?#25351;?#20986;来。迭代20次后的速度模型如图8所示,反演此?#24065;?#22522;本
收?#30149;?br />

图9为误差泛函随着迭代次数增加的下降曲线,误差泛函取以10为底的对
数,前15次误差泛函下?#21040;?#24555;,后面的迭代中误差泛函下降十?#21482;?#24930;,说明反
演已收敛至极值处。本实验采用的规则化因子分别为εd=0.02,
另外一种策略是在迭代过程中逐步下降规则化因子,这使得在保证反
演稳定的同?#22791;?#22810;的模型?#38468;詒换指?#20986;来,但规则化因子也不能下?#26723;?#26576;个指定
的阈值一下,以防止反演再次出现不稳定。本实验说明本发明层析对于强非均质
性模型具有较强的适应力。

图10为采用有限差分法对新疆岩性模型进行正演的炮道集(显示部分)。以反
演第20次迭代的结果作为速度输入,实施Kirchhoff叠前深度偏移。偏移叠?#27833;?br />像如图11所示,可以看出基本的构造形态已得到较好地?#19978;瘢?#20063;进一步验证了
本文算法的精度。

综上,本发明实施例提供了一种地下速度模型的构建方法,该构件方法中利
用共反射面元叠加算子来提取、观测地震波场的运动学波场属性参数,并用于层
析速度反演,避免了叠前地震同相轴的人工/半自动拾取,理论上较反射层析更加
适应于低信噪比叠前地震数据;同时,除反射地震走时外,将同相轴斜率与波前
曲率用来约束地下速度模型,缓解了层析?#19978;?#30340;射线多路?#26029;?#35937;,具有更高的精
度与反演稳定性。此外,本发明属于层析?#19978;?#30340;范畴,只需要实施一次数据拾取
操作,因?#31169;?#20559;移速度分析具有更高的计算效率。

虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明
而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人
?#20445;?#22312;不脱离本发明所公开的精神和?#27573;?#30340;前提下,可以在实施的形式上及?#38468;?br />上作任何的修改与变化,但本发明的专利保护?#27573;В?#20173;须以所附的权利要求书所
界定的?#27573;?#20026;准。

关于本文
本文标题:一种地下速度模型的构建方法.pdf
链接地址:http://www.pqiex.tw/p-5994791.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

[email protected] 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备17046363号-1 
 


收起
展开
平码五不中公式规律 山东十一选五任四遗漏 彩票中奖都是没人信 福建31选7开奖号码查询今天 3d缩水怎么杀组选号 篮球场地标准尺寸画法 棋牌注册送50元50提现 21点技巧安卓版 516棋牌平台 马会一尾中特 组六