一种改进的地震反射层析成像方法

针对复杂介质的地震反射走时层析成像存在数据拾取困难问题,本文提出了一种新的地震反射层析成像速度模型建立方法,该方法用速度和地震射线走时描述模型,用地震反射波走时、地震波在源点和接收点处的传播方向信息反演模型.为提高反演的稳定性和计算效率,引入了Hamilton函数描述射线,在相空间计算反演所需的射线路径和目标函数对模型参数的导数,对理论模型和实际地震资料进行了试算,试算表明该方法对复杂介质具有较强的适应能力.     

基于改进Runge定理的调和延拓及其相应的边值问题

利用改进的Runge定理引入了外大地水准面和外重力位等新概念 ,建立了关于外扰动位的边值问题 ,该问题同时具备了Stokes问题和Mododensky问题的优点 .作为一个整体 ,讨论了外正高的计算和地面重力的归算 ,给出了在O(T2 )量级的精度下确定外大地水准面、地面以及地球外部重力场的方法 .     

改进的角变换零偏移距成像方法研究

本文基于椭圆展开成像原理, 引入地震波入射角和地层倾角两个参数推导出了角变换零偏移距成像方法, 该方法只使用两个角度实现零偏移距成像, 但成像存在强干涉噪声问题.针对该问题, 并从实际生产应用方面考虑, 推导出了分别使用地震波入射角、地层倾角、速度和地层倾角进行成像的三种改进方法, 通过对比分析, 优选实现了含有速度和地层倾角的改进的角变换零偏移距成像方法, 通过数学分析、数值模型和实际资料处理证实了新方法的有效性: 新方法没有地下反射层为水平的前提假设条件, 地下反射层可以倾斜或弯曲; 均匀介质条件下新方法可以求取真实的地下速度和地层倾角, 实现真正的零偏移距成像处理, 这和CMP类方法有本质区别; 对复杂地下介质, 新方法也可以实现较好的零偏移距成像处理.

     

新丰江6.1级水库地震前后的波速比变化

广东省新丰江水库区于1962年3月19日发生6.1级地震。本文分析了1961年以来十四年间,特别是主震前后纵、横波速比Vp/Vs的变化过程。结果表明:波速比Vp/Vs在主震前出现了明显的负异常,总异常时间约为十一个月,于震前一个半月返回正常值,在主震震源附近的深水峡谷区异常幅度最大,异常时间和空间范围比一般构造地震的要小.后期几次强余震前也有异常变化。虚波速度在主震前存在某些异常特征。上述波速比的变化可以用扩容模式较好地解释。 波速比值是由两个台站的记录得到的。对虚波速度的计算采用了联合地震参数的方法。     

基于粒子群优化的多目标体中梯电阻率异常反演

在采用中梯装置的电阻率剖面法应用中,特别是在环境、水文和工程等领域,经常遇到需要对多个异常目标体进行快速定位以便及时进行相应处理的情况.利用倾斜椭球体来近似模拟这些电阻率异常目标体,使得正演计算可以采用解析表达式来实现,提高正演计算的时效性.相应的地球物理模型也可简化为由椭球体个数、中心点位置、倾角、轴径以及电阻率等参数构成的粒子,多个粒子组成的粒子群在粒子群优化算法的控制下在给定的搜索空间中并行地搜索最优模型.通过对粒子群优化算法参数的合理设计,利用其良好的全局与局部均衡的搜索能力实现对多个异常目标体的同时反演.数值实验结果表明该反演方法能有效实现对多个目标体的同时反演,计算速度快、反演拟合精度较高、同时具有一定的抗噪音能力.快速的多目标体反演,可以实时准确的定量解译中梯剖面法圈定的异常目标体,较好地满足工程等领域的高时效性要求.     

新丰江6.1级水库地震前后的波速比变化

广东省新丰江水库区于1962年3月19日发生6.1级地震。本文分析了1961年以来十四年间,特别是主震前后纵、横波速比V_p/V_s的变化过程。结果表明:波速比V_p/V_s在主震前出现了明显的负异常,总异常时间约为十一个月,于震前一个半月返回正常值,在主震震源附近的深水峡谷区异常幅度最大,异常时间和空间范围比一般构造地震的要小.后期几次强余震前也有异常变化。虚波速度在主震前存在某些异常特征。上述波速比的变化可以用扩容模式较好地解释。 波速比值是由两个台站的记录得到的。对虚波速度的计算采用了联合地震参数的方法。     

三维分块倾斜界面的反演及其应用

介绍了利用反射波走时反演介质层速度和包含直立断层的三维界面的计算方法．对各层界面利用多个倾斜或平界面方程来描述．给出了数值计算实例．结果表明，反演解与真模型很接近，说明计算方法是有效的．利用该方法处理了唐山地震区的ＰｍＰ反射波资料，获得了该区域莫霍界面的三维分块形态．唐山附近莫霍界面错断与该区域１９７６年唐山地震及一系列余震的发生有密切的联系．     

基于重力异常反演的三维地震波速度结构分析方法研究

利用重力异常反演测试三维地震波速度结构,存在解不唯一、可靠性不高的问题。将面波反演充分融合到重力异常反演方程中,降低传统反演方法的非唯一性,并提升可靠性。以川滇地区为例,采用融合后的重力异常反演方法分析三维地震波速度结构。通过速度和密度的关系转换,得到对应的重力异常数据。由于面波频射数据主要对地震波横波速度敏感,因此将重力异常数据和初始横波速度相连,依据地震波速度和岩石密度之间的关系,获取重力异常反演方程,用于分析速度结构。选取21.6°~34.2°N、97.1°~105.9°E范围内的川滇地区活动块体作为实验数据,经过实验分析发现：使用该方法迭代反演川滇地区地壳上地幔顶部横波速度,重力异常数据和面波频射数据的残差值分别是6.24 mGal和0.027 km/s,实际拟合效果较好;分析该地区不同深度切面横波速度发现,在24 km深度处,上地壳中含有相对低速层,在44 km深度处,中下地壳中存在低速层;且该方法分析川滇地区三维地震波速度结构解的分辨率较高。     

基于概率成像技术的低纬度磁异常化极方法

化极转换是磁异常解释的重要基础,为了克服在地磁纬度较低的地区尤其是磁赤道附近化极不稳定的问题,出现了多种化极方法．本文基于概率成像技术提出了一种等效物性的反演方法,实现对地下等效场源的反演成像,取得了对低纬度磁异常稳定化极的效果．化极反演中逐次对剩余异常进行反演成像,实现由概率模型到物性模型的复杂映射,避免了类似反演中需要大型方程组求解等问题,并将概率模型的构制、物性参数的反演和反演评价有机地集成到一起,加速了反演成像的进程,使反演成像速度与目前概率成像的计算速度达到可相比拟的程度．对理论模型和实际资料的计算表明,该方法对低纬度磁异常化极处理是稳定有效的,而且可以较好地压制噪声干扰,能够在噪声干扰条件下进行反演化极．     

基于聚焦合成的显微三维成像系统

本文利用显微镜物镜焦深范围有限,显微三维样本成像时需要一序列图像才能聚焦清晰特点,研究开发利用图像处理技术扩展显微成像焦深范围,应用改进Laplacian聚焦算子实现序列图像的融合显示,同时,提出了一种基于区域小波变换的序列显微图像融合算法并加以实现。用高斯插值算法得到三维显微图像的高度图,设计了一种网格轮廓线算法对高度索引图进行了三维重建显示和测量,并利用立体视觉原理生成了三维物体的左右视图和立体视图。应用此方法可以完成显微目标的表面三维参数的测量分析。     

基于背景噪声的宁夏及邻区瑞利波层析成像*

使用了宁夏及邻区90个地震观测台站资料,反演得到了该区高分辨率的瑞利面波层析成像结果。选用2012年1月至2013年12月的垂直分量连续波形数据,通过对台站对间进行波形互相关和叠加运算,得到各台站对间的经验格林函数,使用CPS(Computer Programs in Seismology)软件测量得到了4005条瑞利面波群速度频散曲线。选取信噪比大于10的经验格林函数对频散曲线进行了筛选,提取了3182条台站对间的频散曲线,反演得到了研究区周期6~50s,分辨率达0.5°×0.5°的瑞利波群速度分布图像。银川盆地在6~26s为明显的低速异常,且盆地内部的低速异常呈现出不均匀的减弱趋势;关中盆地在6~22s为条状低速异常,呈现出东南部的速度略强于西北部的横向不均匀性;在30~50s呈现出在银川断陷盆地及其以南存在一自北东向西南走向的大区域低速异常区,且表现出山地与盆地之间的主要活动构造带存在明显的地貌界线,即六盘山附近的深大断裂是中国东西两大构造单元的分界线。反演结果与研究区地质构造特征、地层地貌具有较好的相关性,且和银川盆地与穿过盆地的人工地震剖面结果吻合,为该地区活动构造带的动力学和强震发震机理研究提供重要依据。     

3S技术支持下鄱阳湖区地质构造调查分析

采取遥感(RS)技术对鄱阳湖区地质构造TM影像特征进行分析,建立湖区不同断裂遥感解译标志.利用人机交互方式解译出湖区不同断裂,借助地理信息系统(GIS)和全球定位系统(GPS)将该区的地质构造信息和相关地质资料进行空间叠加对比和实地验证分析的基础上,对该区地质断裂及活动性断裂构造进行深入研究和遥感影像新构造类型分区.     

青藏高原P波速度层析成像与岩石圈结构

利用中国西部地震台网的数据,通过体波层析成像反演了青藏高原及邻域的三维P波速度结构.根据地壳和上地幔的速度变化和构造特征,重点讨论了下地壳流动、地幔上涌、岩石圈减薄以及与藏北新生代火山岩和藏南裂谷系的关系等问题.分析表明,青藏高原中、下地壳平均速度偏低,低速区主要分布在拉萨和羌塘块体内部,随着深度的增加逐渐扩大到松潘—甘孜块体.上述低速区之间多被高速带分隔,暗示地壳中、下部的韧性变形被限制在特定的区域,不太适于产生贯穿整个青藏高原的大规模横向流动.此外,地幔上涌也并非普遍发生于整个青藏高原,而是集中在羌塘、松潘—甘孜以及喜马拉雅东构造结附近,导致上述区域的岩石圈地幔较薄,并且伴生火山活动和岩浆作用.此外,由于印度大陆岩石圈在向北俯冲,板片下沉过程中引起地幔上涌,热流物质有可能上升进入地壳,这一作用对藏北新生代火山岩和藏南裂谷系的形成以及中、下地壳的韧性变形产生了明显的影响.

     

小江断裂带周边地区三维P波速度结构及其构造意义

作为青藏高原的东南边界,小江断裂带在高原物质的侧向逃逸中发挥着重要的作用.本文利用流动地震台阵及固定台站的走时观测资料,对小江断裂带及周边区域的壳幔三维P波速度结构进行了研究.结果表明,在中上地壳,小江断裂带内部主要为低速异常,其东侧主要为高速异常.在中下地壳,小江断裂带中部为低速异常,北部和南部主要为高速异常,其中北部的高速异常可延伸到地表附近,南部的高速异常可一直延伸到上地幔.我们推测,小江断裂带中部的低速异常与深部热作用有关;北部的高速异常可能是晚古生代地幔柱活动导致大量基性和超基性幔源物质侵入地壳引起的,它的存在对青藏高原物质向南逃逸起到了一定的阻挡作用,可能是导致川滇活动块体北部次级块体快速抬升的重要因素;南部顶界面向北倾斜的高速异常体对川滇活动块体向南滑移起到了进一步的阻挡作用,导致其上覆的中上地壳低速异常区发生较强的变形和强烈的地震活动,同时在上地幔深度范围起到了稳定的作用,使其南部区域的介质受青藏高原物质向南挤出的影响明显减小.     

煤矿井巷间地质构造及其异常多波联合探测技术与应用

影响煤矿层安全高效开采的最主要的地质因素是断层构造及其异常,特别是落差在3~5m的中小断层构造对生产影响最大.为了保证开采工作的顺利进行,通常需在工作面回采前对面内地质条件进行探测与评价.利用地震波CT多波联合探测反演与解释,可精确确定工作面内的地质构造及异常范围.其弯曲射线追踪SIRT法反演可获得收敛性好的探测区域波速结果,以纵波波速结果为主,高速区代表断层构造迹线,低速区代表煤层增厚,同时结合其它参量成图,可有效识别各种地质构造及其特征.3241工作面的探测与回采实例有效地说明了,多波联合探测技术可为面内生产解决大量的实际地质问题.     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.     

Tunnel seismic tomography method for geological prediction and its application

Typical existing methods of tunnel geological prediction include negative apparent velocity, horizontal seismic profile, and the Tunnel Seismic Prediction （TSP） method as this technology is under development at home and abroad. Considering simpler observational methods and data processing, it is hard to accurately determine the seismic velocity of the wall rock in the front of the tunnel face. Therefore, applying these defective methods may result in inaccurate geological inferences which will not provide sufficient evidence for classifying the wall rock characteristics. This paper proposes the Tunnel Seismic Tomography （TST） method using a spatial observation arrangement and migration and travel time inversion image processing to solve the problem of analyzing the velocity structure of wall rock in the front of the tunnel face and realize accurate imaging of the geological framework of the tunnel wall rock. This method is very appropriate for geological prediction under complex geological conditions.     

TST隧道地震CT地质超前预报技术及应用

国内外的隧道地质超前预报技术正在发展之中,4目前应用的负视速度法,水平剖面法,TSP技术等,由于观测方式和资料处理方法比较简单,不能准确确定掌子面前方围岩的波速,从而致使地质界面定位不准确,围岩工程类别划分缺乏依据,亟待改进。本文提出的TST技术采用空间观测方式,资料处理方法采用偏移成像和走时反演成像,解决了掌子面前方围岩速度结构的分析问题,实现了隧道围岩地质结构的精确成像,适合复杂地质条件下的地质超前预报应用。     

云南通海盆地的浅层速度结构

为了研究通海盆地的地震动场地效应,在强地面运动的数值模拟计算中需要提供该地区的浅层速度结构模型.为此,在通海盆地开展了10 km高分辨率反射地震勘探和30个测点的面波联合勘探工作.综合本次探测结果与褐煤普查中的电测深和钻孔资料,得到了通海盆地的基岩埋深分布图.结果表明,在通海盆地中部,新生界厚度为100~300 m,整体表现出西南厚东北薄的趋势,盆地沉积层厚度变化明显受到小江断裂南段分支断裂的控制.浅层速度随深度分布表现出在深度相同的情况下,盆地边缘的速度高于盆地中部.     

俄罗斯贝加尔湖-日本仙台断面地震波速结构及其地质意义

本文以中俄、俄日学者合作所得到的地球物理资料为主,结合其它相关地质-地球物理数据,组构了俄罗斯贝加尔湖-日本仙台(BS)4000 km长断面,用于区域性大尺度地研究东北亚洲地壳结构和一系列地质构造问题.研究BS断面地震波速结果表明:(1)西伯利亚板块和黑龙江板块地壳结构变化较大,并可分为上、中、下部地壳,欧亚板块东部陆缘带地壳结构较简单,基本两分.贝加尔裂谷带下部地壳厚度比松辽盆地的薄约7 km,而上部地壳则相反,前者的比后者的厚约9 km.两个裂谷带在Moho界面之下的波速分布差异也较大.(2)结合前人认识,综合分析认为,贝加尔裂谷带属主动式裂谷,松辽盆地属于混合型裂谷.贝加尔裂谷形成动力主要来自地球构造圈B″层物质上涌所形成的地幔热柱的垂向作用,由BLV带佐证,松辽盆地形成动力主要来自太平洋板块斜向俯冲的中远程效应.(3)日本国所位于的西太平洋岛弧带是多地震带,除了太平洋板块俯冲产生的浅部效应、地壳中断裂与流体的直接作用等因素,本文指出仙台等速块的物性条件是岛弧带的主要不稳定因素.同时指出需要关注日本东海岸深约30~40 km的大级次地震的发生.