A fast integral equation solver for 3D induction well logging in formations with large conductivity contrasts

Simulation of induction logging responses in formations with large conductivity contrasts is an important but challenging problem due to the singularity of a linear system caused by large contrasts. Also, three‐dimensional (3D) analysis of complex geophysical structures usually encounters high computational demands. In this paper, a pre‐corrected fast Fourier transform (pFFT)‐accelerated integral equation method is applied to overcome these difficulties. In the approach, the entire formation is included in the solution domain. The volume integral equation is set up in the region based on the fact that the total field is the summation of the excitation field and the secondary field. The emitted field by the transmitter coil (treated as a magnetic dipole) is regarded as the excitation of the system. Then the method of moments (MoM) is used to solve the integral equation. To reduce the high computational requirements of the MoM, the pFFT method is used to speed up the solution of the matrix equation and reduce the memory requirement as well. The resultant method is capable of computing induction logging problems involving large and complex formations. For problems with high conductivity contrasts, the solution of the matrix equation usually converges very slow or even fails to converge due to the large condition number of the coefficient matrix. To overcome this difficulty, an incomplete LU pre‐conditioner is used to significantly speed up the convergence of the matrix equation, thus further reducing the computation time. Numerical results show that the present method is efficient and flexible for 3D simulation of induction logging and is specifically superior for problems with high conductivity contrasts.     

The mathematical representation of freezing and thawing processes in variably-saturated,non-deformable soils

Recently, there has been a revival in the development of models simulating coupled heat and water transport in cold regions. These models represent significant advances in our ability to simulate the sensitivity of permafrost environments to future climate change. However, there are considerable differences in model formulations arising from the diverse backgrounds of researchers and practitioners in this field. The variability in existing model formulations warrants a review and synthesis of the underlying theory to demonstrate the implicit assumptions and limitations of a particular approach. This contribution examines various forms of the Clapeyron equation, the relationship between the soil moisture curve and soil freezing curve, and processes for developing soil freezing curves and hydraulic conductivity models for partially frozen soils. Where applicable, results from recent laboratory tests are presented to demonstrate the validity of existing theoretical formulations. Identified variations in model formulations form the basis for briefly comparing and contrasting existing models. Several unresolved questions are addressed to highlight the need for further research in this rapidly expanding field.     

Baseflow separation in a small watershed in New Brunswick,Canada, using a recursive digital filter calibrated with the conductivity mass balance method

Baseflow separation is important for obtaining critical parameters for hydrological models. As measuring the baseflow component directly is difficult, various analytical and empirical baseflow separation methods have been developed and tested. The recursive digital filter (RDF) method is commonly used for baseflow separation due to its simplicity and low data requirement. However, parameters used in the RDF method are often determined arbitrarily, resulting in high uncertainty of the estimated baseflow rate. A more accurate method is the conductivity mass balance (CMB) method, which is established based on the differences in physical processes between baseflow and surface runoff. In this research, the output of the CMB method was used to calibrate the parameters of an RDF model, and the calibrated RDF model was used to estimate monthly, seasonal and annual baseflow rate and baseflow index for the past 19 years using streamflow discharge records. The characteristics of the baseflow hydrographs were found to be consistent with the hydrological and hydrogeological conditions of the research area. Research results indicated that the accuracy of the RDF model has been greatly enhanced after being calibrated with the CMB method so that the RDF model can provide more reliable baseflow separation results for a long‐term study. Copyright © 2012 John Wiley & Sons, Ltd.     

西藏南部蛇绿岩套电导率研究

大地电磁(MT)资料显示,青藏高原地壳及地幔中普遍存在着高导层.作为大陆造山带中古洋盆岩石圈残片,蛇绿岩套的电导率测量可为了解古洋盆地区地壳及地幔的电性结构提供极其有用的信息.本研究中,我们在压力为1 GPa或3 GPa下,用交流阻抗谱法测量了采自西藏南部地区的蚀变辉长岩、玄武岩、角闪橄榄岩及方辉橄榄岩四个样品的阻抗谱,并进一步得出样品的电导率,不同样品电导率与温度之间的关系满足Arrhenius关系式.在实验温度范围内,蛇绿岩套电导率的对数logσ位于-6.0~-0.5 S/m之间,且随着温度的增高,不同样品电导率增大约4~5.5个量级.样品在未脱水的情况下,低温段的活化焓变化范围在0.4~0.6 eV之间,高温段的活化焓变化范围为1.7~2.6 eV之间.同时,我们研究了样品中结构水含量及铁含量对实验电导率的影响,验证了样品电导率与铁含量之间呈正比关系.当对样品结构水含量进行归一化后,相同温度下各样品的电导率随铁含量的增加而增大,而对样品铁含量归一化后,相同温度下各样品的电导率随样品中水含量的增加而增大.将实验电导率与藏南地区大地电磁结果进行了对比,发现本研究中各样品高温段实验电导率结果均落在大地电磁结果范围内.     

Regression‐Based Estimation of Lodgepole Pine Forest Age from Landsat Thematic Mapper Data

Abstract

A linear regression‐based model for mapping forest age using Landsat Thematic Mapper data is evaluated in the lodgepole pine forests of Yellowstone National Park. Regression models predicting age (R²=0.62) and a logarithmic transformation of age (R² = 0.90) used a combination of visible, near‐infrared, and middle‐infrared TM bands. Forest age maps produced using the regression method match broad‐scale patterns of forest age within the Yellowstone Central Plateau study area. Per‐pixel estimates of forest age, however, may depart substantially from actual forest age, particularly for older stands, and the maps are most appropriate for depicting regional patterns of forest age.     

中心回线式直升机TEM资料的电导率-深度成像方法

以理论上的中心回线式直升机航空瞬变电磁系统为例,提出一种快速的电导率-深度成像方法,将感生电动势瞬变响应数据变换为地下介质电导率-深度断面。该方法针对成像计算中视电导率非唯一性问题和成像深度难以确定的困难,首先给出一个关于均匀半空间模型电导率?的数据函数?(?)的变换公式,由此建立数据表???来查询视电导率值,可以得到可靠的视电导率值。其次用均匀半空间中的感应电场最大幅值对应深度值来定义视深度,建立视电导率-视深度的数据表,并在2个相邻时间道期间的最大电场深度变化的基础上定义成像深度,从而导出CDI结果。该方法的特点是理论简单、物理意义明确,且利用了数据表查询技术大大提高了计算速度,使得实时CDI处理成为可能。最后用理论模型进行试算,结果表明,电导率值的查询范围宽至10-4～102S/m时,高阻薄层、低阻薄层及二维目标体模型均有较好的成像效果。     

Thermal Properties of a Supraglacial Debris Layer with Respect to Lithology and Grain Size

This paper focuses on the impacts of debris cover on ice melt with regards to lithology and grain size. Ten test plots were established with different debris grain sizes and debris thicknesses consisting of different natural material. For each plot, values of thermal conductivity were determined. The observations revealed a clear dependence of the sub‐debris ice melt on the layer thickness, grain size, porosity and moisture content. For the sand fraction the moisture content played a dominant role. These test fields were water saturated most of the time, resulting in an increased thermal conductivity. Highly porous volcanic material protected the ice much more effectively from melting than similar layer thicknesses of the local mica schist. However, the analysis of thermal diffusivities demonstrated that the vertical moisture distribution of the debris cover must be taken into consideration, with the diffusivity values being significantly lower in deeper layers.     

坏囿矿化度地下水灌溉对民勤土壤环境的影响

通过野外试验和室内分析相结合的方法,研究了民勤绿洲农田不同矿化度(0.8g·L-1、2.0g·L-1和5.0g· L-1)地下水灌溉对土壤环境的影响.结果表明:随着灌溉水矿化度的增加,生长季水分消耗量逐渐降低,土壤总孔隙度和土壤有机质含量逐渐下降,而土壤有效水含量和土壤电导率则逐渐增加.当灌溉水矿化度为0.8g·L-1时,收获后土壤中上层(0～60 cm)电导率比播种前明显降低,而底层(60～90 cm)电导率略有增加,即淡水灌溉的淋溶作用明显;当灌溉水矿化度为2.0g,L-1时,生长季中的灌水依然对土壤盐分有淋溶作用,但明显弱于淡水的淋溶效果,收获后土壤盐分有明显的表聚现象;当灌溉水矿化度为5.0g·L-1时,收获后在0～60 cm深度电导率明显增加,在60～90 cm深度没有明显变化.     

Experimental hydrological forcing to illustrate water flow processes of a subarctic ladder fen peatland

Large peatland complexes dominate the landscape of the James Bay Lowland in subarctic Ontario, Canada. However, there is not a thorough understanding of the hydrological processes occurring in these important systems, particularly how ladder fens connect large domed bogs to the aquatic ecosystems that drain the peatland complex. Ladder fens consist of a pool‐rib topography where flow downgradient is controlled by the peat ribs. Within the ribs, low‐lying preferential flow paths typically enhance the transmission of water, whereas the elevated ridge microforms impede water flow to downgradient aquatic ecosystems. To assess the hydrological connectivity, we study the role of the water table, peat transmissivity, and microtopography of a small ladder fen for 3 summers (2013–2015) in the James Bay Lowland. The system was manipulated with a sustained hydrological forcing (water addition) to the upslope boundary of the fen during 2014 (38 m³/day) and 2015 (30 m³/day). There was an exponential increase in transmissivity towards the peat surface due to extremely high‐hydraulic conductivities within the upper few centimeters of the peat deposit. At the maximum water table, the saturated hydraulic conductivity of the 0.1 m layer of peat below the water table varied depending on peat microtopography (preferential flow paths = 42–598 m/day and ridges = 16–52 m/day), resulting in high‐hydrological connectivity periods. Furthermore, during 2015, there was an abnormally large amount of precipitation (300 mm vs. long‐term average ~ 100 mm) that resulted in complete surface water connectivity of the site. This caused rapid movement of water from the head of system to the outlet (~15 hr) and runoff ratios >1, compared to low‐water table periods (runoff ratio ~ 0.05). This study highlights the profound importance of the transmissivity–water table feedback mechanism in ladder fens, on controlling the water retention and drainage of large peatland complexes.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.