Two-dimensional physical-based inversion of confined and unconfined aquifers under unknown boundary conditions

An inverse method is developed to simultaneously estimate multiple hydraulic conductivities, source/sink strengths, and boundary conditions, for two-dimensional confined and unconfined aquifers under non-pumping or pumping conditions. The method incorporates noisy observed data (hydraulic heads, groundwater fluxes, or well rates) at measurement locations. With a set of hybrid formulations, given sufficient measurement data, the method yields well-posed systems of equations that can be solved efficiently via nonlinear optimization. The solution is stable when measurement errors are increased. The method is successfully tested on problems with regular and irregular geometries, different heterogeneity patterns and variances (maximum K_max/K_min tested is 10,000), and error magnitudes. Under non-pumping conditions, when error-free observed data are used, the estimated conductivities and recharge rates are accurate within 8% of the true values. When data contain increasing errors, the estimated parameters become less accurate, as expected. For problems where the underlying parameter variation is unknown, equivalent conductivities and average recharge rates can be estimated. Under pumping (and/or injection) conditions, a hybrid formulation is developed to address these local source/sink effects, while different types of boundary conditions can also exert significant influences on drawdowns. Local grid refinement near wells is not needed to obtain accurate results, thus inversion is successful with coarse inverse grids, leading to high computation efficiency. Furthermore, flux measurements are not needed for the inversion to succeed; data requirement of the method is thus not much different from that of interpreting classic well tests. Finally, inversion accuracy is not sensitive to the degree of nonlinearity of the flow equations. Performance of the inverse method for confined and unconfined aquifer problems is similar in terms of the accuracy of the estimated parameters, the recovered head fields, and the solver speed.     

Conductivity structure and rheological property of lithosphere in Southern Tibet inferred from super-broadband magnetotelluric sounding

To understand deep lithosphere structure beneath the Qinghai-Tibet Plateau more comprehensively and objectively and to explore important scientific issues,such as characteristics of plateau lithospheric deformation,state of strain,thermal structure,plate (or terrane) movement,and crust-mantle rheology,it is necessary to research the variation of crust-mantle electrical structure in the east-west direction in every geological unit.For this purpose,six super-broadband magnetotelluric (MT) sounding profiles ha...     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

Statistical significance of modal parameters of bridge systems identified from strong motion data

Modal parameters of structural systems have commonly been determined using system identification (SI) methods for damage detection and health monitoring. For determining the deterioration of the integrity of structural systems correctly, modal parameters of a healthy structure have to be obtained with adequate certainty so that these parameters can be used as reliable references for the healthy system to compare with those of the damaged system. In this study, the statistical significance of modal parameters identified using strong motion time histories recorded on two bridge structures is assessed. The confidence intervals of identified modal frequencies and damping ratios are obtained using Monte Carlo simulations and sensitivity analyses in conjunction with eigenrealization algorithm. The dependence of the statistical bounds on model parameters is examined. The effect of using different number of sensors on the statistical significance is evaluated using simulated time history data from a validated finite element model of a bridge. Copyright © 2005 John Wiley & Sons, Ltd.     

二相随机介质有效电导率的广义混合率模型研究(英文)

广义混合率模型目前多应用于岩石的流变、杨氏模量等力学性质的研究,较少应用于多相介质岩石的有效电导率研究。本文利用三维有限元方法计算得到大量二相随机介质模型的有效电导率数据,引入广义混合率的有效电导率模型进行数据拟合,发现广义混合率模型参数J与两相介质电导率比值有关,并首次获得参数 J 与两相介质电导率比值之间的关系式,据此可以快速准确的预测(计算)任意二相介质的有效电导率,其结果较已有的随机介质模型和有效介质理论模型公式更为准确,为精细储集层评价奠定坚实基础。     

Sq的经度效应

本文提供了一个事实:同纬度不同经度台站的S_q曲线可以有很大的差异。对拉萨台的S_q异常作了一些初步的、定性的分析与解释。对八丈台的S_q异常则认为难以解释,值得深入研究。     

Conductivity of NaCl solution at 0.4–5.0 GPa and 25–500 °C

NaCI-H₂O is the most fundamental ternary system in geology. Until now, the measurements of electrical conductivity of NaCl solutions are still little at high pressures (> O.5 GPa) We measured the conductivity of 0.01 m NaCl solution at 0.4–5.0 GPa and 25-500°C. The results are consistent with that of Quist and Marshall (1968) at 0.4 GPa. The conductivity of NaCl solution increases with increasing temperature. The results also show that the conductivity of NaCl solution changes little with increasing pressure below 1.5 GPa and changes rapidly with increasing pressure above 1.5 GPa. The rapid increase of the conductivity of NaCl solution may play an important role in many geological processes (such as the genesis of ore deposits under hydrothermal condition) and other fields. Project supported by the National Natural Science Foundation of China.     

地磁1年周期变化的响应函数

用通日月均值减去夜均值计算的月均值资料和Ｓｏｍｐｉ谱分析法分析了１４个台站的北向分量１年周期变化的振幅向相位，时段为１９８４－１９８８年。分析结果显示用Ｐ＾０ｎ（ｄｎｄ＝１，…，６）复合模式的估算的Ｃ值较合理。各台的Ｃ值比较一致。     

Ressi experimental catchment: Ecohydrological research in the Italian pre-Alps

Ressi is a small (2.4 ha) forested catchment located in the Italian pre-Alps. The site became an experimental catchment to investigate the water fluxes in the soil–plant–atmosphere continuum and the impact of vegetation on runoff generation in 2012. The elevation of the catchment ranges from 598 to 721 m a.s.l. and the climate is humid temperate. The bedrock consists of rhyolites and dacites; the soil is a Cambisol. The catchment is covered by a dense forest, dominated by beech, chestnut, maple, and hazel trees. The field set up includes measurements of the rainfall in an open area, streamflow at the outlet, soil moisture at various depths and locations, and depth to water table in six piezometers at a 5- or 10-min interval. Samples of precipitation, stream water, shallow groundwater and soil water are collected monthly for tracer analysis (stable isotopes (²H and ¹⁸O), electrical conductivity and major ions), and during selected rainfall–runoff events to determine the contribution of the various sources to runoff. Since 2017, soil and plant water samples have been collected to determine the sources of tree transpiration. Data collected in the period 2012–2016 are publicly available. Data collection is ongoing, and the data set is expected to be updated on an annual basis to include the most recent measurements.