基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

褶皱陡坎中相关断层在缩短量计算中的作用——以东秋里塔格背斜为例

通过对东秋里塔格背斜地区内与褶皱陡坎伴生的次级断层进行的地质调查,确定了该次级断层属于伸出向斜的逆冲断层,也给出了次级断层的发生时间晚于褶皱作用起始时间的1个变形实例。褶皱陡坎中发育的次级断层使得陡坎上盘地层沿断层面整体向上迁移,不仅次级断层使得阶地面发生掀斜并增大褶皱陡坎的坡度,而且直接影响了区域缩短增量计算的真实性。计算结果显示,当不考虑次级断层对褶皱陡坎高度的影响时,计算所得缩短增量为51.42m,考虑次级断层对褶皱陡坎高度的影响时,计算得到的地壳缩短量为45.23m。二者相差6.19m,偏差占总缩短增量的13.7%,是一个不可忽视的量值。东秋里塔格背斜北翼和南翼的褶皱陡坎发育于相同岩性的基岩中而且具有相同的形成机制。但是背斜北翼在水平距离不足300m的范围内发育了3级褶皱陡坎,这表明北翼相对于南翼,其活动枢纽带更为紧闭。这是由于北翼存在更加强烈的挤压应力及更加快速的隆升作用。因此,次级断层的研究对于正确认识区域构造演化、了解褶皱与断层相互作用关系具有重要意义。但是仍存在很多问题:1)受观察剖面范围有限、次级断层分布不连续和断距沿延伸方向不断增大等因素的制约,次级断层对缩短增量造成的影响可能被低估,计算结果应为缩短量的最小值。2)次级断层增加的陡坎高度与断层的断距、倾角具有怎样的定量关系?3)若次级断层只发育于活动枢纽带内,又会产生怎样的影响?这些问题还需更多的研究实例来进行更加深入的研究。     

Surfactant Enhance DLLME/FO‐LADS: Assay of Malachite Green Level in Aquatic Environment of Trout Fish

Malachite green (MG), a traditional agent used in aquaculture although is not approved; its low cost and high efficacy make illicit use likely. We developed a small‐scale, simple, and sensitive dispersive liquid–liquid microextraction procedure for the assay of trace amounts of MG in aquatic environment of Trout fish. Fiber optic‐linear array detection spectrophotometry with charge‐coupled device detector benefiting from a microcell was used for this purpose. The method is based on enhancement effect of an anionic surfactant on the extraction of MG in to very fine multidroplets of microextraction solvent which made assisted by disperser solvent. Under the optimum conditions, the enrichment factor 77.5 was obtained from a 5‐mL water sample. The calibration graph was linear up to 5 × 10^?7 mol L^?1 with detection limit of 1 × 10^?8 mol L^?1. The relative standard deviation for seven replicate measurements of 4 × 10^?7 and 5 × 10^?8 mol L^?1 of MG were 3.3 and 4.5%, respectively.     

Growth by extension,and reworking,of a south‐western Kalahari linear dune

Optically stimulated luminescence (OSL) dating studies of linear (longitudinal) dunes have been used extensively to elucidate late Quaternary environments and climates in arid or formerly arid regions, yet understanding of the development of such dunes is incomplete. In particular, conflicting opinions have been presented regarding the propensity of linear dunes to migrate laterally, the degree to which they rework their own sediment during accumulation and whether they form primarily by extension, as opposed to lateral sand movement from adjacent interdunes. This study focuses on this last point, although the importance of the other controversies is discussed in context. A simple linear dune in the south‐western Kalahari, which has a prominent termination on a pan (playa) surface, provides an opportunity to directly test hypotheses of dune extension. Chronostratigraphy along a ~600 m transect along the crest of the dune, constrained by 42 OSL ages, reveals that the dune grew by extension on occasions in the late Pleistocene and early Holocene, but has also been subject to reworking along its length, which has continued until recent times. Dune development by extensional growth is suggested to operate under environmental conditions more conducive to net accumulation, whereas reworking is largely independent of conditions throughout the last ~18 ka, and may represent seasonal fluctuations in the position of the dune crest. The relative significance of these two modes of development is suggested to be a key control on the efficacy of linear dunes as archives of environmental and climatic change. Copyright © 2011 John Wiley & Sons, Ltd.     

位场DFT算法研究

本文推广了经典的抽样定理,并据此导出了函数有限离散傅里叶变换误差方程（简称DFT误差方程,下同）。该方程把有限离散傅里叶变换中固有的离散效应和有限效应表示为确切的数学形式。离散效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数;有限效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数的DFT。 基于DFT误差方程和位函数特点,作者提出了两种位场数值傅里叶变换新算法--移样法和等效源续尾叠样法。移样法可近百倍地提高位场数值傅里叶反变换的精度,等效源续尾叠样法可数十倍地提高正变换精度。两种算法都不需要增加资料长度和取样密度,因而基本不需要增加计算机时间和内存。文中给出了算例。     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

Asymptotic Rayleigh instantaneous unit hydrograph

The instantaneous unit hydrograph for a channel network under general linear routing and conditioned on the network magnitude,N, tends asymptotically, asN grows large, to a Rayleigh probability density function. This behavior is identical to that of the width function of the network, and is proven under the assumption that the network link configuration is topologically random and the link hydraulic and geometric properties are independent and identically distributed random variables. The asymptotic distribution depends only on a scale factor, , where is a mean link wave travel time.     

甚高频时岩石介电常数的测试方法和特性

本文讨论20-200MHz频率范围内平行板电容法测量岩石介电常数和电导率的方法,指出常用的静电场公式的局限性,给出较精确的似稳场动态公式,并讨论20-50MHz频率范围内的岩石介电特性。     

An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.