基于姿态变化的航空频率域电磁法仪器偏置实时校正方法研究

磁偶极子的航空频率域电磁法仪器在飞行测量的过程中由于仪器偏置的存在,且仪器偏置会随着外部气压、温度等环境因素以及收发线圈晃动的影响而呈现非线性变化,使得观测数据出现误差,因此需要对仪器偏置进行校正.而传统的在测线飞行前后将仪器抬至高空的“零场值”标定方法具有成本高、受测区环境影响大以及采用线性插值获取测线飞行过程中仪器偏置的精度低等缺点.本文根据仪器偏置与仪器姿态角变化无关的特性,通过测得仪器的姿态角信息,在满足重叠偶极子模型的条件下,实现对仪器偏置的高精度实时校正.模型仿真结果表明,在30 m常规飞行高度下,该方法实时测得的仪器偏置精度接近于110 m高空测得的精度;校正后仪器偏置的绝对误差与理论二次场的比值即相对误差小于5%,满足反演大地电导率的精度要求.该方法不仅减少了飞行的工作量,降低了飞行成本和飞行难度,而且可更加精确地获得测线飞行过程中仪器偏置的非线性变化值,提高航测数据的观测精度.     

利用2.5维数值模式匹配算法高效高精度建立垂直井眼中多分量阵列感应井眼校正库

本文将应用2.5维数值模式匹配算法研究建立柱状横向同性地层中偏心条件下磁流源并矢Green函数的高效算法,并利用该算法高效计算多分量阵列感应井眼校正库.首先,通过Fourier展开技术将偏心条件下磁流源并矢Green函数的数值模拟转化为一系列的轴对称问题;然后利用模式匹配算法求解轴对称问题,得到柱状介质中偏心条件下磁流源并矢Green函数的半解析解;最后,结合中国石油集团测井有限公司开发的三维感应测井仪器（TDIT）的结构参数,推导出复合线圈系视电导率张量的计算公式.在此基础上,根据井场实际地层资料归纳出理论地层模型相应参数的变化范围,高效高精度建立一套直接应用于工业生产的多分量阵列感应井眼校正库,并利用井眼校正库详细地分析考察各模型参数对TDIT响应的影响.

     

用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应

采用递推矩阵方法计算各向异性介质的磁流源并矢Green函数,并利用上述并矢Green函数对层状各向异性倾斜地层中多分量感应测井的响应进行数值模拟,分析了线圈距、层厚、倾角和围岩对多分量感应测井响应的影响.计算发现,共面视电导率比共轴视电导率变化规律复杂,不能反映地层电导率的真实情况.各向异性地层共轴视电导率随井眼相对倾角的增加而减小,而共面视电导率随井眼相对倾角的增加而增加.仪器的垂向分辨率、围岩各向异性对目的层响应的影响程度取决于线圈距,围岩各向异性对高电导率目的层中测井响应的影响大于对高电阻率目的层中测井响应的影响.     

Integral equation of relationship on the apparent conductivity and true conductivity in induction logging

This paper investigates the relation between the induced electromotive force measured by induction logging tool and the apparent conductivity, and the relation between the apparent conductivity and the formation true conductivity. Assuming the conductivity in Green’s function to be the function of the field point coordinate, the apparent conductivity expression of electric-field intensity is derived using Green’s formula, and the integral equation has been established representing the relationship of the apparent conductivity with the true conductivity under this condition. The integral equation is analyzed and then leads to the conclusion that the equivalent conductivity is the apparent conductivity and the values of the apparent conductivity function contain the true conductivity, and the method derived the true conductivity from the apparent conductivity around the well axis is put forward. The validity of the approach given in this paper is verified through numerical calculation. On the basis of above means, the transmitter coil produces an electric-field distribution in formation at every point where the induction logging tool moves along a wellbore, and a number of the receiving coils obtain an apparent conductivity distribution; this is what induction electric-field logging is.     

Inversion of multi-frequency electromagnetic induction data for 3D characterization of hydraulic conductivity

Electromagnetic induction (EMI) instruments provide rapid, noninvasive, and spatially dense data for characterization of soil and groundwater properties. Data from multi-frequency EMI tools can be inverted to provide quantitative electrical conductivity estimates as a function of depth. In this study, multi-frequency EMI data collected across an abandoned uranium mill site near Naturita, Colorado, USA, are inverted to produce vertical distribution of electrical conductivity (EC) across the site. The relation between measured apparent electrical conductivity (EC_a) and hydraulic conductivity (K) is weak (correlation coefficient of 0.20), whereas the correlation between the depth dependent EC obtained from the inversions, and K is sufficiently strong to be used for hydrologic estimation (correlation coefficient of ? 0.62). Depth-specific EC values were correlated with co-located K measurements to develop a site-specific ln(EC)–ln(K) relation. This petrophysical relation was applied to produce a spatially detailed map of K across the study area. A synthetic example based on EC_a values at the site was used to assess model resolution and correlation loss given variations in depth and/or measurement error. Results from synthetic modeling indicate that optimum correlation with K occurs at ~ 0.5 m followed by a gradual correlation loss of 90% at 2.3 m. These results are consistent with an analysis of depth of investigation (DOI) given the range of frequencies, transmitter–receiver separation, and measurement errors for the field data. DOIs were estimated at 2.0 ± 0.5 m depending on the soil conductivities. A 4-layer model, with varying thicknesses, was used to invert the EC_a to maximize available information within the aquifer region for improved correlations with K. Results show improved correlation between K and the corresponding inverted EC at similar depths, underscoring the importance of inversion in using multi-frequency EMI data for hydrologic estimation.     

Real time 2D inversion of induction logging data

Real time 2D inversion for an induction logging instrument may be achieved using a fast forward modeling and special inversion strategy. The fast forward modeling employs a low-frequency approximation of an induction response known as Doll's geometric factor. Modeling geometric factors is much faster than modeling the electromagnetic field in the frequency domain. To transform real data into the Doll's limit, multi-frequency skin-effect correction is applied. The correction technique involves an asymptotic theory of the integral equation for a 2D boundary value problem. The inversion is based on separating the parameter space into subspaces of lower dimension. Initially, adaptive overlapping windows split logging data into manageable portions. Each window consists of three subwindows: the predictor, corrector and upgrader. Further separation of parameters is introduced by Doll's approximation: the low-frequency response is linear with respect to formation conductivity. This allows us to split inversion for conductivity and geometric parameters. The next level of splitting inversion is achieved by independently determining parameters of the near borehole zone and remote formation areas. This is done by utilizing different subsets of sensors. The inversion does not require initial guess: layers are introduced dynamically, if necessary. The resolution is improved in sequential iterations by adding finer details to the previously obtained models. The final selection of parameters satisfies a variety of a priori constraints formulated as target resistivity distributions. The technique for imposing constraints is based on the analysis of data mapping into the model space. Interpretation of synthetic and real data confirms the viability of the method.     

兰州观象台FSQ和SSY-Ⅱ数字资料与模拟资料的对比

2002年对兰州观象台FSQ与SSY-Ⅱ两套仪器同时进行模拟与数字记录发现,模拟记录与数字记录有着较为明显的差异。为了对比模拟记录与数字记录的优劣,对FSQ的两套观测资料进行了年趋势变化、调和分析、日均值的相对噪声水平等的对比分析,对SSY-Ⅱ的资料进行了年趋势变化、调和分析、Nakai检验等的对比分析。通过对比发现模拟记录与数字记录的差异,是由于对这两种资料的处理不一致造成的。     

Behavior of electric potential fields over random conductivity,layered earth models

Solutions are derived for the potential distributions over one-layer and two-layer random conductivity earth models for the direct current resistivity method. Monopole and Wenner arrays are considered. It is shown that the random potential due to a spatially stationary Gaussian random conductivity is also Gaussian but not spatially stationary. The ensemble and sample statistics of the random kernel and apparent resistivity are examined. Representative curves of these functions are presented to show the effects of a random conductivity profile on them. It is shown that the effect of the random nature of the conductivity profile on an apparent resistivity sounding curve is greatest for small electrode spacings.     

The application of electromagnetic techniques to groundwater recharge investigations

The use of electrical methods for estimating spatial patterns of groundwater recharge was evaluated at a field site in southeastern Australia. Here, recharge increased from less than 0.2 mm year⁻¹ under native Eucalyptus vegetation, to between 1 and 14 mm year⁻¹ under dryland agriculture. This increase in recharge results in progressive leaching of salts in the soil profile. Differences in recharge can be estimated from differences in depth of leaching. The estimated recharge rates are correlated with soil texture, with higher recharge rates generally occurring through sandier soils. The relationships of recharge to salt content and soil texture both contribute to lower apparent electrical conductivities for higher recharge rates.

The effect of recharge rate on measured apparent electrical conductivities was modelled for various geophysical devices (including frequency-domain (FEM) and time-domain (TEM) electromagnetic instruments and direct current resistivity). The soil-texture effect was shown to have a greater effect than the solute leaching effect in determining the correlation between recharge and apparent electrical conductivity. Analysis of sensitivity to geological noise showed that variations in soil type below 2 m could disguise any correlation.

Correlations between recharge rate, measured at core sites from chloride tracer techniques, and apparent electrical conductivity, measured with FEM electromagnetic devices, supported the conclusions of the model. For DC resistivity and TEM methods, correlations between recharge and apparent electrical conductivity were not significant, although for resistivity this may be due partly to the small number of measurements made. The FEM device most sensitive to variations in recharge had an operating frequency of 9.8 kHz. At lower frequencies the sensitivity is reduced, as the instruments are sensing too deeply. The poor correlations for TEM, as compared with FEM, are due probably to the relatively deeper penetration of the TEM instrument used in the study, rather than any inherent differences between the techniques.

Because the major reason for the correlation between recharge and apparent electrical conductivity is soil texture, in this area the geophysical devices are mostly mapping soil type.     

FHD仪偏置线圈电流反向对D分量基线值的影响及修正方法

在地磁连续观测中,需要通过分析计算观测数据对仪器的工作状态作出正确的判断,介绍了判断FHD仪偏置线圈电流是否反向的2种方法,研究了反向前后基线值的变化情况并提出了修正方法。     

一体化低功耗宽频带数字地震仪研制

流动观测台网与固定观测台网的结合是当前地震观测技术系统发展的一个趋势. 针对流动观测时地震计与地震数据采集器互相独立、 携带不方便、 功耗高的问题, 自主研发了一款适合流动观测的集地震信号提取、 数据采集、 记录和服务为一体的数字地震仪. 该仪器具有频带宽(60 s—80 Hz)、 动态范围高(>140 dB)、 功耗低(0.6 W)、 携带方便(整机重量在15 kg左右, 包括供电系统、 GPS天线和包装箱)等特点. 详细介绍了该地震仪的外观结构、 整体硬件结构、 低功耗处理技术和所采用的灵敏度校正、 标准方位和正交校正技术. 对仪器的主要性能参数指标进行了严格的测试, 并给出了具体的测试结果. 该仪器研制完成后, 投入到了青海玉树M_S7.1地震震后流动观测中. 从半年的实际使用结果来看, 该仪器能够满足流动观测的要求.     

坑(井)-地多参数电磁接收系统

坑（井）-地多参数电磁接收系统在地面、坑道、井中三维空间观测天然场源及人工源电磁信号,观测装置接近或穿过矿体,借助动源、大功率发射可获得更加明显的异常值,旨在加大探测深度的基础上提高深部分辨率,同时获取电阻率、极化率、复阻抗等参数.为坑（井）-地电磁成像方法现场数据采集提供硬件设备支持,针对坑道、井中电磁观测的特殊需求,开发了由不极化电极、非接触电极、小型感应式线圈、微型三分量音频磁传感器、地面-坑道电磁接收机、井中电磁接收机等组成的坑（井）-地多参数电磁接收系统,解决了传统不极化电极坑道硬岩表面接地困难和传统感应式磁传感器难以适应坑道狭窄空间的问题;通过多层次的室内、弱干扰条件下野外测试验证了接收系统的各项技术指标,接收机通道全频段噪声水平接近10 nV/rt（Hz）,微型三分量音频磁传感器体积压缩至32 cm×32 cm×32 cm;强干扰条件下的矿山试验验证了接收系统的可靠性、适用性及先进性.

     

强震动观测仪器面临的机遇和挑战

本文介绍了国内外强震仪、力平衡加速度计、烈度计、MEMS加速度仪等强震动观测仪器的发展历史与现状,概述了强震动记录中典型异常波形产生的原因和力平衡加速度计的仪器响应误差及校正,建议了解决强震仪器缺陷的措施,分析了强震仪在强震动观测发展中面临的机遇和挑战。研究结果如下:① 进行强震动记录异常波形的研究是有针对性改进现有强震仪存在问题的很好途径;② 早期数字强震仪（力平衡加速度计频带范围0至30Hz）获取的强震动记录应进行仪器校正;③ 应跟踪强震动观测新方法和相关领域新技术,发展MEMS加速度仪和光纤加速度仪等新型仪器。     

Relationships between magnetic susceptibility and heavy metals in urban topsoils in the arid region of Isfahan, central Iran

Electromagnetic induction (EMI) method results are shown for vertical magnetic dipole (VMD) configuration by using the EM38 equipment. Performance in the location of metallic pipes and electrical cables is compared as a function of instrumental drift correction by linear and quadratic adjusting under controlled conditions. Metallic pipes and electrical cables are buried at the IAG/USP shallow geophysical test site in São Paulo City, Brazil. Results show that apparent electrical conductivity and magnetic susceptibility data were affected by ambient temperature variation. In order to obtain better contrast between background and metallic targets it was necessary to correct the drift. This correction was accomplished by using linear and quadratic relation between conductivity/susceptibility and temperature intending comparative studies. The correction of temperature drift by using a quadratic relation was effective, showing that all metallic targets were located as well deeper targets were also improved.     

Calibration shifts in a LaCoste-and-Romberg gravimeter: comparison with a Scintrex CG-3M

The reliability and ruggedness of LaCoste-and-Romberg gravity meters has meant that their use has extended over recent years from gravity mapping to the investigation of small gravity changes (> 10 μGal) with time. When using the instrument at the limit of its capability, limitations become apparent. Some instruments drift more than others and some are prone to excessive tares. However, with regular servicing these problems are minimized. Here we demonstrate clear evidence for calibration changes of the order of 1000 ppm for LaCoste-and-Romberg meter G513, which, for gravity differences of over 50 mGal, represents a highly significant induced error of over 50 μGal. Data collected over 15 years using nearly 4000 mGal of the 7000 mGal range of the instrument illustrates that long periods of stability may be interspersed with periods lasting several months while the calibration changes by up to 0.1%. This is comparable with the calibration drift of a new Scintrex CG-3M instrument, although the process causing the change may not be the same.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Fhe application of electromagnetic-induction on the measurement of sea ice thickness in the Antarctic

As an important component of the cryosphere, sea ice is very sensitive to climate change. The study of sea ice physics needs accurate sea ice thickness. This paper presents an electromagnetic induction （EM） technique which can be used to measure the sea ice thickness distribution efficiently and its successful application in the Antarctic Neila Fjord. Based on the electrical properties of sea ice and seawater and the application of electromagnetic field theory, this technique can accurately detect the distance between the EM instrument and the ice/water interface to measure the sea ice thickness. Analyzing the apparent conductivity data obtained by the electromagnetic induction technique and drill-hole measurements at same location allows the construction of a transform equation for the apparent conductivity and sea ice thickness. The verification of the calculated sea ice thickness using this equation indicates that the electromagnetic induction technique is able to determine reliable sea ice thickness with an average relative error of only 5.5%. The ice thickness profiles show the sea ice distribution in Neila Fjord is basically level with a thickness of 0.8 - 1.4 m.     

GSM-19T质子旋进磁力仪比测实验分析

为了检测GSM-19T质子旋进磁力仪主要技术指标及参数的性能,2019年10月17—21日河北省地震局红山基准台以GSM-19W Overhauser磁力仪为标准仪器,对40台GSM-19T质子旋进磁力仪地磁总强度F,通过同步连续采样观测和仪器差进行了比测,计算其差值和均方差。结果表明,大部分测量仪器与标准仪器具有较好的一致性和稳定性,地磁总强度F均方差和仪器差均较小,各项技术指标参数符合大纲测试要求,测量仪器比测合格,可为基站观测和野外流动观测提供准确可靠的数据和仪器保障;小部分测量仪器有漂移,地磁总强度F均方差稍大。     

Distributed snow accumulation and ablation modeling in the American River basin

A spatially distributed snow model procedure for estimating snow melt, snow water equivalent and snow cover area is formulated and tested with data from the American River basin in California’s Sierra Nevada. An adaptation of the operational National Weather Service snow accumulation and ablation model is used for each model grid cell forced by spatially distributed precipitation and temperature data. The model was implemented with 6-hourly time steps on 1 km² grid cells for the snow season of 1999–2003. Temperature is spatially interpolated using the prevailing lapse rate and digital terrain elevation data. Precipitation is spatially interpolated using regional climatological analyses obtained from PRISM. Parameters that control snow melt are distributed using ground surface aspect. The model simulations are compared with data from 12 snow-sensors located in the basin and the daily 500-m snow cover extent product from the MODIS/Terra satellite mission. The results show that the distribution of snow pack over the area is generally captured. The snow pack quantity compared to snow gauges is well estimated in high elevations with increasing uncertainty in the snow pack at lower elevations. Sensitivity and uncertainty analyses indicate that the significant input uncertainty for precipitation and temperature is primarily responsible for model errors in lower elevations and near the snow line. The model is suitable for producing spatially resolved realistic snow pack simulations when forced with operationally available observed or predicted data.     

The dual-beam mini-DOAS technique—measurements of volcanic gas emission,plume height and plume speed with a single instrument

The largest error in determining volcanic gas fluxes using ground based optical remote sensing instruments is typically the determination of the plume speed, and in the case of fixed scanning instruments also the plume height. We here present a newly developed technique capable of measuring plume height, plume speed and gas flux using one single instrument by simultaneously collecting scattered sunlight in two directions. The angle between the two measurement directions is fixed, removing the need for time consuming in-field calibrations. The plume height and gas flux is measured by traversing the plume and the plume speed is measured by performing a stationary measurement underneath the plume. The instrument was tested in a field campaign in May 2005 at Mt. Etna, Italy, where the measured results are compared to wind fields derived from a meso-scale meteorological model (MM5). The test and comparison show that the instrument is functioning and capable of estimating wind speed at the plume height.