首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 171 毫秒
1.
视复电阻率频谱的一种近似反演方法   总被引:5,自引:0,他引:5  
基于稀释系数不随频率变化的假设,导出了复电阻率视频谱与真频谱及柯尔-柯尔频散视参数与真参数之间的近似关系式。通过对理论视频谱进行的反演和计算证明,即使在稀释系数幅值随频率的相对变化达100%的情况下,选用特征频率的稀释系数值,所导出的视参数与真参数的近似关系式仍有实用价值。从而为频谱激电法资料的反演,计算极化体真参数提供了一种简便而有效的近似方法。  相似文献   

2.
一种反演频谱激电法视频谱求取真参数的方法   总被引:4,自引:0,他引:4  
在激电效应引起的电阻率变化范围内,将稀释系数假设为常数的条件下,提出了利用归一化视谱和归一化真谱之间近似关系式,直接反演视谱求极化体真参数的方法。大量的反演算例表明,若固定m_2等于极化体的真充电率值,这种反演方法求得的参数β~'_2、τ~'_2和C~'_2值的稳定性很好。反演求得之c~'_2即为极化体的固有频率相关系数c_2:β~'_2约等于稀释系数实分量频谱的拐点值R_cβ_2;τ_2值虽然保持稳定,但和极化体的真时间常数值不同,不过,引入适当的系数可估算出真时间常数τ_2。  相似文献   

3.
本文根据数值模拟结果讨论了同时存在激电和电磁效应时,频谱激电法的异常性质.指出了低阻和高阻极化体上激电和电磁响应的空间分布形态和频率特性.验证了按两个柯尔- 柯尔模型的组合反演视复电阻率频谱,分离激电和电磁效应的合理性.并利用这种方法反演理论频谱的结果,讨论了激电视频谱参数的异常性态.文中还讨论了覆盖层和非各向同性对频谱激电法异常的影响.  相似文献   

4.
这里研究层状极化介质激发极化效应特征。首先,在忽略电磁效应的情况下,用柯尔~柯尔模型表征极化介质复电阻率,计算三层不同地电断面复电阻率的振幅和相位随极距变化的曲线;然后,分析极化层不同充电率和不同时间常数、不同激发频率和极化层不同埋深,对测深曲线的影响规律,并提出了具有实际指导意义的极化层埋深与极距的对应关系。这些工作都能够为野外开展复电阻率测量和资料解释起到参考作用。  相似文献   

5.
王书民 《物探与化探》2004,28(4):317-319
利用复电阻率柯尔—柯尔模型得出的激电相位与极化率的关系表达式,计算并绘制了相位随极化率的变化曲线,通过对激电相位随极化率的变化曲线特征分析,总结出了其近似关系,提出了激电相位近似的解释方法,并取得了实际应用效果。  相似文献   

6.
大地电磁测深实际测量中,视电阻率的静位移是常见现象,有时会遇到连片静位移的区域,能否对其正确识别与校正将直接影响到反演解释结果的准确度。针对出现连续静位移区的大地电磁测深剖面,相邻测点间同极化模式的视电阻率曲线可参考性较低。本文基于岩性分析,通过实验,初步总结出一套针对连片区的静态位移校正方法——首枝频点统计求平均法。本方法先对研究区浅地表同一地层分布区的畸变测点进行统计,再统计每个测点(相同极化模式)前三个频点的平均值,将这些畸变点的平均值再进一步做平均,得到的值将会作为相应畸变点视电阻率曲线首枝(第二个频点)的预置位置。该方法不改变测点的曲线形态,只改变测点曲线首枝位置,使其达到统计出的预置位置。将该方法应用在扎鲁特地区大地电磁测深实测资料的静校正实验中,取得了较好的应用效果。这也说明,大地电磁测深的剖面连片静位移区资料处理,在没有大量的测井和物性约束条件下,使用首枝频点统计求平均法具有一定的合理性和参考性。  相似文献   

7.
在频率未达相位峰值频率时,不极化水平均匀大地和二层大地的感应耦合谱的每个频段都可近似地用一个Cole-Cole谱拟合。极化水平均匀大地和底层极化的二层大地的总视复电阻率谱的每个频段都可用乘法或加法组合的两个Cole-Cole谱很好拟合,反演得到的激电谱都近似等于激电效应的贡献,但只有乘法组合反演得到的感应耦合才能代表感应耦合的贡献。  相似文献   

8.
在前人工作基础上,对电测深法一维正演解析解递推电阻率转换函数求偏导,给出了电测深法视电阻率对介质电阻率偏导数的解析计算方法;根据Seigel体激发极化理论,采用视电阻率对介质电阻率偏导数矩阵的线性组合,计算出了电测深法的视极化率,其结果与"等效电阻率法"计算视极化率最大相对误差小于1.1%。在一次正演计算过程中就可以同时得到视电阻率、视极化率以及视电阻率对介质电阻率的偏导数矩阵,提高了计算效率,也为后续的多层地电模型电测深法数据优化反演计算提供了关键性技术。  相似文献   

9.
极化大地复电阻率正反演   总被引:3,自引:0,他引:3  
论述了层状极化大地复电阻率正、反演算法。给出了层状极化大地条件下三极和对称四极装置视复电阻率的解析表达式,并引入线性滤波方法实现表达式的求解。正演算例计算结果表明;随着频率的降低,复电阻率相位峰值增大;随着极化层层厚增加,复电阻率虚部平缓变化后连续递增。介绍了递推反演算法的实现过程,从理论模型试算结果看,测深曲线较好地反映了极化层Cole-Cole参数随深度的变化规律,其曲线的首支和尾支趋近于极化层的真Cole-Cole参数,表明了反演计算的可行性。  相似文献   

10.
本文研究了用有限元数值方法模拟柯尔—柯尔复电阻率模型异常规律的问题。计算例子的精度满足实际要求。用数值方法计算了位于极化半空间中的直立板和位于不极化半空间中的水平板的复电阻率剖面响应曲线和相位谱。数值结果表明:在复电阻率振幅测量中,若围岩和矿体的电性差异不大,则工作频率的选择十分重要,否则可能观测不到异常。  相似文献   

11.
由视谱直接反演真谱参数的实用方法   总被引:2,自引:0,他引:2  
张桂青  罗延钟 《地球科学》1991,16(2):229-239
  相似文献   

12.
A numerical experiment has been applied to explore the potentialities and the limitations of the inversion of IP-affected TEM responses of a two-layer earth with a polarizable layer above (Model 1) or below (Model 2) a nonpolarizable layer. The IP effects have been incorporated into the models via a Cole–Cole complex frequency-dependent conductivity. One of us computed synthetic in-loop and coincident-loop transient responses with added Gaussian noise, and the other performed single and joint inversion of the two sets of pseudo-experimental data. Model 1 turns out to be advantageous over Model 2 in TEM applications and gives a good fit of the Cole–Cole parameters (chargeability, IP time constant, and exponent) even in the absence of a priori information. In the case of Model 2, the lack of a priori information causes problems as to recognize which layer is polarizable, and the fit of the Cole–Cole parameters is generally worse. The layer thicknesses and resistivities are rather accurate in both groups of models, irrespective of whether a priori information is available. As the upper layer increases in thickness (H1), the fit of its parameters ever improves in both models while the parameters of the lower layer, on the contrary, contain a greater error. Joint inversion of in-loop and coincident-loop transients improves the fit in most cases.Relative rms error (σrel) does not depend on the upper layer thickness for Model 1 but decreases as H1 increases in the case of Model 2. The error in joint inversion is times that in single inversion, which means that additional criteria other than σrel may be useful to estimate the inversion quality.  相似文献   

13.
将Cole—Cole模型表示的地下异常体各剖分单元复电阻率,引入到体积分方程求解均匀大地三维电磁场计算中,实现电偶源地面激发、接收复电阻率三维电磁场正演模拟。  相似文献   

14.
二维视电阻率断面的快速最小二乘反演   总被引:12,自引:0,他引:12  
快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方  相似文献   

15.
The non-inductive galvanic disturbances due to surficial bodies, lying smaller than high frequency skin depth, cause serious interpretational errors in magnetotelluric data. These frequency independent distortions result in a quasi-static shift between the apparent resistivity curves known as static shift. Two-dimensional modelling studies, for the effects of surficial bodies on magnetotelluric interpretation, show that the transverse electric (TE) mode apparent resistivity curves are hardly affected compared to the transverse magnetic (TM) mode curves, facilitating the correction by using a curve shifting method to match low frequency asymptotes. But in the case of field data the problem is rather complicated because of the random distribution of geometry and conductivity of near surface inhomogeneities. Here we present the use of deep resistivity sounding (DRS) data to constrain MT static shift. Direct current sensitivity studies show that the behaviour of MT static shift can be estimated using DC resistivity measurements close to the MT sounding station to appreciable depths. The distorted data set is corrected using the MT response for DRS model and further subject to joint inversion with DRS data. Joint inversion leads to better estimation of MT parameters compared to the separate inversion of data sets.  相似文献   

16.
油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。  相似文献   

17.
A numerical experiment was applied to explore the potentialities and limitations of joint inversion of IP-affected transients measured with different loop configurations above a uniform half-space with a Cole-Cole complex conductivity. One of us calculated 200 m × 50 m and 50 m × 50 m loop responses of a uniform polarizable conductor with varied Cole-Cole parameters and imposed synthetic Gaussian noise that simulated measurement errors. Then the generated pseudo-experimental data passed to the other co-author who performed single and joint 1D inversion twice: first being unaware of the “true” underlying models and then after being told that they all were represented by a uniform polarizable earth. More than a half of the fitted models provided a good idea of the true models though misfit was quite large in some cases. The fit was better in single inversion with a priori information available, and improved further through joint inversion of central-loop and coincident-loop responses. Joint inversion with a priori information known was of good quality even at a chargeability as low as 0.02. The standard error in joint inversion was times the measurement error and depended mainly on fitting errors for smaller-loop data. The reason is that the smaller-loop transients included a non-monotonous interval where the signal changed rapidly under the effect of fast-decaying induced polarization.  相似文献   

18.
An evolutionary inversion approach is suggested for the interpretation of nuclear and resistivity logs measured by direct-push tools in shallow unsaturated sediments. The efficiency of formation evaluation is improved by estimating simultaneously (1) the petrophysical properties that vary rapidly along a drill hole with depth and (2) the zone parameters that can be treated as constant, in one inversion procedure. In the workflow, the fractional volumes of water, air, matrix and clay are estimated in adjacent depths by linearized inversion, whereas the clay and matrix properties are updated using a float-encoded genetic meta-algorithm. The proposed inversion method provides an objective estimate of the zone parameters that appear in the tool response equations applied to solve the forward problem, which can significantly increase the reliability of the petrophysical model as opposed to setting these parameters arbitrarily. The global optimization meta-algorithm not only assures the best fit between the measured and calculated data but also gives a reliable solution, practically independent of the initial model, as laboratory data are unnecessary in the inversion procedure. The feasibility test uses engineering geophysical sounding logs observed in an unsaturated loessy-sandy formation in Hungary. The multi-borehole extension of the inversion technique is developed to determine the petrophysical properties and their estimation errors along a profile of drill holes. The genetic meta-algorithmic inversion method is recommended for hydrogeophysical logging applications of various kinds to automatically extract the volumetric ratios of rock and fluid constituents as well as the most important zone parameters in a reliable inversion procedure.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号