基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

Stability of source mechanisms inverted from P‐wave amplitude microseismic monitoring data acquired at the surface

We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots.     

基于在线高分辨率示踪技术的岩溶泉污染来源及含水介质特征分析——以重庆黔江区鱼泉坎为例

在线高分辨率示踪技术以其精度高、成本低的特点在岩溶地下水文系统研究方面得到较为广泛的应用。重庆市黔江区鱼泉坎岩溶泉的污染问题导致蓬东乡、五里乡及其周边3 000居民饮水困难,为了确定水源地主要的污染源,选取荧光素钠、天来宝2种示踪剂,应用在线高分辨率示踪技术确定了重庆某畜禽养殖有限公司的化粪池与鱼泉坎泉之间存在水力联系。两种示踪剂的回收率均较低,分别为:荧光素钠3.51%,天来宝2.56%,且两回收曲线存在一定差别,可能与天来宝的水溶性较荧光素钠低有关。基于荧光素钠初现时间和首个峰值出现时间计算的地下水最大流速与平均流速分别为27.09 m/h、22.32 m/h,反映鱼泉坎岩溶孔隙、裂隙较小,化粪池污染渗漏缓慢、运移时间长。示踪剂质量浓度历时曲线表现为多峰型,推断化粪池至鱼泉坎泉段至少存在4条较大的径流通道,表明示踪段表层岩溶泉裂隙发育的相对均匀性。欲恢复鱼泉坎泉的水质,须对该养殖场进行防渗处理或勒令关闭。      

稀有气体同位素地球化学应用前景

稀有气体有4个源区,包括大气、地下水、地壳和地幔;其中地幔稀有气体又可分为上地幔和下地幔2个源区。不同源区的稀有气体在同位素组成上截然不同,因而可以利用稀有气体同位素进行有关地质和地球化学过程的示踪。稀有气体同位素组成不仅对查明成矿流体的性质和来源及矿床成因具有重要作用,而且还可以为正确理解大规模成矿作用的地球动力学背景和全面弄清大型一超大型矿集区的时空分布规律提供重要线索。     

长江流域沿江镉异常示踪与追源的战略与战术

正在进行的多目标地球化学调查成果显示,长江流域存在全流域的Cd异常。长江流域Cd异常示踪与追源研究的长期目标是查明长江流域沿江各主要支流(汇水面积大于5000km^{2)Cd等重金属元素的物质来源、迁移形式和输入/输出通量、分辨自然源与人为源各自所占份额,建立沿江各支流Cd时空演化模型,监测它的未来发展趋势,对潜在生态效应进行预警预测;短期目标是针对Cd异常的重点地区,如长江源头、三峡库区、湘江流域、江淮流域、长江三角洲及流域内的4大淡水湖泊,查明Cd异常的来源,重建Cd异常形成的地球化学记录,评估可预见的将来(如10～50年)Cd异常的生态效应。}     

宜万铁路野三关隧道“8.5”突水事故成因分析

2007年8月5日宜(昌)万(州)铁路野三关隧道发生了重大的岩溶突水事故.在岩溶水文地质调查的基础上,对野三关隧道所处流域两个主要岩溶洼地以NaCl为示踪剂进行了充水水源和充水途径的示踪试验.两次示踪试验结果查明水洞坪岩溶洼地是"8.5"突水事故的主要水源,茅口组灰岩中发育的周家包暗河被F18断裂切割,将暗河水导入隧道形成突水事故.稻子坪的水在隧道突水段西以渗-涌水的形式进入隧道.水洞坪示踪实验得出二叠系栖霞-茅口组岩溶含水层中的地下水流速为50 m/h.稻子坪示踪试验得到茅口组灰岩含水层中的地下水流速为45.5 m/h;大冶组灰岩含水层中的地下水流速为10.8 m/h.本示踪结果为宜万铁路野三关隧道的后期恢复施工和防渗工作提供了依据,并为岩溶地区深埋隧道的设计和施工提供了试验和分析工作的实例参考.     

三维复杂介质的块状建模和试射射线追踪

为了解决三维复杂介质的射线追踪,本文改变了传统的三维层状地层的建模描述方式,提出了块状结构的建模描述方法,结合三角形面片来描述地质界面,可以构造非常复杂的三维地质模型．为了满足射线追踪的需要,本文对模型界面内的法向量进行光滑处理,光滑后的法向量在界面内是连续变化的．在块状模型的基础上,本文运用三角形的面积坐标,提出了几种试射角度的修正方法：细分三角形法、分割三角形法和子三角形法,计算表明子三角形法最好．文中给出了三维块状模型和射线追踪实例．     

三维复杂介质中转换波走时快速计算

复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.     

Traveling wavefront ray tracing and its applications

The current discussion of whether marine geodesy is something entirely new or part of oceanography or geodesy can be clarified by re‐examining our notions about the basic business of geodesy. This paper contends that the same basic geodetic services, which were needed for millenia to chart and control new territory, are being adapted now to the marine environment to suit modern accuracy standards. In a brief historical review it is shown that conceptually, the oceanic regions were always an integral part of geodetic concern. The oceans were certainly part of the world for the ancient map makers and for Pythagoras’ spherical earth model. Notions connected with the marine geoid were implied in ancient speculations. Distances and directions at sea, even depths, were determined to meet the requirements of the times. With the modern sophistication in geodetic theory and measuring techniques, these ideas became more refined and demanding in turn, but they were there all the time. For obvious reasons, land geodesy developed faster; marine geodesy will catch up now, because only now there is the technical capability and also an urgent economic motivation.