二次电流场多次叠加概率成像

基于地下电流场的积分公式,离散化的二次电流场被分解为电性不连续界面上的一系列点电荷电场的代数和.单位正点电荷电场被引入作为空间扫描函数（SDS）,积累电荷出现的概率（COP）函数定义为二次电流场与SDS的互相关.为对概率成像结果进行定量分析解释,提出了规范的积累电荷出现的概率（NCOP）函数.通过应用有限元算法对2D地电模型进行二次电流场合成,实现了二次电流场的多次叠加概率成像.结果表明对均匀半空间中赋存地质异常体的电性结构,概率成像方法对地下异常体的空间位置有较好的指示作用.     

Control of sediment dynamics by vegetation as a key function driving biogeomorphic succession within fluvial corridors

Riparian vegetation responds to hydrogeomorphic disturbances and environmental changes and also controls these changes. Here, we propose that the control of sediment erosion and deposition by riparian vegetation is a key geomorphological and ecological (i.e. biogeomorphic) function within fluvial corridors. In a 3 year study, we investigated the correlations between riparian vegetation and hydrogeomorphic dynamics along a transverse gradient from the main channel to the floodplain of the River Tech, France. Sediment erosion and deposition rates varied significantly along the transverse gradient as a function of the vegetation biovolume intercepting water flow. These effects, combined with the extremely strong mechanical resistance of pioneer woody structures and strong resilience of pioneer labile herbaceous communities, Populus nigra and Salix spp., explain the propensity of biogeomorphic succession (i.e. the synergy between vegetation succession and landform construction) to progress between destructive floods. This geomorphological function newly identified as an ‘ecosystem function’ per se encompasses the coupling of habitat and landform creation, maintenance and change with fundamental ecosystem structural changes in space and in time. Three different biogeomorphic functions, all related to the concept of ecosystem engineering, were identified: (i) the function of pioneer herbaceous communities to retain fine sediment and diaspores in the exposed zones of the active tract near the water resource, facilitating recruitment of further herbaceous and Salicacea species; (ii) the function of woody vegetation to drive the construction of forested islands and floodplains; and (iii) the function of stabilised riparian forests to act as ‘diversity reservoirs’ which can support regeneration after destructive floods. Overall, this study based on empirical data points to the fundamental importance of sediment flow control by pioneer riparian vegetation in defining fluvial ecosystem and landform organisation in time and in space. Copyright © 2009 John Wiley & Sons, Ltd.     

川西地区壳幔结构与汶川MS8.0级地震的孕震背景

收集了四川数字地震台网记录的57个远震事件,并从宽频带数字化三分量地震记录中计算出了马尔康（MEK）、都江堰（YZP）、中江（JJS）、江油（ZJG）、广元（YTS）、康定（GDS）、汉源（XJP）、雅安（MDS）、峨眉山（EMS）、沐川（WMP）、仁寿（YGD）、荣县（HMS）等12个台站下方的远震P波径向接收函数.另外,引入地震勘探中的动校正技术,将各台的接收函数校正到67°的参考震中距处,然后对接收函数进行叠加以增强信号,并把叠加接收函数作为台站下方的平均接收函数.最后,利用台站下的平均接收函数反演得到S波速度结构.反演结果表明：以锦屏山—龙门山断裂为界,其西侧地壳厚达70 km,而东侧仅为50 km左右,Moho面在断裂下方形成了一个陡坎;川西地区的地壳速度结构与川中地区差异较大,主要表现在川西地区的中地壳存在厚度为8~22 km的低速层.在都江堰、雅安一带,低速层的厚度最大,其厚度在20~22 km之间,其上地壳为一个坚硬固体,在区域构造应力场作用下,形成了孕育大地震的构造环境.     

接收函数方法估计Moho倾斜地区的地壳速度比

造山带地区的复合力系作用往往使Moho界面发生变形,局部表现为倾斜状态.为了得到这些区域精确地壳速度比结构,本文基于H-κ方法发展了H-κ-θ方法.该方法不仅考虑了倾斜Moho层的响应,同时利用径向和切向接收函数信息,增加了对扫描的约束.利用该方法对青藏高原东南缘地壳厚度和速度比结构进行研究,结果表明：研究区内地壳明显存在不均匀性,松潘—甘孜地体平均地壳厚度约为60 km,四川盆地西缘约为47 km,扬子地台约为43 km,三江块体和扬子地台东南缘已接近正常地壳厚度;松潘—甘孜地体与扬子地台相邻部位地壳平均地震波速度比（V_p/V_s）普遍偏高,且四川盆地西侧发现一绕盆地边缘的弧形高V_p/V_s异常区（>1.88）,我们推测该异常可能由青藏高原向东逃逸的地壳流体受到高强度的四川盆地阻挡,在其西侧堆积所致.     

提取视震源时间函数的PLD方法及其对2005年克什米尔MW7.6地震的应用

通过重构用于确定视震源时间函数有效持续时间的判别函数,对提取视震源时间函数的PLD方法进行了改进;利用合成资料和实际资料,验证了改进后PLD方法的可行性和稳定性.将PLD方法应用于2005年克什米尔M_W7.6地震及其11个余震的1887条记录,在84个台站处获得了这次地震的视震源时间函数.分别平均从不同台站的P波、S波、Rayleigh波和 Love波中得到的视震源时间函数,获取了主震的平均视震源时间函数.对视震源时间函数的分析表明,2005年克什米尔M_W7.6地震的持续时间大约为25 s,这是一次“急始型”地震,总体上表现为圆盘形破裂.但有迹象表明,破裂在初期有向西北方向发展的单侧传播趋势.     

基于精确权函数的仪器响应失真校正方法与算例

当前,经仪器响应失真校正后,加速度计有效观测频带通常为0至1倍自振频率,这限制了强震仪应用领域的延拓。有鉴于此,本文提出一个基于精确权函数的仪器响应失真校正方法并给出一个验证算例。分析表明:较之现有方法,新方法兼有高精度和高效率的优点,并可以将力平衡式加速度计的有效观测频带至少延拓3倍,对于拓展强震仪的应用领域具有重要意义。     

Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free‐falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of 1 m × 1 m, the expected number of received free‐falling raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re‐detachment amount. The re‐detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free‐falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re‐detachment amount were small parts of the total splash amount. Their proportions were 0·15% and 2·6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil‐splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright © 2007 John Wiley & Sons, Ltd.     

井间电磁测量的2．5维层析成像方法

利用正则化最小二乘反演方法实现了井间电磁测量数据的层析成像,对井间地层电阻率进行了重建。在成像算法中,我们假设了井间电磁的激发与接收采用电磁偶极子源,井间介质仅在二维（xoz）平面内变化。在数值模拟中,通过对构造走向（y方向）的Fourier变换,将三维电磁场问题转化为一系列二维问题,用等参有限元方法在波数域求解,使实际地层模型的处理得以实现。对于波数域中每个波数对应的电磁场方程采用等参有限元求解,并用高斯积分将波数域解变换为空间域电磁场。利用源与接收器电磁场的互易原理,实现了电磁场响应对电导率分布灵敏度的快速计算。针对正演模拟中源点的奇异性,我们采用具有一定面积的伪艿函数表达源电流分布,使数值解精度得到提高。用层状介质的解析解与数值计算结果的对比,验证了模拟算法的精度。用介质扰动产生的电磁场变化检验互易性定理计算灵敏度的有效性。对简单块状模型、斜向裂缝带模型及“大”字模型的模拟数据成像结果表明,本文介绍的层析成像方法是正确有效的。     

S-wave crustal and upper mantle's velocity structure in the eastern Tibetan Plateau——Deep environment of lower crustal flow

A teleseismic profile consisting of 26 stations was deployed along 30°N latitude in the eastern Tibetan Plateau. By use of the inversion of P-wave receiver function, the S-wave velocity structures at depth from surface to 80 km beneath the profile have been determined. The inversion results reveal that there is significant lateral variation of the crustal structure between the tectonic blocks on the profile. From Linzhi north of the eastern Himalayan Syntaxis, the crust is gradually thickened in NE direction; the crustal thickness reaches to the maximum value (~72 km) at the Bangong-Nujiang suture, and then decreased to 65 km in the Qiangtang block, to 57―64 km in the Bayan Har block, and to 40―45 km in the Sichuan Basin. The eastern segment of the teleseismic profile (to the east of Batang) coincides geographically with the Zhubalong-Zizhong deep seismic sounding profile carried out in 2000, and the S-wave velocity structure determined from receiver functions is consistent with the P-wave velocity structure obtained by deep seismic sounding in respect of the depths of Moho and major crustal interfaces. In the Qiangtang and the Bayan Har blocks, the lower velocity layer is widespread in the lower crust (at depth of 30―60 km) along the profile, while there is a normal velocity distribution in lower crust in the Sichuan Basin. On an average, the crustal velocity ratio (Poisson ratio) in tectonic blocks on the profile is 1.73 (σ = 0.247) in the Lhasa block, 1.78 (σ = 0.269) in the Banggong-Nujiang suture, 1.80 (σ = 0.275) in the Qiangtang block, 1.86 (σ = 0.294) in the Bayan Har blocks, and 1.77 (σ = 0.265) in the Yangtze block, respectively. The Qiangtang and the Bayan Har blocks are characterized by lower S-wave velocity anomaly in lower crust, complicated Moho transition, and higher crustal Poisson ratio, indicating that there is a hot and weak medium in lower crust. These are considered as the deep environment of lower crustal flow in the eastern Tibetan Plateau. Flowage of the ductile material in lower crust may be attributable to the variation of the gravitational potential energy in upper crust from higher on the plateau to lower off plateau.     

Statistical model on the surface elevation of waves with breaking

In the surface wind drift layer with constant momentum flux, two sets of the consistent surface eleva- tion expressions with breaking and occurrence conditions for breaking are deduced from the first in- tegrals of the energy and vortex variations and the kinetic and mathematic breaking criterions, then the expression of the surface elevation with wave breaking is established by using the Heaviside function. On the basis of the form of the sea surface elevation with wave breaking and the understanding of small slope sea waves, a triple composite function of real sea waves is presented including the func- tions for the breaking, weak-nonlinear and basic waves. The expression of the triple composite func- tion and the normal distribution of basic waves are the expected theoretical model for surface elevation statistics.