Pg回折波的上地壳三维层析成象

提出利用人工地震Ｐｇ回折波走时重建区域上地壳三维速度分布的方法。正问题利用经典最小二乘迭代求解垂向速度梯度模型的参数值及相应走时残差。反问题采用泛函空间的最小二乘准则来实现。对速度不需要预先进行参数化处理，可以计算反演区内任意点的速度值。数值计算结果表明了该方法的可靠性。处理了滦县地震区人工地震Ｐｇ走时资料，获得了该区上地壳顶部的三维速度分布。     

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

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

三维复杂层状介质中多震相走时联合反演成像

采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;另外,考虑到同时反演速度模型和反射界面起伏中不同参数变化对...     

横向均匀介质中地震走时的正演和反演

横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态

研究利用反射波的走时反演三维弯曲界面和介质层速度的计算方法．各层界面利用分段非完全三次多项式描述．正问题采用一种快速的三维射线追踪方法．反演过程采用变阻尼最小二乘法．数值模拟结果表明，解很快收敛到真模型．处理了通过唐山震区的实测资料，重建了该震区莫霍界面三维构造形态，并揭示了该区域构造与地震活动的关系．      

联合利用走时与波形反演技术研究地壳三维速度结构(Ⅱ)——应用

使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演．所用资料为Ｓ波和Ｐ波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件．为建立初始模型,又利用天然地震构成了准二维剖面．在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正．以３４°～４２°Ｎ,９４°～１１２°Ｅ作为研究区域,在该区域中收集了１９８６年以来大量地震的Ｓ波和Ｐ波到时差资料,７条人工地震二维速度剖面资料和２个数字化地震台的几个地震的三分向记录资料．对这些资料进行了处理,最后得出了０～２５ｋｍ深度不同截面的速度分布,并对所得结果进行了分析．     

天水台影区地震震相特征讨论

利用天水地震台1986～1989年的地震资料及《中国地震台站临时报告》,试算出了通道波πg、Lg1、Lg2及其速度,用最小二乘法初步求得πg、Lg1、Lg2波在本区域的走时关系式,并对πg、Lg1、Lg2波震相在该台的记录特征做了讨论。     

二维无限频率初至走时反演井间地震实际资料

使用Zelt和Barton的方法,通过一个计算效率高的有限差分求解eikonal方程,正演计算走时和射线路径.使用最小二乘QR分解法,求解稀疏线性系统方程组.使用正则化层析反演,结合用户给定的最小的、最平坦和最平滑的扰动限制,每一个加权因子随深度变化.结合数据残差和模型粗糙度的最小化,为数据残差提供一个最平滑的近似模型.该反演方法为非线性反演,需要一个初始模型,在每一次迭代时,需要计算新的射线路径.使用二维初至走时数据,对某油田二维井间地震实际资料进行无限频率初至走时层析反演.将反演所得到的速度与井的测井速度曲线相比较,二者吻合程度较高,表明该反演方法所得速度的分辨率比较高.证实了二维无限频率初至走时层析反演可以为全波形反演提供一个分辨率较高的长波长速度模型,从而为全波形反演井间地震实际资料提供了一个比较可靠的初始速度模型.     

井间地震层析成象

多数井间地震走时层析成象法是将震源到接收器间的射线路径近似为直线,但这只有当研究区域速度变化较小时才有效.实际上许多地区都可观测到速度变化达到10—20%以上,引起显著的射线弯曲.虽然在研究井间成象问题时有人已考虑到这一非线性效应.但多数结果并不令人满意.本文提出一种二维射线追踪的迭代反演方法并成功地应用于野外实测数据. 本方法是先进行迭代射线追踪然后来修改速度模型.每次迭代都把在当前模型中经射线追踪计算得到的理论走时与观测走时之差通过线性方程组与未知的速度扰动联系起来,并通过阻尼最小二乘法求出用来修改模型的速度扰动.迭代一直进行到理论走时与观测资料拟合到给定误差,或走时拟合无法再改进时为止.理论模拟表明本方法经几次迭代后就可收敛到正确解. 利用这一方法处理了结晶岩地区的实测井间地震实验数据.地震波的频率范围是1—6.6 kHz,允许的分辨率为几米.得到的速度图象与其它地质和地球物理资料十分吻合.根据求解的速度模型计算出的马斯洛夫(Maslov)理论地震图与实际波形图相当一致,这也有力地验证了本反演方法的有效性.     

黄土高原及邻区地壳P波速度结构

利用地震波走时联合反演算法(改进型最短路径算法)进行三维弯曲地震射线追踪正演,以及共轭梯度法求解带约束的阻尼最小二乘问题进行反演,同时更新速度模型和地震震中位置,结合地方震和区域地震走时资料得到了黄土高原(含汾渭断陷盆地)及邻区地壳三维P波速度结构.其横向变化结果表明,研究区地壳内的P波高速异常区与其内的地震活动构造带相一致,地震多发生在P波高速异常区的边缘或高、低速异常区的交汇处.秦岭山区和鄂尔多斯块体东南区为P波低速异常区.而垂向变化结果则表明研究区存在低速异常区.     

黄土高原及邻区地壳P波速度结构

利用地震波走时联合反演算法(改进型最短路径算法)进行三维弯曲地震射线追踪正演, 以及共轭梯度法求解带约束的阻尼最小二乘问题进行反演, 同时更新速度模型和地震震中位置, 结合地方震和区域地震走时资料得到了黄土高原(含汾渭断陷盆地)及邻区地壳三维P波速度结构. 其横向变化结果表明, 研究区地壳内的P波高速异常区与其内的地震活动构造带相一致, 地震多发生在P波高速异常区的边缘或高、 低速异常区的交汇处. 秦岭山区和鄂尔多斯块体东南区为P波低速异常区. 而垂向变化结果则表明研究区存在低速异常区.      

辽宁地区近震走时与震中距间关系研究

本文通过对辽宁数字地震台网控制区内近11年的地震观测数据进行处理分析,应用最小二乘法求回归方程等一系列相关公式进行计算,最终依据辽宁地区近震走时与震中距之间关系,归纳成表。为编制辽宁地区近震走时表提供可靠的数据依据。     

基于多震相联合反演江苏地区一维速度模型研究

利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。     

联合利用走时与波形反演技术研究地壳三维速度结构(Ⅰ)——理论与方法

论述地震波走时反演技术中的阻尼最小二乘方法。为了提高解的唯一性,我们把人工地震资料的解释结果作为该方法的特定约束条件。同时还论述了波形反演技术中的遗传算法．     

利用射线追踪作井间层析成象

许多井间地震走时的层析成象解释方法是将震源和接收器之间的射线路径近似为直线,而不管射线路径上速度结构的不均匀性效应。当速度没有多大变化时,这样的近似是有效的。在许多地区观察到的速度横向变化为10—20％,甚至更多,这引起了射线相当大的弯曲。有些工作已经考虑了这种非线性效应,但在应用于井间地震问题时,还没有出现许多成功的实例。在这里我们提出了一个基于二维射线追踪的迭代反演方法,并成功地应用于野外资料。这个解释方法是由迭代射线追踪和更新速度模型来实现的。在每次迭代中,通过一个线性方程组将资料和由射线追踪得到的当前模型的走时差值与未知的速度扰动量联系起来。用阻尼最小二乘法求得这个速度扰动量,并因此更改模型。这个迭代过程一直进行到合成走时和野外资料间的拟合达到给定的误差范围之内,或者直到走时拟合不再有改进为止。我们把这个方法应用于结晶岩中井间试验所得的野外资料(Wong,1983)。地震图的频率范围是1—6.6kHz,允许速度结构的分辨率达到几米的数量级,所得的速度图象与其它地质和地球物理数据符合得很好。     

一种改进的线性走时插值射线追踪算法

线性走时插值法（LTI）在走时的计算中,由于射线方向考虑不全,计算得到的节点走时不一定最小,导致追踪的射线路径无法满足最小走时.针对这一问题,本文提出了一种改进的射线追踪算法,通过采用多方向的循环计算,得到所有计算节点的最小走时,使追踪到的射线路径能真正满足最小走时,以确保射线追踪的精度.模拟实验结果表明,在介质速度变化剧烈的结构中,该算法与传统的LTI算法相比,有效地提高了射线追踪的精度.     

三角网格剖分下速度与反射界面的同时反演

为了解决复杂速度模型中的走时正、反演问题,例如:含不规则起伏地表、不规则地下波阻抗界面、以及不规则速度异常体的复杂地学模型,本文采用三角网格单元模型参数化下的分区多步改进型最短路径算法,实现了多震相地震射线的追踪计算,结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合同时反演成像的方法技术.当界面起伏较大时出现散射,从而造成散射点所在区域射线密度过密,导致该区域内速度和界面的过度更新.为了克服上述问题,我们在同时反演中引入了射线密度的概念,从而有效地解决了上述过度更新问题.数值模拟实验表明:采用三角网格单元进行模型参数化,可保证在复杂模型中的正演计算具有较高的计算精度;同时反演中可以准确地刻画不规则异常体和不规则反射界面.因此,本文提出的走时成像方法技术具有较广的实用价值.     

半监督最小二乘支持向量机的研究及其在海上油田储层预测中的应用（英文）

勘探开发初期海上油田钻井少、井间距离大,在应用地震多属性分析技术预测储层参数的过程中,直接采用监督最小二乘支持向量机算法预测精度较低。本文将最小二乘支持向量机与半监督学习理论结合,提出基于最小二乘支持向量机协同训练的半监督回归模型,并在模型训练过程中引入矩阵迭代求逆的方法,提高模型训练速度。利用UCI数据集实验研究,对比了半监督与监督最小二乘支持向量机模型,结果表明,半监督学习机制能够有效地提高最小二乘支持向量机的泛化性能,且随着训练样本的减小,效果更加明显;同时对比了半监督最小二乘支持向量机与半监督k-临近算法,结果显示,在小样本建模中,半监督最小二乘支持向量机有着更高的预测精度。最终将半监督最小二乘支持向量机运用于锦州工区,预测该区的砂体及储层孔隙度的分布,获得了较好的地质效果。     

利用最小二乘配置在球面上整体解算GPS应变场的方法及应用

通过对位移与应变微分公式的分析,给出了球坐标系下应用最小二乘配置法由GPS数据整体解算应变场的方法.首先,对模拟数据的理论应变场与采用全部采样数据和随机抽取50%限定采样数据的应变场结果的对比分析,初步验证了最小二乘配置球面整体解应变场方法的有效性.然后,利用1999~2001和2001~2004两期中国大陆区域网GPS速度场数据给出了基于最小二乘配置球面解获得的应变率场结果,概要分析了此期间中国大陆应变率场的动态变化.最后,讨论了最小二乘配置在球面上计算应变场需要注意的问题和建议.