远距离井间电磁成像测井快速反演算法

远距离井间电磁成像方法可以实现500 m以上地层剖面成像,在资料反演中反演地层的精度以及速度是研究的重点和难点.本文提出利用格林函数方法取代差分近似方法计算雅克比矩阵的快速算法,通过与全空间差分近似法对比,验证了快速雅克比矩阵算法的正确性,格林函数方法极大地提高了反演成像效率;通过建立地质模型进行仿真,验证了格林函数方法求取雅克比矩阵的高斯牛顿反演算法完全可应用于远距离井间电磁成像,实现了对二维地层电参数以及边界位置的确定.反演速度是差分近似法的数十倍.     

小波与曲波变换探地雷达数据去噪对比分析

探地雷达勘探工程目标及观测工程环境越加复杂给其精确的数据处理带来极大挑战,高效的探地雷达数据去噪算法是当前关注的重要研究领域.基于阈值去噪思想,小波变换和曲波变换去噪算法在探地雷达数据去噪应用中受到限制,有必要开展上述两种去噪算法适用性和实用性系统评价及改进.基于传统高阶相关统计阈值和块状复数域阈值函数,本文开展了小波变换及曲波变换去噪算法在合成含噪数据去噪效果对比分析;提出窗口高阶相关统计阈值小波变换去噪算法、探讨了块状复数域阈值函数取值变化对曲波变换去噪效果影响规律.通过对实测数据去噪分析,验证了窗口高阶相关统计阈值小波变换和估计块状复数域阈值函数曲波变换去噪算法的可行性及有效性.     

层状饱和场地中圆柱形基础的阻抗函数

本文利用间接边界元法,建立层状饱和场地中三维土-基础动力相互作用模型,研究了孔隙中流体的存在、场地动力特性和圆柱形基础的高度对基础阻抗函数的影响。数值结果表明,孔隙中流体的存在对基础阻抗函数有明显的影响,尤其是水平和竖向分量;场地动力特性也对基础阻抗函数有明显的影响,层状饱和场地中阻抗函数围绕相应均匀半空间结果振荡,且基岩与土层刚度比越大或土层厚度越小,振荡的幅度越大,土层厚度越小,振荡的周期越小;另外,基础高度也对阻抗函数有明显的影响,随基础高度增大,阻抗函数逐渐增大。     

新疆天山地区壳幔S波速度结构特征及变形分析

天山地区地质构造复杂,地震活动频繁,其壳幔变形和深部结构一直受到学者们的高度关注.然而,由于天山地区地震台站资料较少,致使壳幔变形研究结果与解释存在诸多争议.本研究利用在天山地区（40°N-46°N,78°E-92°E）新布设的11个流动宽频带地震台站和该地区39个固定台站的观测资料,采用接收函数与面波联合反演方法,获得了研究区地壳厚度及壳幔S波速度结构.反演结果显示天山地区（41.5°N-44°N,78°E-88°E）平均地壳厚度为56 km,塔里木盆地（40°N-41.5°N,79°E-90°E）、准噶尔盆地（44°N-46°N,82°E-90°E）和吐鲁番盆地（42°N-43°N,88°E-90°E）具有较厚的沉积层,地壳平均厚度为43 km、53 km和46 km,整体表现为天山厚、盆地相对较薄的特征;在研究区南天山的最高峰（42°N,80.5°E）及北天山的最高峰（43.5°N,86°E）附近,中下地壳存在较厚的低速层,我们认为在强烈挤压作用下低速、低强度的中下地壳强烈变形可能是导致该区域快速隆升的主要原因.在研究区中部,位于塔里木盆地与准噶尔盆地之间的天山地区,中下地壳及上地幔均存在低速层,且盆地莫霍面向天山倾斜明显.结合前人的研究成果推测,在南北向构造挤压应力作用下,塔里木盆地与准噶尔盆地发生了向天山造山带方向的双向壳幔层间插入俯冲.在研究区东部,塔里木盆地东北缘与天山东部接触带的地壳内没有明显的低速层,推测应处在早期挤压变形状态,该区域的壳幔边界为缓变的速度梯度带,可能与上地幔热物质侵入或渗透有关.     

空间波数混合域磁异常场积分解三维数值模拟

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.     

肺亚实性结节CT图像定量变化对浸润的判断

目的:定量评估肺亚实性结节（SSN）内非实性成分从基线到随访CT扫描之间的变化,以鉴别诊断有浸润表现的SSN,即微浸润性腺癌（MIA）和浸润性腺癌（IAC）。材料与方法:回顾性纳入有基线和复查CT,并有组织学结果的SSN。测量指标为SSN整体最大直径,非实性成分CT值均值、最大CT值、最小CT值和CT值标准差。以时间为因变量,用线性回归评估SSN测量值从基线到复查扫描的进展。采用受检者工作特征曲线（ROC）和曲线下面积（AUC）评价测量值年化改变对有浸润表现SSN的诊断能力。结果:71名患者（男32,女39）中纳入78个SSN,其中良性19个,原位腺癌10个,微浸润腺癌36个,浸润腺癌13个。所有SSN分类的直径均与时间呈显著正相关（P<0.05）。有浸润表现SSN的非实性成分CT值标准差与时间显著正相关（P=0.037）。良性和原位腺癌中,所有4个CT值指标与时间无相关性（P>0.05）。ROC分析提示CT值标准差年化改变可以明显诊断有浸润表现的SSN（P=0.014）,诊断敏感性、特异性和AUC分别为68.7％,77.8％和0.747。结论:在SSN的复查过程中,良性、原位腺癌、MIA和IAC的直径都有增加,评估CT值标准差的变化能够鉴别诊断有浸润表现的SSN（MIA和IAC）。      

基于三维应变格林函数反演中小震震源机制

随着复杂速度结构反演的发展和高性能计算能力的提升,基于高精度3D介质模型计算格林函数反演震源机制更具可行性。中小地震因具有更好的覆盖和近似点源效应,在区域结构成像中有着广泛的应用。基于波形类的反演方法如波动方程层析成像\,全波形反演都需要震源机制解,而传统的震源机制反演方法不能很好地应用于中小地震。采用有限差分法构建应变格林张量(Strain Green Tensor,SGT)数据库,将合成波形和实际波形按震相截窗并滤波到不同的频带范围,先通过最小化互相关走时差来进行震源重定位,再通过最小化波形残差反演震源机制。通过合成数据测试验证方法的正确性,随后将该方法应用于青藏高原东部边缘龙门山地区,反演一系列M_W3.4～5.7的中小地震震源机制。由于应变格林张量数据库可预先构建,该方法可以应用于(近)实时震源机制解反演。     

元谋地电场观测数据小波阈值去噪

为了解决元谋地电场观测数据的噪声干扰问题,采用小波阈值去噪的方法,根据元谋地电场观测数据特征,选取适当的阈值及不同的小波函数、分解层数进行研究。综合各小波函数去噪效果,通过对相对误差、均方根误差、信噪比等指标进行评价,确定适合元谋地电场观测数据去噪的阈值为软阈值函数、分层阈值;最优小波函数分别为Sym4、Sym5、db6和db8;分解层数分别为5、7、9。最后应用研究结果处理实际地电场观测数据,取得了较好的去噪效果。     

上海区域量规函数研究

量规函数是用来校准补偿地震波随距离的衰减的,本文对上海区域台网的量规函数进行重新计算,同时和全国通用的量规函数加以对比,观察是否需要对上海区域量规函数进行修正。并使用统计平均震级和M_L震级定义两种方法来计算量规函数。      

碰撞非相干散射理论谱宏观法和微观法的比较研究

非相干散射理论谱是理解非相干散射雷达探测的理论基础,其求解过程分为微观法和宏观法.在考虑碰撞不考虑磁场的条件下,本文使用这两种方法对理论谱分别进行求解,给出了关键求解过程,并对等离子体的归一化导纳和极化率以及非相干散射的微分散射截面和谱密度函数进行分析,得到了它们各自之间的数学关系,概括性分析了两种方法在理论选择、最终结论及多种因素下求解的异同点.