2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

基于走时的保幅偏移方法

振幅随偏移距变化是描述储层特征的重要方法之一,保幅偏移方法就是使偏移剖面能够反映出振幅随偏移距的变化.本论文中的保幅偏移是以走时为基础,主要的方法是采用走时的双曲线展开法,通过走时的二阶空间导数来确定波前曲率.该方法通过建立在大网格上的走时表来确定插值系数,将大网格插值成为较为精细的网格,这样就节省了数据的存储空间.对于相同的网格密度,通过插值来计算走时表比采用程函方程有限差分法直接计算走时要节省5至6倍的时间.走时的插值系数还可以用来计算几何扩散因子、权函数,不仅提高了成像质量,还大大节省了计算时间.     

非椭圆各向异性介质高阶优化走时计算方法及应用

随着地震勘探的深入,对勘探精度的要求越来越高,目前,提高勘探精度的主要方法之一就是提高地震波走时计算的精度,地震走时精度直接影响如动校正、静校正和偏移成像质量.根据地质模型假设条件的不同,有不同的走时计算方法,常用的有各向同性介质假设和各向异性介质假设条件下的走时计算方法,各向异性假设具有比各向同性假设更好的适用性,具有更高的走时计算精度,能够满足大偏移距和高陡构造情况下成像的要求.本文根据Alkhalifah的VTI介质中走时计算方法,以Fomel的TI介质中群速度的计算为基础,利用各向异性介质中Snell定律,提出了在非椭圆各向异性介质中的高阶优化走时计算方法,通过理论和实际数据的应用,表明该方法克服了常规各向异性走时方法中不能够计算大偏深比情况下地震波走时的缺点,具有走时计算精度高、便于实际应用等特点,对叠加和偏移成像有广阔的应用.     

层状各向异性介质转换波克希霍夫叠前时间偏移

在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.     

远震相对走时数据快速计算方法及应用

远震层析成像技术已成为研究地球深部速度结构最为主要的方法,而其所用数据通常为初至震相的相对走时残差值.本文针对信噪比差的远震波形资料提出一种计算远震相对走时残差数据的快速方法,整个计算过程全部由计算机自动完成.该方法主要有两大特点:(1)利用互相关方法计算波形间的偏移时间时引入波形相位信息,有效压制噪声,增强有用震相的识别;(2)改进常规方法中计算相对走时残差的理论公式,避免求解超定方程组,提高计算精度.为验证新方法的可行性和实用性,以长江中下游地区远震波形资料为例,采用常规方法和新方法分别计算相对走时残差数据,对比结果显示:新方法比常规方法的工作效率和数据精度都高.     

基于高斯束传播算子的成像域走时层析成像方法

层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性.     

多步法快速射线追踪与圆台型高精度走时外插计算

为更好地适应复杂构造的地震偏移成像,本文提出了一套快速射线追踪算法和一种高精度的走时外插计算方法.采用线性多步法的预测-校正公式求解射线追踪方程组,与传统的四阶Runge-Kutta法相比,提高了计算效率.在网格节点上的走时计算中,应用一种基于圆台的外插方法,该方法以射线的方向为轴确定圆台,将轴上的走时外插到圆台内的网格节点上.与传统的矩形体外插方法相比,圆台走时外插方法提高了计算精度,且具有更好的稳定性.另外,该方法利用稀疏分布的射线即可获得高精度的走时表,节省计算量,对复杂构造的偏移成像非常有利,尤其是三维偏移.最后通过逆散射偏移成像算例,验证了算法的有效性和适用性.     

角度域叠前时间偏移:非均匀介质中的单程波方法

从单程波方程出发,推导出角度域叠前时间偏移的走时、入射波与反射波夹角、成像幅值计算方法;构建了可直接生成角度域成像道集的叠前时间偏移方法与偏移流程.文中定量分析了速度梯度对走时、角度、幅值的影响,给出了可更好考虑介质非均匀性的角度、幅值计算方法.理论模型以及实际数据验证了本文方法的有效性.     

基于快速推进与波前构建的联合3D走时计算（英文）

3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点,但快速推进法在震源附近区域的计算精度不高,降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时,在剩余区域使用快速推进法计算走时,由于模型中绝大多数网格节点走时是通过快速推进法计算的,故新方法保留了快速推进法高效的特点,同时由于震源附近网格节点走时精度的提高,整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。     

地震波走时场模拟的快速推进法和快速扫描法比较研究

地震波走时信息在叠前偏移、叠前速度分析、地震层析成像、走时反演及地震定位等中都有重要应用.快速推进法因其理论完善、精确灵活,无条件稳定,近年来已在走时计算领域得到广泛应用.快速扫描法作为求解一阶非线性双曲型偏微分方程的高效方法,已在图像处理、计算机视图、控制论等领域得到有效应用,且在走时计算方面有所应用且展现了广泛的应用前景.本文介绍了两种方法的基本原理且(通过均匀介质模型、局部低速体模型和Marmousi模型)把两种方法做了详细对比.研究结果表明:1)基于逆风差分格式的快速推进法和快速扫描法对纵横向速度变化很大的不均匀介质依然有很好的稳定性和适用性,均可以准确地计算地震波初至走时;2)对于相同的模型和在相同的计算条件下,两种方法的精度相当,但快速扫描法所耗的CPU时间较快速推进法明显减少,效率显著提高.     

江西省地震局虚拟专用网(VPN)建设

叙述了在FreeBSD环境下构建基于PPTP的VPN网关的过程和技术，分析了基本概念、实现方式以及在江西省地震信息网的实际应用等。     

山西大同数字遥测地震台网记录的震相特征

从地球及地壳的分层结构出发，使用较成熟的理论，从地震波传播的基本原理入手，结合多年的震相观察分析工作实践，选取了山西大同数字遥测地震台网近年来的一些典型震相进行了分析，得出本台网记录到震中距小于100km的地方震震相有Pg，P11，Sg，S11震中距在100km～500km的近震震相有Pn，Pg，P11，Sn，Sg,S11对比分析了在本台网内及周边地区发生的爆破与天然地震在震相上的差异性，并对其特征作了较详细的描述。     

兰州地震研究所发展的三段历史回顾

对兰州地震研究所的三段历史作回顾,分别为国家地震局、中国地震局兰州地震研究所的体制沿革变迁,中国科学院地球物理研究所;中国科学院兰州地球物理研究所观象台早期创建;兰州地震研究所硕士点早期阶段。     

山东数字地震台网子台背景噪声位移值的分布特征

利用山东数字地震台网的实时记录波形资料,采用在频率域积分方法,计算出各台站背景噪声的位移值。然后进行统计分析,得出了山东地区背景噪声位移值的空间分布和时间变化特征：台站背景噪声幅值与不同地区地质环境、岩性密切相关,位于第四系黄土覆盖层较厚的鲁西聊城、菏泽地区,鲁西北德州地区和鲁北东营、滨州地区各子台的背景噪声位移值较大,鲁中、鲁南山区和胶东半岛地区的基岩台站背景噪声位移值较小;同时,白天噪声值大于夜间噪声值,但深井观测台站的背景噪声值相对稳定。此项研究为绘制台网监测能力图提供基础数据,并为计算各台站的场地响应及震级修正值等测震学研究提供基础资料,促进台网的改建和发展。     

裂隙岩体渗透性空间分布的指示克里格估值

本文较详细介绍了指示克立格估值的基本方法和原理，以指示变异函数为基本工具分析了向家坝坝址区裂隙岩体渗透性空间分布的结构特征，表明该址区裂隙岩体渗透性存在明显的各向异性特征；在此基础上用指示克立格法对未采样点处进行估值，估值结果显示坝址区裂隙岩体渗透性存在明显的壳状渗透结构特征.通过这一实际应用，表明指示克立格法可以较好地描述裂隙岩体渗透性的空间分布规律.     

对齐误差对测试地震计自噪声的影响

本文在理论上研究了对齐误差对测试地震计自噪声的影响。指出对齐误差影响到2台仪器记录间的计算互功率谱密度, 进而影响到被测仪器记录间相干函数计算, 最终影响到被测仪器的自噪声计算。发现当存在对齐误差时, 在信噪比较高频段计算得到的自噪声水平比仪器的实际自噪声水平要高, 在一般的实际测试中可能高20 dB以上。为了消除对齐误差对自噪声测试的影响, 提出通过投影逆变换把观测记录变换到标准坐标系下, 再用变换后的观测记录来计算地震计自噪声。给出了基于以上思路和传感参数对比测试方法的地震计自噪声计算处理流程, 合成数据、 实际测试数据的处理检验计算结果表明, 本文方法可以有效克服地震计自噪声测试中对齐误差的影响。     

萧山-球川断裂富阳—建德段第四纪活动性研究

萧山-球川断裂是浙江地区1条大规模的北东向断裂,在该断裂附近曾发生过多次破坏性地震。本文通过地质调查、电法勘探和探槽开挖等方法,对萧山-球川断裂富阳—建德段第四纪活动性进行了研究。萧山-球川断裂富阳—建德段的遥感影像线性特征清楚,对地貌的控制作用较为明显,对山前的第四系发育有明显影响。通过对断裂露头剖面的分析,认为萧山-球川断裂富阳—建德段活动性质以走滑兼具逆冲为主。通过本次开挖的富阳峙山村探槽,结合ESR年龄测定,判定该断裂在第四纪早、中期有过活动,但未断错上覆中更新统上部地层,其最新活动时代为早、中更新世。     

层状地基土双自由度集总参数模型的参数识别

在土-结构的动力相互作用简化分析中,常采用双自由度集中参数模型简化土层.这种模型适合各向同性的单一土层,而实际工程场地的地基土多是各向异性的成层土层.为了将这一模型应用于此类土层,本文根据实际地震动记录,提出了一种变阻尼最小二乘反分析方法,可用来识别模型中的各个参数值.文中提出的方法具有适用范围广、计算精度较高的优点,并已通过编制FORTRAN程序加以实施和应用.     

山西夏县中心地震台温泉水化学特征

采集夏县中心地震台温泉水及周边水点样品进行水化学组分特征分析。利用矩形图、Na-K-Mg三角图等方法,初步分析温泉水的水质类型、水-岩平衡状态、热储温度以及循环深度等,并结合氢氧同位素组成特征,初步分析温泉水补给来源。研究表明,夏县中心地震台温泉水化类型属于Na-Cl·SO₄型,水-岩反应属于部分成熟水,热储温度为148.8℃,循环深度为3.91 km,补给源主要为大气降水,温泉为断裂型温泉。以上结果可为该台流体异常分析提供基础研究资料,为夏县地区水化研究奠定基础。     

倾斜仪时号·光源灯控器的改进

为解决地震台光记录倾斜仪灯光系统故障率高的问题 ,研制了一种时号·光源灯控器。采用了数字钟报时及定闹信号对直接延时时号灯进行控制的措施 ,使仪器的无故障运行期及蓄电池的使用寿命得以延长 ,明显提高了观测质量。