面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

长江中下游浅水湖泊悬浮颗粒物吸收系数测定方法探讨

对两种水体悬浮颗粒物吸收系数测定方法及相关计算进行对比研究.通过长江中下游湖泊典型藻类的实验室培养,利用T方法和T-R方法分别对藻类颗粒物、藻类泥沙混合悬浊液进行吸收系数测定.通过颗粒物光谱吸收系数与叶绿素a之间的相关性关系,对比了两种方法的测量稳定性.通过对不同比例的藻类和无机悬浮颗粒物(ISS)的混合悬浊液进行分析,获得了不同浊度水体悬浮物吸收光谱的变化情况.结果表明,在纯藻或者泥沙含量较少的水体进行颗粒物吸收系数光谱测定时,T方法和T-R方法均可以采用,并且均具有较高的测定精度.然而,在泥沙含量相对较高的浑浊水体,应尽量选取T-R方法进行颗粒物吸收光谱的测定,以提高测定精度.长江中下游浅水湖泊由于底泥易受风浪影响发生再悬浮,因此在颗粒物吸收系数光谱测定中,当水体中ISS含量超过30 mg/L时,应选择T-R方法.     

CHARACTERISTICS OF FOCAL MECHANISMS AND STRESS FIELD OF YUNNAN AREA

Using the seismic waveform data recorded by regional seismic network of Yunnan and Sichuan and the method of CAP, we calculate and obtain the focal mechanism of 268 earthquakes with the magnitude of M_L≥4.0 occurring in Yunnan during Jan. 1999 to Aug. 2014; then, we analyze the types and the regional feature of the focal mechanism of earthquakes in Yunnan, on the basis of the focal mechanism of 109 earthquakes analyzed by Harvard University. Based on the data of the above focal mechanism solutions, we adopt the method of damped regional-scale stress inversion to calculate the best-fitting tectonic stress tensor of every grid in Yunnan; and adopt the method of maximum principal stress to calculate the direction of maximum horizontal principal stress in Yunnan. The result shows that: (1)the strike-slip type is the most principal type of the earthquake focus in the study area and the second is the normal faulting type; while, the reverse-fault type is relatively small. The spatial distribution of focal mechanism is obvious. This reflects that the dynamic source and acting force are different in different parts of the study area. (2)The direction of the stress field in Yunnan shows a certain spatial continuity. Maximum horizontal principal compressive stress is mainly clockwise from north to south and counterclockwise from the west to the east. The direction of stress field shows inhomogeneity in space. There exist two stress conversion zones respectively in EW and NS direction. The inversion result of stress field shows that the stress field in Yunnan is complex and the principal stress direction changes greatly; and there are obvious differences in different regions.     

Development of a distributed hydrological model to facilitate watershed management

To facilitate precise and cost-effective watershed management, a simple yet spatially and temporally distributed hydrological model (DHM-WM) was developed. The DHM-WM is based on the Mishra-Singh version of the curve number method, with several modifications: The spatial distribution of soil moisture was considered in moisture updating; the travel time of surface runoff was calculated on a grid cell basis for routing; a simple tile flow module was included as an option. The DHM-WM was tested on a tile-drained agricultural watershed in Indiana, USA. The model with the tile flow module performed well in the study area, providing a balanced water budget and reasonable flow partitioning. The daily coefficient of determination and Nash-Sutcliffe coefficient were 0.58 and 0.56, for the calibration period, and 0.63 and 0.62 for the validation period. The DHM-WM also provides detailed information about the source areas of flow components, the travel time and pathways of surface runoff.

EDITOR A. Castellarin; ASSOCIATE EDITOR F.-J. Chang     

滇西龙陵地区地壳电性结构及其对大瑞铁路地质选线影响研究

高黎贡山地处印度板块与欧亚板块碰撞缝合带附近的横断山脉南段,是大理一瑞丽铁路(大瑞线)的必经之地,地形起伏大、构造复杂、活动性强,高黎贡山隧道作为全线控制性工程之一,其地质选线的最大困难就是对隧道深部构造环境的了解,特别是缺少对与地热、地震等联系紧密的深部地质构造的认识.为此,本文以大地电磁方法为手段,以高黎贡山隧道为主要研究对象,通过对滇西龙陵地区高黎贡山隧道越岭段两条大地电磁剖面数据的处理解释对研究区的地壳电性结构特征进行了勘探研究.结合区域地质构造特征与主要工程地质问题之间关系的分析,根据隧道主要断层地质条件设计了三维垂直断层模型,利用三维有限元开展正演模拟研究发现,测点点距、位置与横向分辨率密切相关,点距越密,分辨率越高,测点位于断层在地表投影位置能有效提高分辨率.采用大地电磁阻抗张量分解技术对两条剖面上各测点的二维偏离度和电性走向进行了计算和分析,对剖面视电阻率和阻抗相位数据进行了二维NLCG联合反演研究,揭示了沿剖面的腾冲地块、龙陵一瑞丽断裂带及保山地块10 km深度的电性结构特征及相互关系.结果表明:剖面CD电性结构呈现区域构造的三分性,腾冲地块电性结构成层性较好,保山地块成层性较差,两者均以中高阻电性特征为主,中间夹龙陵一瑞丽断裂带,电性结构反映从3 km深度以下存在几乎近于直立延伸的低阻带,推测为班公湖一怒江缝合带滇西段丁青一怒江缝合带的反映;剖面AB共划分了6条与工程密切相关的深部隐伏断裂,结合地震地质、地表地质及龙陵地震深部背景研究,推测F7-3断裂为1975龙陵7.3级地震断裂;从地表黄草坝断裂开始向下延伸,有一条发育最大深度约为4 km的低阻通道,推测为地热断裂深循环通道,其与黄草坝断裂共同控制研究区地下热水的补给、径流和排泄条件,在高黎贡山隧道线位位置形成了一个相对低温通道,为隧道方案成立的关键工程地质条件.勘探结果表明:滇西龙陵地区地壳电性结构有效的反映了高黎贡山隧道深部隐伏断裂和地热断裂深循环通道等深部构造特征,为大瑞线隧道工程地质选线提供了深部地质背景依据.     

基于动态图像变形的PP与PS波层位直接匹配

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分:首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.     

A computational complexity analysis of Tienstra’s solution to equality-constrained adjustment

Tienstra’s method was developed to solve parameter adjustment with linear equality constraints, which has been otherwise often carried out by directly applying the conventional (or standard) method of Lagrange multipliers to quadratic optimization problems with a positive definite matrix. We analyze the computational complexity of the celebrated Tienstra’s method and compare it with that of the method of Lagrange multipliers. We show that Tienstra’s method is not only statistically elegant but also interestingly of significant computational advantage over the method of Lagrange multipliers to solve weighted least squares adjustment with linear equality constraints, with the saving of computational cost by a minimum of about 38% to a maximum of 100%. Tienstra’s method can be very important for large scale problems of adjustment and inversion with linear equality constraints.     

波浪作用下锚链系泊浮式防波堤动力响应的数值模型研究

浮式防波堤充分利用波能在水深方向的分布特性,在满足工程消浪要求的同时对海域水沙交换影响较小,且能够快速布置,在某些实际工程有一定应用前景。为了深入了解波浪作用下浮式防波堤的动力响应,基于OpenFOAM标准求解器olaFlow,在刚体运动求解计算中植入锚链求解模块MOODY(mooring cable dynamics),实现了基于重叠网格方法的浮体运动与锚链受力耦合求解,建立了锚链系泊浮式防波堤动力响应的二维数值模型。利用该数值模型对锚链系泊单方箱浮式防波堤在波浪作用下的透射系数、运动响应、锚链张力进行了模拟,并和相关试验结果进行了比较。结果表明,模型能够准确模拟二维波浪和浮式防波堤的相互作用,并用于三维模型的改进。     

基于EMD-AR谱的3C 273光变周期分析

光变是Blazar最显著的观测特征之一,有的光变呈准周期性. Blazar周期光变可能为其中心黑洞和吸积盘等内部物理结构提供重要信息.因此,周期分析对于讨论Blazar的光变有重要意义.结合经验模态分解(Empirical Mode Decomposition, EMD)和自回归(Auto-regressive, AR)模型谱估计方法的优点,提出了一种基于EMD-AR谱的组合光变周期分析新方法.先对观测数据进行EMD分解,得到各阶模态分量,并计算其与原始光变曲线的相关系数,再取相关度高的分量求和后估算AR谱,最后将此功率谱的周期分析方法用于天体的光变数据分析.讨论了EMD-AR谱方法的原理和应用步骤,并用此方法对类星体3C 273在1887—2016年的观测数据进行了分析,得到3C 273的长时标周期分别为:21.23、13.51、11.02、5.51、4.69、3.79、2.76 yr,与文献中常用方法的分析结果有很好的一致性;短时标周期分别为:(30±1)、(15±0.3)、(7.5±0.2)、(10±0.1)、(5±0.6)、(6±0.4)和(3±0.5) min,此结果未见相关...