Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

泰青威天然气管道临朐段水毁灾害分布特征与风险评价

对泰青威天然气管道临朐段地质灾害发育特征进行调查研究发现,临朐段地质灾害类型主要为水毁灾害,具体可分为坡面水毁、河沟道水毁和台田地水毁。野外地质灾害实地调查临朐段管道沿线共发现地质灾害15处,其中坡面水毁点6处,河沟道水毁点4处,台田地水毁点5处。采用定性与半定量相结合的评价方法对其进行地质灾害风险评价,结果表明: 地质灾害风险等级较高的有1处,占6.66%; 风险等级中等的有4处,占26.67%; 风险等级较低的有10处,占66.67%。根据管道沿线地形地貌、地质灾害发育密度、风险等级等因素,划分地质灾害中易发区61 km,低易发区10 km,管道沿线以地质灾害中易发区为主。最后,针对不同类型、不同风险等级地质灾害提出了相应的防治消减措施,为后期管道安全运营和风险整治决策提供了有效的技术依据。     

用接收函数方法研究中国及邻区上地幔间断面的埋藏深度

使用中国及邻区的18个数字化地震台的宽频带三分向远震记录图,筛选计算得到了263条径向接收函数,并拾取了在41Okm和660km间断面上的透射转换震相P4S、P660S与P波的到时差.通过调整各台站之下的速度模型中这两个间断面的深度,使理论接收函数的P4S、P660S与P波的到时差与观测值相符,最终确定这两个间断面的实际埋藏深度.结果表明,各台站之下4km间断面的平均深度为403km,具有明显的横向不均匀性和区域性,其中北京台、西安台、兰州台和恩施台等的埋藏深度较浅,琼中台、CHTO等台较深.660km间断面的平均埋藏深度为663km,也存在明显的横向不均匀性,其中牡丹江台和MAJO台的较深,琼中台和HYB台的较浅.     

区域土壤环境地球化学研究——异常成因判别·环境质量·污染程度评价的思路与方法

从概念内涵和应用目标出发,认为土壤环境地球化学质量评价应该分为质量现状和污染程度两个方面.在此基础上探讨了区域地球化学资料在土壤环境现状质量评价、异常成因识别、污染元素累积速率和土壤污染程度评价研究的基本思路和方法.     

塔中地区良里塔格组海平面变化对高频层序和沉积演化的控制

【研究目的】碳酸盐岩地层易受海平面变化影响而发育高频层序,但是反映海平面变化的地化指标的精度不足,导致高频层序对海平面振荡的响应认识不够深入,高频海平面变化对台地边缘礁滩沉积的高频层序和沉积演化的控制尚不明确。【研究方法】论文研究选取塔里木盆地塔中地区某钻井的上奥陶统良里塔格组连续取心资料,通过密集取样开展微相分析和碳、氧同位素测试分析。【研究结果】该井良里塔格组中上部发育7种微相类型,沉积演化分析表明良里塔格组沉积中—晚期,沉积环境由潮坪向生物礁-颗粒滩再到较深水的开阔台地演化。碳、氧同位素测试显示δ¹³C为0.5993‰~1.6228‰（均值1.139‰）,δ¹⁸O值为-8.3608‰~-5.1452‰（均值为-6.790‰）;δ¹³C和δ¹⁸O的振荡变化与微相变化和沉积旋回对应良好。古海洋条件分析表明,良里塔格组沉积时期气候温热,礁、滩体发育层段样品的古温度值最高,代表着高的碳酸盐产率;Z值所反映的古盐度在底部潮坪沉积段最高,在高频旋回的顶部古盐度均明显减小,指示可能存在大气淡水的影响。与塔北地区南缘和巴楚地区的δ¹³C和δ¹⁸O存在一定差异,主要受控于沉积环境和古水深。【结论】良里塔格组中—上部存在至少3个完整的海平面升降变化周期,内部又包括至少2个次级周期;不同级别的海平面振荡主要受古气候变化驱动,并控制了高频层序的发育。海平面先缓慢上升再持续稳定或小规模振荡,是连续厚层的生物碎屑和砂屑颗粒滩沉积发育的有利条件;良里塔格组沉积晚期海平面快速上升,沉积环境向滩间洼地和较深水的开阔台地转变。在高频层序顶部海平面的小规模下降导致的早期岩溶作用是礁滩体储层质量改善的重要因素。     

平原地区机载激光雷达数据的抽稀算法分析

目前,机载激光雷达点云数据在测绘行业中的应用还存在较多的瓶颈。为了使机载激光雷达点云数据更好地服务等值线等数据的生产,发挥其高效和高精度的优势,本文归纳、总结了国内外现有的LiDAR点云数据抽稀算法,并通过对比分析现有LiDAR点云数据抽稀算法存在的优缺点,如系统抽稀、格网抽稀、TIN抽稀和坡度抽稀等算法,结合平原地区激光点云在实际生产中的应用,研究了更适合平原地区点云数据的抽稀方法,通过大量的数据测试和试生产。结果表明,该方法可以在应用项目精度约束下保证数据质量,减少了后期数据处理应用的难度,提升了后续成果数据的质量,提高了作业生产效率,对机载激光雷达点云数据在测绘行业中的应用推广具有重要的现实意义。     

基于稳定性耦合分析法的余王扁边坡稳定性分析

边坡稳定性分析的极限平衡法与有限元法(FEM)的耦合分析法,首先利用有限元法分析获得边坡岩土体的整体应力场,在此基础上利用极限平衡法进行边坡的稳定系数求解。该方法既反映了边坡的稳定和变形之间的关系,又克服了极限平衡法与有限元法的不足,使二者的优点相互补充,获得的稳定系数基于极限平衡理论体系,可以同传统的稳定系数评价体系接轨。以西安市雁塔区余王扁削坡后边坡为例,用稳定性耦合分析法对其进行了稳定性分析,并把分析结果与各种极限平衡法的计算结果进行了比较,证明了耦合分析方法的可靠性和可行性。     

基于Janbu法的边坡整体稳定性滑动面搜索新方法

对于边坡整体稳定分析,找到了一种用随机角来搜索随机滑动面的新方法。在此过程中,采用简化Janbu法计算安全系数。这种新方法产生随机滑动面的过程是：先假定滑动面与边坡上、下缘的交点,在某一初始方向的情况下,以下缘交点为基点开始每随机一个角度方向,使其与对应的滑动体土条划分竖线相交,直到随机角度方向与上、下缘交点的连线重合,将这些交点依次连接,即可得到一条随机曲线（或随机滑动面）。通过随机模拟曲线,发现随机角物理意义明确,具有以下两个独特优点：①只有一个随机条件,随机因素易于控制,便于计算机编程;②为保证其随机曲线为一光滑曲线,可对随机角进行适当优化,并且在条分数足够大时,能够模拟任何曲线。针对以往随机搜索精度不高的问题,采用近似曲线来代替随机曲线的方法加以改进,并且将近似曲线的生成进行了简化。经各种实例对比分析后,得出此方法计算得出的安全系数与以往研究成果颇为接近,得到的最危险滑动面亦与实际情况相符,因而可说明此方法能够适用于边坡整体稳定性分析。     

航空物探勘查金属矿产时提取弱信息的解释方法

摘 要　 结构逻辑法是一种提取航空物探局部弱信息的新的解释方法。在介绍该方法原理的 基础上‚探讨了利用这种方法提取航空物探局部弱信息的效果及其在矿产资源快速勘查评价 中的作用。     

既有路基下浅埋隧道开挖引起地层的位移及应力解析解

在既有公路下进行浅埋隧道开挖会引起地层产生位移,导致地表沉降。为了保证施工安全,快速精确地计算地层位移及应力,有效地控制地表沉降,对浅埋隧道的设计与施工均十分重要。在镜像法原理的基础上,引入等效模量当层法,推导出了在有上覆不同材料下浅埋圆形隧道开挖引起地层位移、应力及地表沉降的解析分布规律。计算结果表明,既有路基下开挖隧道引起地层产生的位移和应力解析解能够很好地体现上层高弹性模量路基对沉降及应力的扩散作用,可以求得隧道开挖后洞周地层包括上覆路基及下卧层中任意位置处的位移、应力以及地表沉降,为在双层地基材料中开挖浅埋隧道引起地层位移及地表沉降的计算提供了一种新的解析方法。