CN109061731A - 面波频散与体波谱比联合反演浅层速度的全局优化方法 - Google Patents
面波频散与体波谱比联合反演浅层速度的全局优化方法 Download PDFInfo
- Publication number
- CN109061731A CN109061731A CN201811087270.5A CN201811087270A CN109061731A CN 109061731 A CN109061731 A CN 109061731A CN 201811087270 A CN201811087270 A CN 201811087270A CN 109061731 A CN109061731 A CN 109061731A
- Authority
- CN
- China
- Prior art keywords
- wave
- shallow
- curve
- bulk
- population
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 230000002068 genetic effect Effects 0.000 claims abstract description 29
- 238000002922 simulated annealing Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 239000002689 soil Substances 0.000 claims description 26
- 239000011435 rock Substances 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 229910002056 binary alloy Inorganic materials 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000000137 annealing Methods 0.000 claims description 3
- 238000001816 cooling Methods 0.000 claims description 3
- 238000013517 stratification Methods 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 230000035772 mutation Effects 0.000 claims description 2
- 230000008521 reorganization Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims 1
- 238000011160 research Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect 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
本发明涉及一种面波频散与体波谱比联合反演浅层速度的全局优化方法,包括以下步骤:1)给定反演的初始参数,构造以二进制编码表示的未知向量种群;2)确定面波频散曲线的理论计算方法和体波HVSR曲线的理论计算方法;3)建立能同时拟合面波频散曲线和体波HVSR曲线的目标函数;4)进行遗传模拟退火的全局优化计算,获得全局最优解;5)根据全局最优解即可获得的浅层速度结构。本发明结合了遗传和模拟退火两种方法的优势,可以提高反演的效率和准确性,可以广泛应用于测线难以展开的城市地区,也能用于背景噪声、测震或强震记录等多种数据源。
Description
技术领域
本发明涉及一种浅层地震勘探领域中的近地表浅层速度结构的获取方法,特别是关于一种基于面波频散与体波谱比联合反演浅层速度的全局优化方法。
背景技术
浅层速度结构是开展场地地震影响分析并评价场地地震效应的基础资料。获取浅层速度结构的常规方法有钻孔测井法、主动源探测法等。钻孔测井受到费用、钻探深度限制,对一些空间展布范围广的工程只能在关键控制点布设,主动源探测则受场地条件限制较大,对土层结构的分辨能力也难以满足工程需要。因此,研究和开发方便、经济的浅层速度结构探查方法,一直是地震工程以及浅层地震勘探领域的重要内容。
目前,地震学利用实际观测记录数据的反演仍多用于深部结构的研究,用于反演地震工程较为关注的局部场地土层(一般为百米或数百米深度的浅地层)速度的技术手段并不多见。单独利用面波频散曲线反演地下速度结构已有比较成熟的方法,但这种方法多解性较强,需要增加有效的约束条件。近年来,从噪声水平与竖向谱比方法扩展而来的水平与竖向谱比(Horizontal-to-Vertical Spectral Ratio,以下简称HVSR)也是一种研究场地固有特性的有效方法,其仅需单个台的三分量地表观测记录,数据获取容易且成本较低,引起了国内外研究者的重视与持续研究。已有的研究表明,体波HVSR曲线仅与场地S波、P波传递函数以及基岩S波、P波波速等场地本身的特性有关,因而可用于反演土层速度结构,但体波HVSR方法用于反演浅层速度结构时同样面临多解性较强的问题。由于面波频散曲线和体波HVSR曲线都能体现地层速度结构特征,且二者原理不同,若互为补充可以非常有效地增加反演时的约束条件,有效解决反演多解性问题。另外,场地浅层速度结构的反演计算实质是一个土层场地模型空间搜索的全局优化问题,这一问题的解决较为依赖反演方法的选择,当模型搜索空间的复杂程度增大时,目前常用的搜索算法如蒙特卡洛方法收敛速度变慢,计算效率降低。
相比较而言,遗传算法和模拟退火算法是更智能的自适应全局优化方法,其概率搜索高效,在波速结构反演中得到广泛应用。遗传算法虽然具有隐并行性,全局搜索能力强,但容易收敛到局部最优解,即收敛早熟,模拟退火算法局部搜索能力强,但全局搜索能力却较差,若能将二者结合起来,就能取长补短极大提高反演的效率与准确性。
如前所述,反演浅层速度结构的方法大多采用单独的面波频散曲线反演或者单独的体波HVSR反演,联合应用面波频散曲线和体波HVSR曲线反演的研究尚十分少见。在联合反演方法方面,现有的多为单独的遗传算法或者模拟退火算法等,遗传算法存在收敛早熟、模拟退火算法存在搜索效率低的问题。随着我国浅地层探测技术的不断进步和发展,钻孔测井方法费用较高、主动源探测在城市地区难以开展的缺点已经日益突出,急需寻求方便、经济、稳定可行的浅层速度结构探查方法。传统方法已经难以适应大面积地区以及城市地区浅层速度结构探测的要求。
发明内容
针对上述问题,本发明的目的是提供一种面波频散曲线与体波HVSR曲线联合反演浅层速度结构的全局混合优化方法,该方法结合了遗传和模拟退火两种方法的优势,可以提高反演的效率和准确性。
为实现上述目的,本发明采用以下技术方案:一种面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,包括以下步骤:1)给定反演的初始参数,构造以二进制编码表示的未知向量种群;2)确定面波频散曲线的理论计算方法和体波HVSR曲线的理论计算方法;3)建立能同时拟合面波频散曲线和体波HVSR曲线的目标函数;4)进行遗传模拟退火的全局优化计算,获得全局最优解;5)根据全局最优解即可获得的浅层速度结构。
优选的,所述步骤1)具体包括:给定的初始参数包括根据反演目标场地的有限信息得到的浅地层的各土层密度、P波和S波品质因子、P波和S波波速之间的经验关系,将S波波速和土层厚度作为未知量来构造初始种群,初始种群以二进制进行编码以便于后续遗传模拟退火方法的种群进化。
优选的,所述步骤2)具体包括:面波频散曲线的理论计算公式如下:
其中,ω为地震波的角频率;VR为地震波的相速度;n为岩土体的层数;Vpn和Vsn分别为第n层岩土体的P波波速和S波波速;V、W、R和S分别为岩土体介质特性相关系数;
体波HVSR曲线的理论计算公式如下:
其中,TFS(f)和TFp(f)分别为S波和P波传递函数;α和β分别为基岩处P波波速和S波波速。
优选的,所述步骤3)中,目标函数的表达式如下:
其中,X={x1,x2,x3,…,xj}为地下土层结构分层模型的参数向量,j=1,2,3,…,n;xj={Vsj,Vpj,hj,μj,ρj,Qsj,Qpj},Vsj,Vpj,hj,μj,ρj,Qsj,Qpj分别为第j层岩土体的S波波速、P波波速、厚度、泊松比、密度、S波品质因子、P波品质因子;f为地震波的频率;DCt(f)和DCo(f)分别为理论和观测面波频散曲线;HVSRt(f)和HVSRo(f)分别为理论和观测S波HVSR曲线;A为面波频散曲线最大值与体波HVSR曲线最大值的比值。
优选的,所述步骤4)具体包括:⑴利用式(3)计算初始种群的目标函数;⑵设定模拟退火初始温度、最低温度、最大遗传代数和降温比例系数,同时设定的目标函数阈值;⑶若目标函数最小值小于设定的目标函数阈值且模拟退火计算温度大于设定的模拟退火初始温度,则根据计算的初始种群的目标函数分配适应度值;⑷进行遗传算法的选择操作,即根据选择适应度值较高的个体进入下一代,淘汰适应度低的个体;⑸进行遗传算法的重组操作,即将群体中个体的二进制码串随机两两配对并随机确定一个码串的交叉位置,互换配对的两个个体在交叉位置之后的二进制串,重组成新的群体;⑹进行遗传算法的变异操作,随机确定新的群体中个体的变异位置,按给定的概率将变异位置上的二进制码进行0、1互换,形成子代群体;⑺利用式(3)计算子代群体的目标函数;⑻进行模拟退火操作,如果子代个体的目标函数小于父代个体,则保留该子代个体;如果子代个体目标函数大于父代个体,则以Boltzman接受概率接受该子代个体;⑼生成经过模拟退火操作后新的子代种群;⑽循环进行⑷~⑼步骤,直到遗传代数达到设定的最大遗传代数;⑾循环进行⑶~⑽步骤,直到模拟退火计算温度达到设定的最低温度且目标函数达到设定的最小值,输出S波波速和土层厚度的全局最优的解。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于同时使用了地表观测得到的面波与体波信息,比常规的单独利用面波频散曲线的反演和单独利用体波HVSR的反演方法增加了更多的约束条件,能够更加准确地获得浅地层速度结构。2、本发明方法采用的面波频散曲线和体波HVSR曲线计算方法均采用最新的研究成果,计算精度高,克服了常规方法频散曲线和HVSR曲线计算不准确的缺点。3、本发明方法采用了遗传模拟退火混合全局优化方法,计算速度快、精度高,有效避免了遗传算法的早熟问题和模拟退火搜索效率低的缺点。4、本发明是一种低成本的浅地层速度结构探查方法,可以广泛应用于测线难以展开的城市地区,也能用于背景噪声、测震或强震记录等多种数据源。
附图说明
图1是本发明的流程示意图;
图2是本发明反演获得的模型计算得到的面波频散曲线与目标频散曲线的比较;
图3是本发明反演获得的模型计算得到的体波HVSR曲线与目标HVSR曲线的比较;
图4是本发明反演获得的浅层速度结构与目标速度结构的比较。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提出了一种面波频散与体波谱比联合反演浅层速度的全局优化方法,包括以下步骤:
1)给定反演的初始参数,构造以二进制编码表示的未知向量种群。
给定的初始参数包括根据反演目标场地的有限信息得到的浅地层的各土层密度、P波和S波品质因子、P波和S波波速之间的经验关系,将S波波速和土层厚度作为未知量来构造初始种群,初始种群以二进制进行编码以便于后续遗传模拟退火方法的种群进化。
2)确定面波频散曲线的理论计算方法和体波HVSR曲线的理论计算方法。
面波频散曲线的理论计算采用Knopoff改进的Thomson-Haskell方法,该方法采用地下均匀水平固体介质分层模型考虑角频率为ω、相速度为VR的地震波在n层岩土体中的传播,面波频散曲线的理论计算公式如下:
其中,Vpn和Vsn分别为第n层岩土体的P波波速和S波波速;V、W、R和S分别为岩土体介质特性相关系数。
体波HVSR曲线的理论计算公式如下:
其中,TFS(f)和TFp(f)分别为S波和P波传递函数;α和β分别为基岩处P波波速和S波波速。
3)建立能同时拟合面波频散曲线和体波HVSR曲线的目标函数。
目标函数的表达式如下:
其中,X={x1,x2,x3,…,xj}为地下土层结构分层模型的参数向量,j=1,2,3,…,n;xj={Vsj,Vpj,hj,μj,ρj,Qsj,Qpj},Vsj,Vpj,hj,μj,ρj,Qsj,Qpj分别为第j层岩土体的S波波速、P波波速、厚度、泊松比、密度、S波品质因子、P波品质因子;f为地震波频率;DCt(f)和DCo(f)分别为理论和观测面波频散曲线;HVSRt(f)和HVSRo(f)分别为理论和观测S波HVSR曲线;A为面波频散曲线最大值与体波HVSR曲线最大值的比值。当n较大时,Φ(X)为多参数、多极值的非线性函数。
4)进行遗传模拟退火的全局优化计算,具体计算过程如下:
⑴利用式(3)计算初始种群的目标函数;
⑵设定模拟退火初始温度、最低温度、最大遗传代数和降温比例系数,同时设定的目标函数阈值;
⑶若目标函数最小值小于设定的目标函数阈值且模拟退火计算温度大于设定的模拟退火初始温度,则根据计算的初始种群的目标函数分配适应度值;
⑷进行遗传算法的选择操作,即根据选择适应度值较高的个体进入下一代,淘汰适应度低的个体;
⑸进行遗传算法的重组操作,即将群体中个体的二进制码串随机两两配对并随机确定一个码串的交叉位置,互换配对的两个个体在交叉位置之后的二进制串,重组成新的群体;
⑹进行遗传算法的变异操作,随机确定新的群体中个体的变异位置,按给定的概率将变异位置上的二进制码进行0、1互换,形成子代群体;
⑺利用式(3)计算子代群体的目标函数;
⑻进行模拟退火操作,如果子代个体的目标函数小于父代个体,则保留该子代个体;如果子代个体目标函数大于父代个体,则以Boltzman接受概率接受该子代个体;
⑼生成经过模拟退火操作后新的子代种群;
⑽循环进行⑷~⑼步骤,直到遗传代数达到设定的最大遗传代数;
⑾循环进行⑶~⑽步骤,直到模拟退火计算温度达到设定的最低温度且目标函数达到设定的最小值,输出S波波速和土层厚度的全局最优的解。
5)根据S波波速和土层厚度的全局最优解即可获得的浅层速度结构。
图2和图3分别显示了利用本发明反演获得的模型计算得到的面波频散曲线与目标频散曲线的比较,以及反演获得的模型计算得到的体波HVSR曲线与目标HVSR曲线的比较,图4给出了本发明反演获得的浅层速度结构与目标速度结构的比较。这进一步说明本发明方法可以根据观测的面波频散和体波HVSR较为准确地得到地下浅层速度结构,且该浅层速度结构得到的理论面波频散和体波HVSR与观测吻合较好。从图中可以看出,采用本发明反演策略获得的浅层速度结构能有效避免反演的多解性问题。
上述各实施例仅用于说明本发明,其中的一些方法和步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,包括以下步骤:
1)给定反演的初始参数,构造以二进制编码表示的未知向量种群;
2)确定面波频散曲线的理论计算方法和体波HVSR曲线的理论计算方法;
3)建立能同时拟合面波频散曲线和体波HVSR曲线的目标函数;
4)进行遗传模拟退火的全局优化计算,获得全局最优解;
5)根据全局最优解即可获得的浅层速度结构。
2.根据权利要求1所述的面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,所述步骤1)具体包括:
给定的初始参数包括根据反演目标场地的有限信息得到的浅地层的各土层密度、P波和S波品质因子、P波和S波波速之间的经验关系,将S波波速和土层厚度作为未知量来构造初始种群,初始种群以二进制进行编码以便于后续遗传模拟退火方法的种群进化。
3.根据权利要求1所述的面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,所述步骤2)具体包括:
面波频散曲线的理论计算公式如下:
其中,ω为地震波的角频率;VR为地震波的相速度;n为岩土体的层数;Vpn和Vsn分别为第n层岩土体的P波波速和S波波速;V、W、R和S分别为岩土体介质特性相关系数;
体波HVSR曲线的理论计算公式如下:
其中,TFS(f)和TFp(f)分别为S波和P波传递函数;α和β分别为基岩处P波波速和S波波速。
4.根据权利要求1所述的面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,所述步骤3)中,标函数的表达式如下:
其中,X={x1,x2,x3,…,xj}为地下土层结构分层模型的参数向量,j=1,2,3,…,n;xj={Vsj,Vpj,hj,μj,ρj,Qsj,Qpj},Vsj,Vpj,hj,μj,ρj,Qsj,Qpj分别为第j层岩土体的S波波速、P波波速、厚度、泊松比、密度、S波品质因子、P波品质因子;f为地震波的频率;DCt(f)和DCo(f)分别为理论和观测面波频散曲线;HVSRt(f)和HVSRo(f)分别为理论和观测S波HVSR曲线;A为面波频散曲线最大值与体波HVSR曲线最大值的比值。
5.根据权利要求1所述的面波频散与体波谱比联合反演浅层速度的全局优化方法,其特征在于,所述步骤4)具体包括:
⑴利用式(3)计算初始种群的目标函数;
⑵设定模拟退火初始温度、最低温度、最大遗传代数和降温比例系数,同时设定的目标函数阈值;
⑶若目标函数最小值小于设定的目标函数阈值且模拟退火计算温度大于设定的模拟退火初始温度,则根据计算的初始种群的目标函数分配适应度值;
⑷进行遗传算法的选择操作,即根据选择适应度值较高的个体进入下一代,淘汰适应度低的个体;
⑸进行遗传算法的重组操作,即将群体中个体的二进制码串随机两两配对并随机确定一个码串的交叉位置,互换配对的两个个体在交叉位置之后的二进制串,重组成新的群体;
⑹进行遗传算法的变异操作,随机确定新的群体中个体的变异位置,按给定的概率将变异位置上的二进制码进行0、1互换,形成子代群体;
⑺利用式(3)计算子代群体的目标函数;
⑻进行模拟退火操作,如果子代个体的目标函数小于父代个体,则保留该子代个体;如果子代个体目标函数大于父代个体,则以Boltzman接受概率接受该子代个体;
⑼生成经过模拟退火操作后新的子代种群;
⑽循环进行⑷~⑼步骤,直到遗传代数达到设定的最大遗传代数;
⑾循环进行⑶~⑽步骤,直到模拟退火计算温度达到设定的最低温度且目标函数达到设定的最小值,输出S波波速和土层厚度的全局最优的解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811087270.5A CN109061731B (zh) | 2018-09-18 | 2018-09-18 | 面波频散与体波谱比联合反演浅层速度的全局优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811087270.5A CN109061731B (zh) | 2018-09-18 | 2018-09-18 | 面波频散与体波谱比联合反演浅层速度的全局优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109061731A true CN109061731A (zh) | 2018-12-21 |
CN109061731B CN109061731B (zh) | 2019-10-25 |
Family
ID=64762913
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811087270.5A Active CN109061731B (zh) | 2018-09-18 | 2018-09-18 | 面波频散与体波谱比联合反演浅层速度的全局优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109061731B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109799530A (zh) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | 用于地震面波勘探的瑞雷波频散曲线反演方法 |
CN110426740A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110426739A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探检测方法、装置和存储介质 |
CN110426741A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110568484A (zh) * | 2019-08-02 | 2019-12-13 | 中铁第四勘察设计院集团有限公司 | 一种反演方法、装置及存储介质 |
CN110879412A (zh) * | 2019-10-31 | 2020-03-13 | 南方科技大学 | 地下横波速度反演方法、装置、计算设备及存储介质 |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
CN112327364A (zh) * | 2020-11-02 | 2021-02-05 | 中国石油大学(华东) | 一种改进型一步法三维crs叠加方法 |
CN112698407A (zh) * | 2020-11-23 | 2021-04-23 | 中国石油天然气集团有限公司 | 一种快速反演声波测井弯曲波频散曲线的方法 |
CN112731558A (zh) * | 2020-12-16 | 2021-04-30 | 中国科学技术大学 | 一种地震面波与接收函数的联合反演方法及装置 |
CN112904407A (zh) * | 2020-10-14 | 2021-06-04 | 东华理工大学 | 复杂地形与干扰条件下的微动勘探方法 |
CN112987090A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
WO2021227236A1 (zh) * | 2020-05-12 | 2021-11-18 | 山东大学 | 基于盾构机施工噪声的多波场地震探测方法与系统 |
CN116520402A (zh) * | 2023-06-02 | 2023-08-01 | 大连理工大学 | 岩石半空间场地下考虑体波与面波的多震相波场反演方法 |
CN116660974A (zh) * | 2023-04-11 | 2023-08-29 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
CN116819622A (zh) * | 2023-08-30 | 2023-09-29 | 北京工业大学 | 土层三维速度结构的背景噪声水平竖向谱比联合反演方法 |
CN118039022A (zh) * | 2024-04-11 | 2024-05-14 | 国能龙源环保有限公司 | 二氧化硫浓度预测方法、装置及计算机设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
-
2018
- 2018-09-18 CN CN201811087270.5A patent/CN109061731B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
Non-Patent Citations (3)
Title |
---|
崔建文: "一种改进的全局优化算法及其在面波频散曲线反演中的应用", 《地球物理学报》 * |
胥鸿睿: "鄂尔多斯块体东缘横波速度结构及各向异性研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
荣棉水: "基于单台加速度纪录的混合全局优化HVSR反演场地浅层速度结构", 《地球物理学报》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109799530A (zh) * | 2018-12-25 | 2019-05-24 | 核工业北京地质研究院 | 用于地震面波勘探的瑞雷波频散曲线反演方法 |
CN110426741A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110426739A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探检测方法、装置和存储介质 |
CN110568484B (zh) * | 2019-08-02 | 2021-07-16 | 中铁第四勘察设计院集团有限公司 | 一种反演方法、装置及存储介质 |
CN110568484A (zh) * | 2019-08-02 | 2019-12-13 | 中铁第四勘察设计院集团有限公司 | 一种反演方法、装置及存储介质 |
CN110426740A (zh) * | 2019-08-02 | 2019-11-08 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110426740B (zh) * | 2019-08-02 | 2022-02-15 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110426741B (zh) * | 2019-08-02 | 2021-07-16 | 中铁第四勘察设计院集团有限公司 | 一种地震噪音成像勘探方法、装置和存储介质 |
CN110879412A (zh) * | 2019-10-31 | 2020-03-13 | 南方科技大学 | 地下横波速度反演方法、装置、计算设备及存储介质 |
CN112987090B (zh) * | 2019-12-02 | 2024-05-28 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
CN112987090A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油天然气集团有限公司 | 面波频散曲线拾取方法及装置 |
US11644590B2 (en) | 2020-05-12 | 2023-05-09 | Shandong University | Multi-wavefield seismic detection method and system based on construction noise of shield machine |
WO2021227236A1 (zh) * | 2020-05-12 | 2021-11-18 | 山东大学 | 基于盾构机施工噪声的多波场地震探测方法与系统 |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
CN112904407A (zh) * | 2020-10-14 | 2021-06-04 | 东华理工大学 | 复杂地形与干扰条件下的微动勘探方法 |
CN112904407B (zh) * | 2020-10-14 | 2023-05-30 | 东华理工大学 | 复杂地形与干扰条件下的微动勘探方法 |
CN112327364A (zh) * | 2020-11-02 | 2021-02-05 | 中国石油大学(华东) | 一种改进型一步法三维crs叠加方法 |
CN112698407A (zh) * | 2020-11-23 | 2021-04-23 | 中国石油天然气集团有限公司 | 一种快速反演声波测井弯曲波频散曲线的方法 |
CN112731558A (zh) * | 2020-12-16 | 2021-04-30 | 中国科学技术大学 | 一种地震面波与接收函数的联合反演方法及装置 |
CN116660974B (zh) * | 2023-04-11 | 2024-02-23 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
CN116660974A (zh) * | 2023-04-11 | 2023-08-29 | 中国地震局地球物理研究所 | 一种基于结构耦合约束的体波和面波三维联合反演方法 |
CN116520402B (zh) * | 2023-06-02 | 2023-11-28 | 大连理工大学 | 岩石半空间场地下考虑体波与面波的多震相波场反演方法 |
CN116520402A (zh) * | 2023-06-02 | 2023-08-01 | 大连理工大学 | 岩石半空间场地下考虑体波与面波的多震相波场反演方法 |
CN116819622A (zh) * | 2023-08-30 | 2023-09-29 | 北京工业大学 | 土层三维速度结构的背景噪声水平竖向谱比联合反演方法 |
CN116819622B (zh) * | 2023-08-30 | 2023-11-21 | 北京工业大学 | 土层三维速度结构的背景噪声水平竖向谱比联合反演方法 |
CN118039022A (zh) * | 2024-04-11 | 2024-05-14 | 国能龙源环保有限公司 | 二氧化硫浓度预测方法、装置及计算机设备 |
CN118039022B (zh) * | 2024-04-11 | 2024-07-23 | 国能龙源环保有限公司 | 二氧化硫浓度预测方法、装置及计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN109061731B (zh) | 2019-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061731B (zh) | 面波频散与体波谱比联合反演浅层速度的全局优化方法 | |
Zhang et al. | Parameter prediction of hydraulic fracture for tight reservoir based on micro-seismic and history matching | |
Guyaguler et al. | Optimization of well placement in a Gulf of Mexico waterflooding project | |
CN104615862B (zh) | 基于进化算法的高含水油田确定井位的方法 | |
CN102880903B (zh) | 混沌量子粒子群模糊神经网络叠前非线性流体识别方法 | |
CN109543237A (zh) | 基于ga-bp神经网络的基坑位移预测方法 | |
Kalogirou et al. | Artificial neural networks for the generation of geothermal maps of ground temperature at various depths by considering land configuration | |
CN105093278B (zh) | 基于激发主能量优化算法的全波形反演梯度算子提取方法 | |
Clevis et al. | Geoarchaeological simulation of meandering river deposits and settlement distributions: A three‐dimensional approach | |
CN113610945B (zh) | 一种基于混合神经网络的地应力曲线预测方法 | |
CN109343125A (zh) | 一种基于探地雷达的红壤关键带地下结构空间预测方法 | |
CN104422969A (zh) | 一种减小电磁测深反演结果非唯一性的方法 | |
CN105242328A (zh) | 古热岩石圈厚度的确定方法及装置 | |
Mahmod et al. | Analysis of groundwater flow in arid areas with limited hydrogeological data using the Grey Model: a case study of the Nubian Sandstone, Kharga Oasis, Egypt | |
Yadav et al. | Feedforward neural network for joint inversion of geophysical data to identify geothermal sweet spots in Gandhar, Gujarat, India | |
CN102288996A (zh) | 叠前afi反演定量预测含气概率的方法 | |
CN107939372A (zh) | 针对小断块油藏的最优井位部署方法和装置 | |
Wang et al. | Use of retrospective optimization for placement of oil wells under uncertainty | |
Bekhit | Sustainable groundwater management in coastal aquifer of Sinai using evolutionary algorithms | |
Yuhao et al. | Research on Prediction of Ground Settlement of Deep Foundation Pit Based on Improved PSO-BP Neural Network | |
CN113323618A (zh) | 一种特低渗透油藏人工裂缝水窜治理方法及设备 | |
CN106595798B (zh) | 一种使用高频地下水位数据自动求取含水层参数的方法 | |
Li et al. | Back-analysis of pavement thickness based on PSO-GA hybrid algorithms | |
CN112507615A (zh) | 一种陆相致密储层岩相智能识别与可视化方法 | |
CN111751886A (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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Rong Mianshui Inventor after: Tang Hui Inventor before: Rong Mianshui |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |