城市浅层地震勘探过障碍变观设计

在城市浅层地震勘探数据采集中,地震测线经常会遇到河流、桥、交叉路口等地表障碍。为了尽量减小地表障碍造成的影响,需要对观测系统作变观设计。作者分析总结了城市浅层地震勘探中常见障碍以及过障碍变观模式,并开发了浅层地震过障碍变观设计软件,利用此软件可在采集现场快捷地进行过障碍变观。将这种方法应用于实际工作中,有效减小了资料缺失造成的影响,提高了资料的信噪比。该软件能够根据实际工作环境和设备条件在野外现场对观测系统进行灵活变观,设计出合理的过障碍观测系统,有效地提高了变观设计的效率。     

Processes of initial collision and suturing between India and Asia

The initial collision between Indian and Asian continents marked the starting point for transformation of land-sea thermal contrast, uplift of the Tibet-Himalaya orogen, and climate change in Asia. In this paper, we review the published literatures from the past 30 years in order to draw consensus on the processes of initial collision and suturing that took place between the Indian and Asian plates. Following a comparison of the different methods that have been used to constrain the initial timing of collision, we propose that the tectono-sedimentary response in the peripheral foreland basin provides the most sensitive index of this event, and that paleomagnetism presents independent evidence as an alternative, reliable, and quantitative research method. In contrast to previous studies that have suggested collision between India and Asia started in Pakistan between ca. 55 Ma and 50 Ma and progressively closed eastwards, more recent researches have indicated that this major event first occurred in the center of the Yarlung Tsangpo suture zone (YTSZ) between ca. 65 Ma and 63 Ma and then spreading both eastwards and westwards. While continental collision is a complicated process, including the processes of deformation, sedimentation, metamorphism, and magmatism, different researchers have tended to define the nature of this event based on their own understanding, an intuitive bias that has meant that its initial timing has remained controversial for decades. Here, we recommend the use of reconstructions of each geological event within the orogenic evolution sequence as this will allow interpretation of collision timing on the basis of multidisciplinary methods.     

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     

宾川地区的地壳厚度和泊松比分布特征研究

This paper presents the changes of crust thickness and Poissons ratios distribution in the Binchuan region,where the first air-gun transmitting station and its a small dense array were deployed.From September 2011 to January 2014,more than 239 teleseismic events of M≥6.0 were recorded in 16 stations in the Binchuan region.Their P-wave receiver functions were analyzed respectively.The high spatial resolution result shows that the average crust thickness of Binchuan region is 45.3km,it follows the rule of"deeper in the north and east part,shallower in the south and west part."The deepest region is in Xiaoyindian Station;the crust thickness is 47.9km;the shallowest region is in Paiying Station,it has the thickness of 42.1km.It shows that the deeper Moho surface nearby the Chenghai fault and shallower nearby the Honghe fault;the isoline distribution of thickness changes greatly nearby the Chenghai fault and slowly nearby the Honghe fault.From the distribution of Poissons ratios,it is unevenly in the study area with a great difference from the north part to the south part,which shows a characteristic of"lower in the north,higher in the south".The Poissons ratio nearby the Honghe fault is medium too high(0.26≤σ≤0.29);lower nearby the Chenghai fault(σ≤0.26).This paper concludes the possible reason of different characteristic between Poissons ratio and crust thickness is thicker in the upper crust in the Binchuan region.     

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

高黎贡山地处印度板块与欧亚板块碰撞缝合带附近的横断山脉南段,是大理一瑞丽铁路(大瑞线)的必经之地,地形起伏大、构造复杂、活动性强,高黎贡山隧道作为全线控制性工程之一,其地质选线的最大困难就是对隧道深部构造环境的了解,特别是缺少对与地热、地震等联系紧密的深部地质构造的认识.为此,本文以大地电磁方法为手段,以高黎贡山隧道为主要研究对象,通过对滇西龙陵地区高黎贡山隧道越岭段两条大地电磁剖面数据的处理解释对研究区的地壳电性结构特征进行了勘探研究.结合区域地质构造特征与主要工程地质问题之间关系的分析,根据隧道主要断层地质条件设计了三维垂直断层模型,利用三维有限元开展正演模拟研究发现,测点点距、位置与横向分辨率密切相关,点距越密,分辨率越高,测点位于断层在地表投影位置能有效提高分辨率.采用大地电磁阻抗张量分解技术对两条剖面上各测点的二维偏离度和电性走向进行了计算和分析,对剖面视电阻率和阻抗相位数据进行了二维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波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

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.     

加勒比板块多学科研究综述及其对地震预测的启示

2010年1月12日海地M_W7.0地震打破了Enriquillo断层近240余年的地震平静期, 可能意味着加勒比地区新的活跃周期的到来。 本文通过收集该地区地质构造演化、 地球物理场、 发震构造、 震源物理、 运动学特征及应力场演化等方面的研究成果, 分析由于海地地震的发生所带来的区域地震危险性等相关认识的改变。 通过对多学科资料的梳理, 分析利用多学科研究成果开展加勒比地区地震危险性分析的基本思路, 以期为未来的中国大陆区域动力学模型构建、 地震预测预报研究提供基础参考。