共查询到20条相似文献，搜索用时 46 毫秒
1.
地磁测深研究的周期范围通常为10~5~10~7s,缺少反映浅部电性结构的短周期信息,而C-响应受浅部电阻率影响明显,因此本文提出在反演中增加浅部(约200 km)电阻率约束以提高深部反演的稳定性和可靠性.在磁层环状电流满足P_1~0假设的条件下,球坐标系中一维导电薄球层状地球的C-响应和电导率分布关系由边界条件通过递推的方法计算得到.反演采用有限内存拟牛顿(L-BFGS)法;浅部电阻率约束通过将目标函数对模型参数的梯度设为零来实现;通过置信区间分析评价约束反演结果的可靠性.合成数据的无约束反演虽然最终的拟合效果很好,但浅部电阻率受初始模型影响,差异较大;采用浅部约束后,反演结果对初始模型依赖性明显减小,同时还能显著提高200~600 km范围内反演结果的准确性.对全球近地轨道卫星观测的C-响应数据约束反演后结果与前人一致,表现为地幔电导率整体上随着深度的增加而增加.参数置信区间分析表明,由于约束反演加入了浅部信息,电阻率的变化范围更加紧致,说明反演结果更加可靠.因此,有必要通过其他地球物理方法,如长周期大地电磁测深等获得浅部电阻率分布,作为先验信息参加反演,进行浅部约束的C-响应反演,获得更可靠的一维全深度电性结构,为地磁测深数据解释奠定基础. 相似文献
2.
3.
全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构. 相似文献
4.
W. D. Parkinson 《Pure and Applied Geophysics》1987,125(2-3):459-463
Diurnal variations supply only four frequencies with sufficient power to be useful, giving an overall frequency ratio of only 4. Also only permanent observatories supply sufficient accurate data, and their distribution on the Earth is very irregular. More importantly, the assumption implicit in deep conductivity determinations, that the conductivity is a function of depth only, breaks down because of near surface conductivity anomalies. 相似文献
5.
大地电磁测深理论与数据处理解释均假定平面电磁波垂直入射地下空间,但随着研究尺度的逐渐扩大,使得因地球弧度产生的影响难以忽略.此时,传统笛卡尔坐标体系及平面波场源不再适用于大尺度的大地电磁数据正反演解释.本文提出并实现了一种基于球坐标系的大地电磁交错网格有限差分三维正演算法,并对电场进行极向-环向分解,结合球谐函数和贝塞尔函数构建了可取代平面波的场源模型.首先利用经纬度信息构建三维地电模型,将场源设置于模型空间正上方,然后通过直接求解球坐标系下麦克斯韦方程来获得大地电磁响应.在此基础上,本文设计了球坐标下具有不同分辨率的多个三维地电模型,阐述了由球体模型到笛卡尔模型的转换方法,详细对比了两种坐标体系在计算效率、所求得的电场和视电阻率方面的差异.结果表明二者差异度主要与电性横向分布和地图投影方法有关,与周期并不存在明显的单调递增关系. 相似文献
6.
V. V. Spichak 《Izvestiya Physics of the Solid Earth》2011,47(1):72-76
Possibilities for three-dimensional (3D) magnetotelluric (MT) sounding of local objects contained in the Earth’s crust are estimated in a case study of the magma chamber of the Vesuvius volcano. Stochastic inversion of the model MT data by the Markov Chain Monte Carlo (MCMC) method has shown that the most efficient approach is not simultaneous but successive estimation of the geometry and the depth of the anomaly and the assessment of the conductivity distribution within the anomalous region. A zone of equivalence is revealed between the a priori estimate of the depth of the anomalous zone and the a posteriori distribution of electric conductivity within it. Based on the present estimation and previous results, an algorithm for determination of the parameters of local crustal anomaly is proposed. 相似文献
7.
A program for identifying magnetic clouds in patrol satellite data, which recorded the interplanetary medium parameters near
the magnetosphere, has been developed based on the cloud model in the form of a force-free cylindrical flux tube. The program
makes it possible to also determine the entire magnetic field distribution in a cloud that approaches the Earth, using the
initial satellite measurements. For this purpose, a model cloud (which has the maximal correlation coefficient with an analyzed
cloud with respect to three magnetic field vector components and minimal rms deviations of the magnetic field and velocity
components) is selected from the preliminarily created database including 2 million model clouds. The obtained magnetic field
distribution in a cloud will make it possible to predict the intensity of a magnetic storm that this cloud will cause. 相似文献
8.
This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data. 相似文献
9.
A procedure is suggested of a more effective and faster computation of the impedance, the transfer function and amplitudes of the induced field in a spherically symmetric model of the electrical conductivity. The existing induction data have been supplemented by about 80 new values derived from the analysis of daily means. The fit of the existing 1-D models of the electrical conductivity of the mantle to the set of induction data is investigated. The characteristic equation for the free electromagnetic oscillations of a radially inhomogeneous Earth is derived and its possible importance in solving the inverse problem of electric conductivity is pointed out. 相似文献
10.
The spatial distributions of electric fields and currents in the Earth’s atmosphere are calculated. Electric potential distributions
typical of substorms and quiet geomagnetic conditions are specified in the ionosphere. The Earth is treated as a perfect conductor.
The atmosphere is considered as a spherical layer with a given height dependence of electrical conductivity. With the chosen
conductivity model and an ionospheric potential of 300 kV with respect to the Earth, the electric field near the ground is
vertical and reaches 110 Vm−1. With the 60-kV potential difference in the polar cap of the ionosphere, the electric field disturbances with a vertical
component of up to 13 V m−1 can occur in the atmosphere. These disturbances are maximal near the ground. If the horizontal scales of field nonuniformity
are over 100 km, the vertical component of the electric field near the ground can be calculated with the one-dimensional model.
The field and current distributions in the upper atmosphere can be obtained only from the three-dimensional model. The numerical
method for solving electrical conductivity problems makes it possible to take into account conductivity inhomogeneities and
the ground relief. 相似文献
11.
传统的均匀分布异常体模型,不能准确描述地下介质的不规则性变化,采用随机介质进行替代,将使电导率这一物性参数更加接近于实际的电导率分布.本文结合地质结构中断裂带特征,采用Von Kármán函数建立三维随机介质模型,通过讨论Hurst指数与自相关长度对三维随机电导率模型建模产生的影响,进行参数优化建立所需的随机介质模型.利用一维傅里叶变换建立三维随机变化的条状断裂带异常体,三维傅里叶变换建立三维随机变化的背景围岩,准确地表征了油气藏断裂带特征.基于时域有限差分方法,实现了磁源激励下的三维随机介质航空时域电磁响应数值模拟.采用均匀半空间模型验证了数值模拟的正确性,分析了随机断裂带与均匀断裂带的电磁响应特征,结果表明随机断裂带可以准确描述地下介质的分布特征,而且与断裂带垂直方向的电磁响应特征清晰地描述了断裂带的倾向、走向与位置,为断裂带结构探测提供了理论依据和技术指导,三维随机断裂带模拟方法同样适用于其它三维随机介质的数值模拟. 相似文献
12.
传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用. 相似文献
13.
A spherical approximation makes the basis for a majority of formulas in physical geodesy. However, the present-day accuracy
in determining the disturbing potential requires an ellipsoidal approximation. The paper deals with constructing Green’s function
for an ellipsoidal Earth by an ellipsoidal harmonic expansion and using it for determining the disturbing potential. From
the result obtained the part that corresponds to the spherical approximation has been extracted. Green’s function is known
to depend just on the geometry of the surface where boundary values are given. Thus, it can be calculated irrespective of
the gravity data completeness. No changes of gravity data have an effect on Green’s function and they can be easily taken
into account if the function has already been constructed. Such a method, therefore, can be useful in determining the disturbing
potential of an ellipsoidal Earth. 相似文献
14.
XinZhuan Guo 《中国科学:地球科学(英文版)》2016,59(4):696-706
Hydrous minerals are important water carriers in the crust and the mantle, especially in the subduction zone. With the recent development of the experimental technique, studies of the electrical conductivity of hydrous silicate minerals under controlled temperature, pressure and oxygen fugacity, have helped to constrain the water distribution in the Earth’s interior. This paper introduces high pressure and temperature experimental study of electrical conductivity measurement of hydrous minerals such as serpentine, talc, brucite, phase A, super hydrous phase B and phase D, and assesses the data quality of the above minerals. The dehydration effect and the pressure effect on the bulk conductivity of the hydrous minerals are specifically emphasized. The conduction mechanism of hydrous minerals and the electrical structure of the subduction zone are discussed based on the available conductivity data. Finally, the potential research fields of the electrical conductivity of hydrous minerals is presented. 相似文献
15.
A Douglas-Jones Predictor-Corrector Program for Simulating One-Dimensional Unsaturated Flow in Soil 总被引：1，自引：0，他引：1
C. Babajimopoulos 《Ground water》1991,29(2):267-270
Abstract. A fully documented program to represent one-dimensional unsaturated flow in soil is described. The program is based on a Douglas-Jones finite-difference implicit method to solve the Richards equation. An implicit linearization scheme is used to estimate the hydraulic conductivity and specific moisture capacity functions. Predicted moisture content profiles compared with two Galerkin finite-element solutions and field observations on a Panoche clay loam soil show very good agreement. 相似文献
16.
The current functions calculated for a spherical Earth model with real near-surface conductivity distribution are presented. The system is excited by a uniform external field of three different polarizations. 相似文献
17.
A two‐and‐half dimensional model‐based inversion algorithm for the reconstruction of geometry and conductivity of unknown regions using marine controlled‐source electromagnetic (CSEM) data is presented. In the model‐based inversion, the inversion domain is described by the so‐called regional conductivity model and both geometry and material parameters associated with this model are reconstructed in the inversion process. This method has the advantage of using a priori information such as the background conductivity distribution, structural information extracted from seismic and/or gravity measurements, and/or inversion results a priori derived from a pixel‐based inversion method. By incorporating this a priori information, the number of unknown parameters to be retrieved becomes significantly reduced. The inversion method is the regularized Gauss‐Newton minimization scheme. The robustness of the inversion is enhanced by adopting nonlinear constraints and applying a quadratic line search algorithm to the optimization process. We also introduce the adjoint formulation to calculate the Jacobian matrix with respect to the geometrical parameters. The model‐based inversion method is validated by using several numerical examples including the inversion of the Troll field data. These results show that the model‐based inversion method can quantitatively reconstruct the shapes and conductivities of reservoirs. 相似文献
18.
The conductivity anisotropy behaviour is described for certain environment in the Earth’s crust and the MT inversion method
for a layered symmetrically anisotropic model is presented. The inversion interpretations of the anisotropic model from the
observational data are helpful to identify the earthquake precusors as indicated by the deep conductivity anisotropic variations,
and also provide some useful information to investigate the stress states and deformation bands in the deep crust of the Earth.
This project is sponsored by the Chinese Joint Seismological Science Foundation. 相似文献
19.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献
20.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献