电偶源瞬变电磁测深研究(六)——瞬变场资料的反演解释

阐述了用广义逆矩阵反演瞬变电磁测深资料的方法和实际效果,重点解决大偏移距电偶源瞬变电磁测深资料的反演。      

Numerical modeling and inverse parameter estimation of the large-scale mass movement Gradenbach in Carinthia (Austria)

This paper deals with the inverse problem of using time-displacement monitoring data to determine the material parameters of a numerical model of a large-scale mass movement. A finite element model for simulating the mechanical behavior is presented for the Gradenbach landslide in Carinthia, Austria. Particular attention is paid to the calibration of the constitutive relationships, which represent a prerequisite for a realistic quantitative analysis. After a short introduction to the concept of model-parameter identification, this paper demonstrates how to apply the proposed model identification strategy to determine model parameters for the Gradenbach example. The impact of the amount of reference data available for the inverse model-parameter analysis is evaluated by means of artificial reference data. Subsequently, the numerical model is calibrated using field measurement data. The results obtained are presented, and the benefits and drawbacks of the proposed concept are evaluated.     

Results of inversion of distorted magnetotelluric sounding curves (numerical experiment)

The paper presents the results of a joint inversion of magnetotelluric sounding (MTS) curves taken at several stations without preliminary selection and normalization and distorted by the presence of lateral electrical-conductivity inhomogeneities in the medium. In the calculations, we used synthetic MTS data for a three-dimensional model. Preparation and interpretation of data are carried out by the Trefftz method using a numerical model of the field and MTS curve distortions. To solve the inverse problem and optimize the subsurface model, we used a nonlinear least-squares method and an iterative process with calculation of the sensitivity matrix and its singular decomposition. The target functional is determined by the discrepancies between the model and synthetic experimental apparent-resistivity curves corresponding to the elements of the impedance tensors on the lateral diagonal. The reliability of the reconstructed subsurface model is characterized by the dispersion of the deviations of its parameters from the parameters of the known model used for the preparation of synthetic experimental data. The joint consideration of distorted apparent-resistivity curves at several stations increases the reliability of interpretation results. The obtained solution to the inverse problem is approximate and can be used as a starting model for more complex algorithms and programs.     

东喜马拉雅构造结岩体冷却的40Ar/39Ar年代学研究

对采自东喜马拉雅构造结核心地段雅鲁藏布大峡谷地区的13件标本中的20件矿物样品进行了系统的常规^40Ar／^39Ar年代学研究。数据显示,样品的（^40Ar／^39Ar）i值均接近尼尔值（295．5±5）,且绝大部分样品的坪年龄与其反等时线年龄在误差范围内一致。从数据统计结果来看,所测样品的^40Ar／^39Ar年龄大都集中在1．3Ma和2．5Ma左右,表明南迦巴瓦地区在上新世中期和更新世早期均经历了快速冷却抬升事件。本次测试的样品采自不同的高程及不同的构造单元,且样品原岩的成因及岩性各异,但沿着大峡谷由北向南不同地段的样品的不同矿物（角闪石、黑云母、白云母、钾长石）的^40Ar／^39Ar年龄相近,而同一样品中不同矿物的^40Ar／^39Ar年龄大小又并非完全按照矿物对氩同位素体系的封闭温度高低来分布,表明该地区在上新世以来的岩体冷却速率很大,以致该地区的矿物对氩同位素体系的封闭过程与处于缓慢冷却环境中的封闭过程明显不同。以本文报道的数据估算,南迦巴瓦地区的岩体在最近3Ma以来的冷却速率达120～240℃／Ma,岩体抬升速率达3．4—6．9mm／a。     

The application of an innovative inverse model for understanding and predicting landslide movements (Salazie cirque landslides,Reunion Island)

The prediction of landslide movement acceleration is a complex problem, among others identified for deep-seated landslides, and represents a crucial step for risk assessment. Within the scope of this problem, the objective of this paper is to explore a modelling method that enables the study of landslide function and facilitates displacement predictions based on a limited data set. An inverse modelling approach is proposed for predicting the temporal evolution of landslide movement based on rainfall and displacement velocities. Initially, the hydrogeology of the studied landslides was conceptualised based on correlative analyses. Subsequently, we applied an inverse model with a Gaussian-exponential transfer function to reproduce the displacements. This method was tested on the Grand Ilet (GI) and Mare-à-Poule-d’Eau (HB) landslides on Reunion Island in the Indian Ocean. We show that the behaviour of landslides can be modelled by inverse models with a bimodal transfer function using a Gaussian-exponential impulse response. The cumulative displacements over 7 years of modelling (2 years of calibration period for GI, and 4 years for HB) were reproduced with an RMSE above 0.9. The characteristics of the bimodal transfer function are directly related to the hydrogeological functioning demonstrated by the correlative analyses: the rapid reaction of a landslide can be associated with the effect of a preferential flow path on groundwater level variations. Thus, this study shows that the inverse model using a Gaussian-exponential transfer function is a powerful tool for predicting deep-seated landslide movements and for studying how they function. Beyond modelling displacements, our approach effectively demonstrates its ability to contribute relevant data for conceptualising the sliding mechanisms and hydrogeology of landslides.     

新疆天山(独山子)-西昆仑(泉水沟)地学断面地震与重力联合反演地壳构造特征

利用滑动平均的布格重力异常剖面结合地震资料和大地电磁资料,参考介质纵波速度与密度的经验转换关系,通过二维均匀介质模型人机交互式计算得到了新疆天山(独山子)-西昆仑(泉水沟)地学断面走廊域中心剖面的地壳构造特征.所得的反演结果表明,准噶尔盆地与北天山地体、中天山地体为一整体;准噶尔盆地和天山构造带的中、下地壳之间存在低密度层,但塔里木盆地不存在低密度层;由于受到南北两侧板块挤压,造成塔里木盆地的地壳结构在横向上似一南陡北缓、中部隆起的不对称"拱形"弯曲;塔里木盆地南缘与西昆仑构造带间呈"V型盆山耦合关系",其下方恰好是莫霍面上隆区.     

地震精确定位问题(综述)

本文综述了地震定位问题的历史沿革和最新进展。着重介绍了用在地震定位问题上的地球物理反演方法,尤其是相对定位方法的原理、优点以及存在的局限性。目的是提供从事这一领域研究工作的人员参考。     

A lithofacies model of upper jurassic productive deposits in the eastern part of the Shirotnoye Priobie region (Western Siberia)

On the basis of an integrated interpretation of well data, core analysis, geophysical well logging, and seismic data a new geological model of the JV₁¹⁻² stratum is suggested for one of the Shirotnoye Priobie areas in Western Siberia. The results of acoustic inverse and the rock elastic property research were used for the validation of lithofacies zonation. Supplementary exploration of the area and production drilling are recommended.     

Attenuation of multiple waves using the continuous wavelet transform (CWT)

This paper describes how the continuous wavelet transform is used to filter multiple waveforms in both time and frequency domains. It is well suited to process the stationary signals, and it shows the signal in both time and frequency scales. This new approach was tested first on synthetic data and then on real data. The results obtained on both cases were good. The method consists of identifying the multiples on which we apply a normal move out using the multiple velocity law. The multiples will be aligned and the primary reflections will not be aligned. This operation allows locating the multiples in the time-scale domain. We compute the continuous wavelet transform (CWT for short) in order to focus on the patterns relative to seismic events. To filter the multiples, we define a zone with frequency and time bounds. These bounds are deduced from the projection of the seismic trace. Then an automatic mask is applied to the pattern to be isolated. Filtering in time–frequency domain is done by keeping only the wavelet coefficients that are outside the mask and assigning zero to the coefficients larger than a threshold amplitude inside the defined zone. The mask shape does not matter, which is not the case in classical filtering, where both the window size and shape play a key role. The mask is defined from three parameters: time, frequency, and the wavelet coefficients. To go back to the time domain, one has to compute the wavelet transform inverse of the trace. This procedure is repeated for all traces. To reset the traces to their initial positions, we apply the dynamic correction inverse with the same velocity law as the multiples. It turns out that the attenuation of multiples by the CWT works fine, in particular, the two identified multiples were quasi eliminated (Fig. 10).     

Straightforward inversion scheme (SIS) for one-dimensional magnetotelluric data

This paper presents a Straightforward Inversion Scheme (SIS) for interpreting one-dimensional magnetotelluric sounding data. The basic steps of SIS are (i) parameterization of the layered model such that the layer thickness, expressed in units of its skin depth, is a constant (α); (ii) expansion of the reflection function at each interface as a power series in parameter u = exp(-2(1 +j)α√f);(iii) development of a recurrence relation between the coefficients of the same powers ofu in the power series of reflection functions of any two successive layers; (iv) estimation of the impedance power series coefficients using regressed minimum norm estimator; and (v) evaluation of layer resistivities and thicknesses using the inverse recurrence relation. The power of SIS is established by inverting four synthetic data sets and two field data sets. The effect of noise is extensively studied on a synthetic data set, deliberately corrupted with increasing levels of Gaussian random noise up to 25%. It is found that the scheme can retrieve broad features of the true model even with noise levels as high as 25%. On the basis of findings of different experiments conducted on SIS, it is concluded that SIS is an efficient, robust algorithm with high resolving power. Further, being linear, it is non-iterative and it dispenses with the requirement of having to choose an initial guess model.     

Partial melting below Tubuai (Austral Islands, French Polynesia)

Process identification diagrams based on trace element data show that mafic lavas from Tubuai, including alkali basalts, basanites, analcitites and nephelinites, result from different degrees of partial melting of an isotopically homogeneous mantle source. Our fractionation-corrected data are consistent with a batch melting model or a dynamic melting model involving a threshold value for melt separation close to 1% and degrees of melting ranging from 5–8% (alkali basalts) to 1.5–3% (nephelinites). The relative source concentration pattern, calculated using an inverse numerical method, shows an enrichment in highly incompatible elements. We propose that the Tubuai lava suite was derived from a two-stage partial melting process. Melting first affected the plume material located within the transition zone between garnet and spinel domains, producing alkali basalts and basanites. Then, the melting zone migrated upwards to the base of the overlying spinel-bearing lithospheric mantle, producing highly silica-undersaturated lavas. The lower lithosphere had previously been enriched by intrusion of pyroxenite veins representing plume-derived melts which percolated away from the main magma conduits. Received: 11 June 1996 / Accepted: 8 January 1997     

逆数据域表面多次波压制方法

基于波动方程压制表面多次波的反馈迭代方法,推导出压制表面多次波的逆数据域方法。该方法利用了一次波和表面多次波在逆数据域的可分离性质。含表面多次波的地震数据变换到逆数据域后,表面多次波的能量在逆数据域形成能量聚焦区域,可通过简单切除能量聚焦区域的信号实现表面多次波压制,避免了基于波动方程预测减去方法中的自适应相减过程,完全损伤不到一次波的能量。单道地震记录和合成炮数据模型试算表明：逆数据域多次被压制方法简单易行,可有效压制表面多次波,并保持一次波振幅不受损失。与常规的基于波动方程的表面多次被压制方法（SRME方法）对比表明,逆数据域方法能更好地应对有效波和多次波同相轴相交的情况。     

岩体水力学基础（七）：岩体水力学参数的确定方法

本文在对参数的内涵分析的基础上,提出了参数的类型及其概念。介绍了岩体渗透系数张量确定的统计学方法、野外抽（压）水试验方法、数学模型反演求解方法,以及应力与渗透系数关系参数的确定方法。     

Lithospheric Refertilization Trends in Xenoliths and Xenocrysts from the Udachnaya Kimberlite (Siberian Craton)

Doklady Earth Sciences - Comprehensive studies of peridotitic xenoliths from the Udachnaya kimberlite (Yakutian diamond province, Siberian craton) confirm that garnet shows inverse correlation of...     

稳定有效的反Q滤波方法

反Q滤波是补偿大地吸收衰减有效方法之一,其稳定性和有效性是反Q滤波的前提条件.笔者研究了一种稳定的反Q滤波方法,并通过理论模型和实际数据与传统的增益限制反Q滤波方法进行了对比.结果表明:稳定反Q滤波更好地恢复了衰减的地震信号,并自动限制了对噪声的过渡放大.同时,还分析了噪声和增益限制因子G_(lim)对稳定反Q滤波的影响.     

Modeling water resources of a highly irrigated alluvial plain (Italy): calibrating soil and groundwater models

Modern and effective water management in large alluvial plains that have intensive agricultural activity requires the integrated modeling of soil and groundwater. The models should be complex enough to properly simulate several, often non-linear, processes, but simple enough to be effectively calibrated with the available data. An operative, practical approach to calibration is proposed, based on three main aspects. First, the coupling of two models built on well-validated algorithms, to simulate (1) the irrigation system and the soil water balance in the unsaturated zone and (2) the groundwater flow. Second, the solution of the inverse problem of groundwater hydrology with the comparison model method to calibrate the model. Third, the use of appropriate criteria and cross-checks (comparison of the calibration results and of the model outputs with hydraulic and hydrogeological data) to choose the final parameter sets that warrant the physical coherence of the model. The approach has been tested by application to a large and intensively irrigated alluvial basin in northern Italy.     

灰色系统GM(n,m)模型参数的递推估计和时变估计

本文给出灰色系统GM(n,m)模型参数的递推最小二乘估计和时变参数估计,前者的优点是可以避免反复求逆矩阵,因而使计算简便,同时克服了文献1中给出的估计方法当观察数据多、矩阵线数高、分配单元增加时,微型机容量不够的弱点,而且随着新数据的到来,通过递推算法自动对系统参数加以修正,提高了模型的自适应性;时变参数估计可以缩小系统的时变性与模型参数非时变性所产生的差异。     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

改进的等参逆变换算法在耦合场分析中的应用

实际工程有限元计算中涉及多个物理场的相关分析或耦合分析时场变量之间的数据传递问题非常重要,往往会遇到等参逆变换的需求。分析并指出现有等参逆变换所存在问题,给出了新的改进算法,即在网格内采用新方法进行等参逆变换,在网格外采用距离加权法进行插值传递,并详细说明本算法实施关键。对一大型拱坝施工过程中的温度场和应力应变场耦合分析问题进行研究,结果表明所提出方法具有简便实用、高效和满足工程精度等特点,可应用于实际工程。