华北中北部地区震源参数和场地响应的联合反演

利用遗传算法反演华北中北部地区27个地震震源谱的低频水平、拐角频率和25个台站的场地响应。具体方法是利用研究区的平均衰减模型对记录谱进行路径校正，然后假定每一个地震的ω^2模型，调整低频水平和拐角频率使所有台站场地响应的标准偏差最小，最后给出参与计算的每个台站的场地响应。结果表明，17个基岩台中11个台的场地响应在1～20Hz范围内接近1，6个基岩台的场地响应部分频段放大或缩小3倍以上。8个井下摆台站的场地响应形态相同，表现为低频放大，高频部分迅速衰减。地震矩与震级的关系为logM0=17．1 1．13ML，地震矩与震源半径的关系为logM0=19．4 0．0052r。     

基于改进遗传算法的BISQ模型多参数反演方法研究

BISQ模型反演在油气资源与海洋地震勘探中具有广阔的应用前景与实用价值.为解决具高维参数与多峰值的BISQ模型参数反演问题,本文在传统的遗传算法基础上添加插入了全局保优机制、Boltzmann生存机制、小生境之DC漂移技术和群体灾变等方法与技术,使改进组合的遗传算法能在高维参数与多峰值的反问题中搜索得到全局最优解.结合海底多频原位沉积物声学测量数据,应用本文方法,得到了与声波传播速度相似性最优模型的各参数值,效果明显.     

利用数字化地震记录进行震源参数和介质衰减特性的联合反演

利用区域台网的数字化地震记录资料，发展了一种体波谱的反演方法，以消除用波谱分析求解震源参数和介质衰减特性过程中的不确定因素的影响．经过数值实验和对实际资料的应用，初步证明了该方法的实用性。利用滦县地区流动台网中的五个台站接收的３５次地震记录的Ｐ波频谱资料，反演了滦县地区的震源动力学参数，得到了５个台站的Ｐ波表现Ｑ值和震源频谱的高频衰减率γ。Ｑ值分布在４０８－８４７之间，均值为５２０。γ值分布在１．５４－３．２２之间，均值为２．４１。利用唐山仪器台阵中的７个台站接收的９次地震记录的Ｓ波频谱资料，反演了唐山滦县地区的震源动力学参数，并得到了７个台站所在的小区域内平均的Ｓ波表观Ｑ值和震源频谱的高频衰减率γ。Ｑ值的均值为５８３。γ值均值为２．５８。     

基于数据与模型联合驱动的波阻抗反演方法

井插值初始模型为基于模型的反演提供的低频信息往往不够准确,导致该模型驱动方法容易出现较大的波阻抗预测误差且建模效率较低.为缓解这些问题,本文利用数据驱动的深度学习反演更加擅长预测低频阻抗的优势,提出一种基于数据与模型联合驱动的波阻抗反演方法.该方法联合地震和测井等数据,先后开展数据驱动和模型驱动的波阻抗反演.首先,数据驱动部分使用井旁地震记录、测井导出的波阻抗曲线以及井插值低频阻抗曲线,搭建以双向门控递归单元为主要模块的波阻抗智能预测网络.其次,该网络预测的波阻抗的低频分量作为数据驱动初始模型,替代井插值初始模型而参与模型驱动部分.最后,模型驱动部分在地震数据匹配和数据驱动初始模型的共同约束下开展基于模型的反演,获得最终的波阻抗结果.合成数据和实际数据测试表明,本文方法相比于单一的数据驱动或模型驱动方法能获得更高分辨率和更高精度的波阻抗反演结果,从而为后续储层预测提供可靠的弹性参数分布.

     

基于非饱和多孔隙介质BISQ模型的储层参数反演

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制．文中基于同时包含这两种力学机制的非饱和多孔隙BISQ模型,利用小生境遗传算法实现了储层参数(孔隙度、渗透率、含流体饱和度等)的反演．结果表明,本文方法在储层参数的反演过程中目标函数收敛快、且具有较强的抗“噪声”干扰性能,当观测噪声ε≤5%时,储层参数的反演精度很高．最后,通过实测数据的反演应用验证了该方法处理储层参数反演问题的有效性和非饱和BISQ模型的准确性．     

密度和速度随机分布共网格模型的重力与地震联合反演

针对重力与地震联合反演存在的问题,结合已有的研究成果,本文研究实现了速度和密度随机分布共网格单元模型的建模技术,以适应密度和速度剧烈变化的复杂模型及联合反演的计算要求.重力正演利用了该网格的二度半体模型,并进一步改进了地震走时的二维射线追踪计算方法,以适用于速度随机分布的网格介质.结合改进的模拟退火算法,实现了这种共网格条件下的重力与地震资料的同步联合反演.模型试验证明了重力与地震联合反演可以准确确定复杂物性界面的密度和速度结构,适用于物性界面不完全一致和物性变化剧烈的复杂模型,并且联合反演结果要优于单独的重力反演.带先验信息约束下的实际资料的联合反演,进一步证明了该方法的适用性和效果,可提高反演精度并减少多解性.     

一种利用频谱激电数据实部和虚部加权反演Cole-Cole模型参数的方法

频谱激电法(SIP)因具有多参数(Cole-Cole参数)解释异常的能力,一直受到电法工作者的关注.一些学者提出了用电场振幅和相位反演Cole-Cole参数的方法,但由于电场振幅和相位对不同Cole-Cole参数的敏感程度不同,往往导致有些Cole-Cole参数较难反演,特别是频率相关系数和时间常数.本文提出了利用视复电阻率的实部和虚部加权反演的方法.该方法采用视复电阻率实部和虚部联合反演的策略,并在反演不同Cole-Cole参数时对实部和虚部的目标函数进行不同的加权.模型试算结果表明:这种反演方法不仅可以很好的反演零频电阻率和极化率,还可以显著提高频率相关系数和时间常数的反演效果.     

邓肯-张模型参数反演的两种不适定问题

利用邓肯-张模型进行参数反演时,不同的参数反演策略会存在不适定问题。结合参数反演的函数图形分析,指出在邓肯-张模型参数反演的过程中,基于φ反演E_t,基于c反演E_t,基于c,φ反演E_t时都会存在两种不适定问题。基于函数导数分析方法,结合莫尔-库仑准则表达式,解释邓肯-张模型参数反演产生两种不适定问题的原因,并提出消除方法。结果表明,邓肯-张模型参数的表达式是产生两种不适定问题的根源。     

基于波动方程的微地震震源位置、震源时间及VTI介质各向异性参数联合反演

对于非常规油气开发,水力压裂监控的效果取决于对微地震事件的分析、解释.准确的微地震震源位置是关乎施工成败的重要因素.微地震震源位置的准确性与多个参数相关,其不仅依赖于微地震事件的激发时间,同时也依赖于储层介质参数信息,因此进行微地震震源位置、震源时间、储层介质参数的联合反演尤为重要.页岩气储层通常表现出较强的各向异性,...     

层状介质参数反演的混合最优化法

波动方程反演是典型的非线性反问题,本文提出了用混合最优化法反演层状介质参数。混合最优化法将广义模拟退火与局部最优化方法结合,能较好地利用两者的优点,在本文中,局部最优化算法采用线性迭代算法。通过数值结果表明反演方法的正确有效性。     

大地电磁与地震资料仿真退火约束联合反演

在对仿真退火算法研究的基础上,将非常快速仿真退火的连续模型空间改为可变的离散化形式,理论模型试算表明,使用本文方法可提高工作效率5-10倍. 为适应联合反演的需要, 对非常快速仿真退火的误差"单判断"检验准则改为误差"双判断"检验准则,避免Metrop olis接受准则的判断失误. 在此基础上,实现了一维及二维起伏地形条件下大地电磁与地震资料多参量信息的联合反演. 理论模型的试算及实例说明了约束联合反演的有效性.     

新丰江水库三维速度结构和震源参数的联合反演

本文对ACH方法做简化处理,先以一个地震为单元形成子矩阵,并将慢度扰动用一平均量来代替;在修正参数时,震源部分按常规进行,对慢度部分则利用代数重建技术将其分配到各个块体上.数值模拟结果表明,该方法能很好地反演震源位置和速度结构.运用新丰江遥测台网的数据进行研究,表明地震强活动区位于大坝附近,特别是在人字石断裂和高寨断裂的交汇处,震源集中沿罗坑—双下、葫芦凹—燕岩和黄竹蒿3条地震带,呈密集的高倾角分布,地震活动主要与北北西向构造和东西向构造有关,而与地表甚为发育的北东东向构造关系不大.速度图像中的低速区对应地表的破碎区,即地震带区域,也是该区重磁异常急剧变化的区域;而高速区内地震较少.上述特征同水的渗透作用有密切关系.     

新丰江水库三维速度结构和震源参数的联合反演

本文对ACH方法做简化处理,先以一个地震为单元形成子矩阵,并将慢度扰动用一平均量来代替;在修正参数时,震源部分按常规进行,对慢度部分则利用代数重建技术将其分配到各个块体上.数值模拟结果表明,该方法能很好地反演震源位置和速度结构.运用新丰江遥测台网的数据进行研究,表明地震强活动区位于大坝附近,特别是在人字石断裂和高寨断裂的交汇处,震源集中沿罗坑-双下、葫芦凹-燕岩和黄竹蒿3条地震带,呈密集的高倾角分布,地震活动主要与北北西向构造和东西向构造有关,而与地表甚为发育的北东东向构造关系不大.速度图像中的低速区对应地表的破碎区,即地震带区域,也是该区重磁异常急剧变化的区域;而高速区内地震较少.上述特征同水的渗透作用有密切关系.     

大地电磁与地震资料仿真退火约束联合反演

在对仿真退火算法研究的基础上,将非常快速仿真退火的连续模型空间改为可变的离散化形式,理论模型试算表明,使用本文方法可提高工作效率5-10倍. 为适应联合反演的需要, 对非常快速仿真退火的误差"单判断"检验准则改为误差"双判断"检验准则,避免Metrop olis接受准则的判断失误. 在此基础上,实现了一维及二维起伏地形条件下大地电磁与地震资料多参量信息的联合反演. 理论模型的试算及实例说明了约束联合反演的有效性.     

电偶源频率电磁测深视电阻率多参量联合反演效果检验

通过大量的数值模拟试验 ,对电偶源频率电磁测深 (FEMS)视电阻率多参量联合反演的效果进行了检验和对比 .结果表明 ,多参量联合反演具有较快的收敛速度和较高的拟合精度     

二维反射地震成象—界面与速度的联合反演

介绍了利用反射走时联合反演地震界面起伏形态和速度横向分布的方法。正问题采用了一种简单而有效的二维介质中两点射线追踪方法，反问题在广义最小二乘意义下求解，描述界面和速度的两类参数不必要先进行分离，数值模拟的结果表明方法是有效的。处理了通过唐山震区的各庄－丰南－丰宁深地震测深剖面，得到了意义的结果。     

用双侧向测井响应求取地层参数的Newton－SVD反演方法

针对测井中经常遇到的完全非均匀的地层模型，利用深、浅双侧向测井的视电阻率曲线作为约束条件，给出了反演侵入半径ｒ_ｉ、原状地层电阻率Ｒ_ｔ、上（下）围岩电阻率Ｒ_{ｓｕ（Ｒｓｄ）的Ｎｅｗｔｏｎ－ＳＶＤ反演方法．数值模拟表明这一方法是可行的，经实际资料处理．说明该方法实用．}     

利用数字地震台网资料联合反演福建地区Q值、场地响应和震源参数

文中利用福建数字地震台网1998年建网以来记录的2.5级以上地震的波形资料,采用Atkinson等提出的方法计算了福建地区地震波非弹性衰减系数、几何扩散系数和台站场地响应。考虑到该方法在计算台站场地响应时得到的结果是相对场地响应,利用Moya等提出的方法对台站场地响应进行重新计算,对震源谱的低频水平和拐角频率进行联合反演,并计算了参与反演的地震的新参数(地震矩、应力降和破裂半径),同时讨论了它们之间的关系。反演结果表明福建地区地壳平均Q值与频率的关系为Q=504.1×f0.332,而场地响应均表现为与频率相关并基本在1附近波动,这与福建数字地震台站基本建立在基岩之上的实际情况相符合。震源参数的计算结果表明,地震矩与震级以及震源谱拐角频率的相关性较好,而地震矩与应力降以及应力降与震源半径之间的依赖关系不很明确,这可能与参与计算的地震数量以及震级范围有关     

大地电磁二维层状模型参数化反演

本文针对目前二维反演存在求解参数过多和求其偏导数困难的两个问题，在近年来关于连续介质二维大地电磁测深资料快速反演的理论基础上，对二维层状模型反演的目标函数，进行线性近似，采用加权法来近似计算二维模型参数的偏导数。从而实现比较快速的二维层状模型的参数化反演。理论模型的反演试验表明，反演迭代可以较快地以足够的精度拟合视电阻率和相位资料，并收敛接近理论模型。