利用位场及其高阶导数的功率谱计算岩层的深度与厚度

本文导出了利用重力与磁异常及其高阶导数的功率谱计算密度体和磁性体上顶与下底埋深的近似公式,还提供了利用垂向或水平导数谱计算深度的方法。理论模型的试验结果表明,由近似公式计算的精度可以满足实际需要。本文最后根据一条重力剖面资料,计算了密度层上、下界面。     

有限长圆柱体磁异常场全空间正演方法

在经典位场理论中,许多简单形体位场异常难以通过积分得到全空间的解析式.圆柱体是一类很重要的理论模型体,常用于模拟圆柱状地质体或非地质体(如管线),但目前还不能用解析公式正演有限长圆柱体在三维空间里的磁异常,而多是采用近似简化为有限长磁偶极子或线模型代替.对于有限长圆柱体,特别是半径相对于上顶埋深较大时,这种近似的误差不可忽略.本文利用共轭复数变量替换法,推导出有限长圆柱体在全空间的引力位一阶、二阶导数,利用Poisson关系得到磁异常正演公式,进而利用有限长圆柱体磁异常正演公式求解管状体的磁异常,得到不同磁化方向、不同大小的管线产生的磁场的特征,并将其推广到截面为椭圆的情况.最后通过模拟计算定量给出了将圆柱体近似为线模型的条件.     

Excitation of hide waves in the core due to transient-spatially periodic stress in the lower mantle

Summary Forced hydromagnetic oscillations in the outer core due to stress at its top are studied by means of a simple linear theoretical model originally proposed byHide [1] for the free oscillation studies. Both the magnetic and the inertial modes are excited. The expressions for the velocity, magnetic and vorticity fields are presented for a variety of stresses. They depend upon the wave number and frequency. In a numerical example wherein the stress varies periodically both in the space and the time, the amplitude of the induced magnetic field is found to be significant for reasonable values of various parameters.     

多层等效源曲面磁异常转换方法

在磁异常数据处理中,利用等效源技术重构磁异常场具有较好的稳定性和较高的计算精度,因而被广泛应用.传统方法是采用设置在近地表的单层等效源拟合实测磁异常数据,尽管拟合精度很高,但向上延拓之后往往会出现较大的拟合误差,即可能存在磁异常信号的"泄漏",尤其在原始数据中存在背景场时更容易出现此种误差.本文提出一种多层等效源技术方案,应用分布于不同深度范围内的等效源模拟实测数据,减少了等效源参数设置的盲目性.理论模型试验表明,采用多层等效源方法重构的磁异常及其梯度与分量,较单层等效源方法具有更高精度,可以吸收更完整的实测磁异常信息.论文详细地讨论了如何优化多层等效源设置、等效源参数选择以及计算方法,通过二维和三维理论模型试验,验证了在复杂条件下多层等效源方法的可行性和适应性,并且将该方法应用于广西某地的实测磁异常数据转换之中,取得了较好的应用效果.     

Magneto-rheological elastomer (MRE) based composite structures for micro-vibration control

Magneto-rheological elastomers (MREs) are used to construct composite structures for micro-vibration control of equipment under stochastic support-motion excitations. The dynamic behavior of MREs as a smart viscoelastic material is characterized by a complex modulus dependent on vibration frequency and controllable by external magnetic fields. Frequency-domain solution methods for stochastic micro-vibration response analysis of the MRE-based structural systems are developed to derive the system frequency-response function matrices and the expressions of the velocity response spectrum. With these equations, the root-mean-square (RMS) velocity responses in terms of the one-third octave frequency band spectrum can be calculated. Further, the optimization problem of the complex moduli of the MRE cores is defined by minimizing the velocity response spectra and the RMS velocity responses through altering the applied magnetic fields. Simulation results illustrate the influences of MRE parameters on the RMS velocity responses and the high response reduction capacities of the MRE-based structures. In addition, the developed frequency-domain analysis methods are applicable to sandwich beam structures with arbitrary cores characterized by complex shear moduli under stochastic excitations described by power spectral density functions, and are valid for a wide frequency range.     

Estimation of lightning and sprite parameters based on observation of sprite-producing lightning power spectra

Spectrograms and ELF power spectra of magnetic variations originated from sprite-producing lightning discharges have been analyzed to extract both parent lightning and sprite parameters. Some of the spectrograms and power spectra have been found to have approximately quasi-oscillatory shape in the frequency range 0–40 Hz with maximum repetition period about 15–20 Hz. A theory predicts that this interesting peculiarity of the power spectra can be due to interference between electromagnetic fields originated from the parent lightning discharge and from the sprite. A smooth envelope of the power spectrum was shown to have a form of damped oscillations with period close to reciprocal value of sprite lag time. A technique of extracting sprite parameters based on the sprite-producing lightning power spectrum is proposed. The lack of the first Schumann resonance and other features occasionally observed in spectral resonance structure were also discussed.     

基于奇异谱分析的重磁位场分离方法

奇异谱分析是一种近年兴起的时间序列分析方法,它利用降秩原理实现信号分离.该方法将数据空间投影到不同特征的子空间中,并用奇异值来表征这些子空间的性质,最后通过截取奇异值实现数据的重构.重磁位场分离可以看成一种多信号叠加的分离问题.不同特征的重磁异常具有不同特征的奇异谱,这是奇异谱分析用于解决位场分离问题的应用基础.本文通过建立理论模型,分析重磁异常的奇异谱特征,得出适用于重磁位场分离的最优参数选择方法,并与传统方法进行比较.对比发现,无论是横向叠加模型、垂向叠加模型还是斜向叠加模型,奇异谱分析都具有很好的分离效果.最后,将奇异谱分析用于鄂东南某矿区的重力资料处理中,实现弱异常的识别和分离.     

Multi-transmitter multi-receiver null coupled systems for inductive detection and characterization of metallic objects

Equivalent dipole polarizabilities are a succinct way to summarize the inductive response of an isolated conductive body at distances greater than the scale of the body. Their estimation requires measurement of secondary magnetic fields due to currents induced in the body by time varying magnetic fields in at least three linearly independent (e.g., orthogonal) directions. Secondary fields due to an object are typically orders of magnitude smaller than the primary inducing fields near the primary field sources (transmitters). Receiver coils may be oriented orthogonal to primary fields from one or two transmitters, nulling their response to those fields, but simultaneously nulling to fields of additional transmitters is problematic. If transmitter coils are constructed symmetrically with respect to inversion in a point, their magnetic fields are symmetric with respect to that point. If receiver coils are operated in pairs symmetric with respect to inversion in the same point, then their differenced output is insensitive to the primary fields of any symmetrically constructed transmitters, allowing nulling to three (or more) transmitters. With a sufficient number of receivers pairs, object equivalent dipole polarizabilities can be estimated in situ from measurements at a single instrument sitting, eliminating effects of inaccurate instrument location on polarizability estimates. The method is illustrated with data from a multi-transmitter multi-receiver system with primary field nulling through differenced receiver pairs, interpreted in terms of principal equivalent dipole polarizabilities as a function of time.     

位场向下延拓的改进迭代维纳滤波法

根据维纳滤波理论导出的位场向下延拓滤波器为最佳下延滤波器,但因其实现需要已知待求位场和噪声的功率谱而在实际应用中受到限制.针对该问题,本文首先提出一种基于位场径向平均功率谱的位场噪声水平估计方法,进而利用偏差准则求取正则化参数,实现位场正则化向下延拓;然后将位场正则化下延结果的功率谱作为待求位场功率谱的估计初值,采用带修正项的迭代维纳滤波方法来更新对待求位场功率谱的估计,最后提出本文的位场向下延拓改进迭代维纳滤波方法.基于理论重力模型数据及航磁实测数据进行了向下延拓对比试验,结果表明,改进迭代法具有较好的收敛性,且下延精度优于Tikhonov正则化法和递增型维纳滤波法.     

LINEAR TRANSFORMATIONS OF GRAVITY AND MAGNETIC FIELDS*

The spectral representation of gravity and magnetic fields shows that the mathematical expressions describing these fields are the result of convolution of factors which depend on the geometry of the causative body, the physical properties of the body and the type of field being observed. If a field is known, it is possible to remove or alter these factors to map other fields or physical parameters which are linearly related to the observed field. The transformations possible are: continuation, reduction to the pole, converting between gravity and magnetic fields, converting between components of measurement, calculation of derivatives, and mapping magnetization and density distribution, relief on interfaces, and vertical thicknesses of layers.     

等效源法三维随机点位场数据处理和转换

为了实现曲面随机点位场数据的曲面延拓和转换,以磁异常位场数据为例,采用一组磁偶极子作为等效源,置于观测面下方的一个曲面上,把观测磁异常作为这组磁偶极子所产生磁异常的边界条件,通过求解线性方程组的方法反演磁偶极子磁矩的大小,再根据反演结果正演所要计算的磁异常数据,实现了曲面随机点磁异常位场数据的向上延拓、向下延拓、求导以及化极处理.在数据量较大时,为了提高反演计算的速度,把磁异常数据和磁偶极子分成若干小块,再利用各块磁异常数据分别反演该块数据下方磁偶极子的磁矩,并通过迭代计算来逐步取得更准确的反演结果.模型试验表明,磁异常位场数据向上延拓的均方根误差小于±2nT,向下延拓和化极也可以取得较高的精度,所提出的分块处理方法提高了延拓和转换的速度,实际资料处理给出了曲面随机点航磁异常数据向下延拓和化极的一个例子.     

INTERPRETATION OF MAGNETIC ANOMALIES USING SPECTRAL ESTIMATION TECHNIQUES*

The depth of causative bodies may be derived from the power spectrum of their magnetic anomalies. After a short review of the theoretical basis of the method, several examples of its application to synthetic cases are shown. Disturbing effects due to improper choice of the sampling interval and to anomalies only partially contained in the examined segment of the profile are studied. The spectral method is then applied to real cases; the tapering of the anomalies, the application of non-linear filters and the effects of anomaly superpositions are investigated. In conclusion, an appraisal of the method and of its possible practical impact has been given.     

Coronal magnetic arcades with a helical structure

A family of force-free magnetic fields, describing multipolar magnetic arcades is described. The family allows for twisted magnetic ropes, where magnetic field lines have a helical shape. These helical structures are specific to filament channels of active regions. In their structure, arcades belong to normal-polarity fields and correspond to prominences of active regions according to available observational data. A potential field is found that is external with respect to an arcade. The constructed force-free fields can be used to further develop the MHD model of filaments that takes into account the gas pressure gradient and gravity.     

Driving Galactic Turbulence by Supernova Explosions

We investigate the general properties of supernova driven interstellar turbulence using local three-dimensional MHD simulations under Galactic conditions. Our model includes the effects of large-scale shear due to Galactic differential rotation, density stratification, compressibility, magnetic fields, heating via supernova explosions and parameterized radiative cooling of the interstellar medium. In addition to investigating isolated supernova explosions we allow for multiple supernovae distributed randomly in the Galactic disc and exponentially in the vertical direction. Single supernova explosions drive a strong shock, the lifetime of which is approximately 2 Myr in our model. This stage is found to be characterized by a kinetic energy spectrum in the diffuse gas with spectral index consistent with k = –2. Large-scale shear and Coriolis force act on the supernova remnant producing some vorticity inside it, but this process was found to be very weak. In the case of multiple supernova explosions, older remnants have an important role causing density fluctuations in the interstellar medium. In this clumpy medium, the propagation velocity of the shock fronts changes due to the changing density, and vorticity is generated. In the absence of these supernova interactions the kinetic energy spectrum shows a relatively wide shock spectrum with spectral index k = –2, but when the supernova interactions become dominant the classical k = –5/3 spectrum is observed.     

水平分层土层系统等效阻尼比的简化计算方法

在实际工程场地中,很多土层可视为水平分层,各层土的物理和力学性质存在差异,其中包括土的振动阻尼比。本文讨论水平分层土层系统的等效阻尼比的近似计算方法,基于5个不同的加权函数推导了10种等效阻尼比的计算公式。通过2个算例,分别以等效阻尼比为参数计算水平分层土层的地震反应,并与准确解相比较,分析了不同等效阻尼比近似计算方法的计算精度。数值结果表明,若等效阻尼比计算方法选择不恰当,会导致土层地震反应的计算结果出现较大误差。针对2种不同类型的水平分层土层,建议采用基于三角形分布的加权函数来计算土层系统的等效阻尼比。     

改进的阻尼法低纬度磁异常化极方法

中低纬度地区的磁异常资料进行化极处理能减少或消除斜磁化影响,对提高磁测资料的利用程度有重要意义.本文在总结分析原有的低纬度磁异常化极过程的基础上,改进了原有化极方法的滤波过程.改进后的阻尼化极方法,采用将振幅谱和相位谱分开进行滤波的过程,过程中只需要考虑对化极因子的振幅谱进行压制即可.通过公式分析表明该方法更加合理,减少了人为因素的影响;经单个模型和组合模型试验表明,此化极方法准确、稳定、可靠,可以有效地突出低纬度地区的异常体信息.     

三度体（物性随深度变化模型）位场波谱的正演计算

本文将均质的任意二维、三维物体位场的波谱解析表达式的研究成果推广到变密度、变磁化强度的更一般的情形。对密度差随深度呈指数函数衰减或线性变化的模型,获得了任意倾斜多边形质量面、斜平行六面体以及一般的多面体等形体的重力谱的解析表达式。它们的结构与均质体相应表达式一样简单,易于计算。以上结果表明,在很一般的条件下,位场波谱具有指数函数和的形式。     

A simplified method for unified buckling and free vibration analysis of pile-supported structures in seismically liquefiable soils

In seismic-prone zones with liquefiable deposit piles are routinely used to support structures (buildings/bridges). In this paper, a unified buckling and dynamic approach is taken to characterize this vibration. The pile–soil system is modelled as Euler–Bernoulli beam resting against an elastic support with axial load and a pile head mass with rotary inertia. The emphasis here is to obtain a simple expression that can be used by practicing engineers to obtain the fundamental frequency of the structure–pile–soil system. An approximate method based on an equivalent single-degree-of-freedom model has been proposed. Natural frequencies obtained from the exact analytical method are compared with approximate results. Proposed expressions are general as they are functions of non-dimensional parameters. It is shown that this simplified method captures the essential design features such as: (a) the continuous reduction of the first natural frequency of the structure–pile–soil system due to progressive reduction of soil stiffness due to liquefaction; (b) the reduction in the axial load-carrying capacity of the pile due to instability caused by liquefaction. The results derived in this paper have the potential to be directly applied in practice due to their simple yet general nature. An example problem has been taken to demonstrate the application of the method.     

ELECTROMAGNETIC FIELD OF SOURCES AT THE SURFACE OF A HOMOGENEOUS CONDUCTING HALFSPACE WITH HORIZONTAL ANISOTROPY: APPLICATION TO FISSURED MEDIA*

The electric and magnetic fields generated by horizontal electric and vertical magnetic dipoles lying on the surface of a conducting medium with horizontal anisotropy are investigated. Full expressions of their Fourier transforms are given, and the fields for a vertical magnetic dipole are calculated numerically. The radial and vertical magnetic components are found to be independent of the receiver-transmitter direction, whereas the other magnetic and electric components strongly vary with this direction. These results give useful criteria for defining the direction and amplitude of anisotropy from ground data; a ground experiment on fissured limestone was found to confirm the expected variations of the various field components. It is believed that this electromagnetic method can be used in order to provide information about the direction and amplitude of electric anisotropy.