Inverse modelling of aquifer parameters in basaltic rock with the help of pumping test method using MODFLOW software

The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.

Wedge Splitting Test and Inverse Analysis on Fracture Behaviour of Fiber Reinforced and Regular High Performance Concretes

The fracture behaviour of three fiber reinforced and regular HPC （high performance concretes） is presented in this paper. Two mixes are based on optimization of HPC whereas the third mix was a commercial mix developed by CONTEC ApS （Denmark）. The wedge splitting test setup with 48 cubical specimens was used experimentally and the cracked non-linear hinge model based on the fictitious crack model was applied for the interpretation of the results. The stress-crack opening relationships were extracted by using inverse analysis algorithm for various multi-linear softening curves. This showed that the refinement of the softening curves reflects in improved accuracy of the WST （wedge splitting test） simulation in comparison with bi-linear softening curves with acceptable increase of computational time. Furthermore, the fracture mechanics parameters such as COD （crack opening displacement）, fracture energy and characteristic length were experimentally determined. Experiments were performed at 1, 3, 7 and 28 days. Fracture energy, Gf, was found to increase with age, while the characteristic length, Lch, was found to decrease.

Time-Lapse Full-Waveform Inversion Using Cross-Correlation-Based Dynamic Time Warping

Offshore carbon capture, utilization, and storage(OCCUS) is regarded as a crucial technology for mitigating greenhouse gas emissions.Quantitative monitoring maps of sealed carbon dioxide are necessary in a comprehensive OCCUS project. A potential high-resolution method for the aforementioned purpose lies in the full-waveform inversion(FWI) of time-lapse seismic data. However, practical applications of FWI are severely restricted by the well-known cycle-skipping problem. A new time-lapse FWI method using cross-correlation-based dynamic time warping(CDTW) is proposed to detect changes in the subsurface property due to carbon dioxide(CO_(2)) injection and address the aforementioned issue. The proposed method, namely CDTW, which combines the advantages of cross-correlation and dynamic time warping, is employed in the automatic estimation of the discrepancy between the seismic signals simulated using the baseline/initial model and those acquired. The proposed FWI method can then back-project the estimated discrepancy to the subsurface space domain, thereby facilitating retrieval of the induced subsurface property change by taking the difference between the inverted baseline and monitor models. Numerical results on pairs of signals prove that CDTW can obtain reliable shifts under amplitude modulation and noise contamination conditions. The performance of CDTW substantially outperforms that of the conventional dynamic time warping method. The proposed time-lapse fullwaveform inversion(FWI) method is applied to the Frio-2 CO_(2) storage model. The baseline and monitor models are inverted from the corresponding time-lapse seismic data. The changes in velocity due to CO_(2) injection are reconstructed by the difference between the baseline and the monitor models.

Fracture network types revealed by well test curves for shale reservoirs in the Sichuan Basin,China

Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productivity and formulating proper development strategies.This study establishes a new well test interpretation model for fractured horizontal wells based on seepage mechanisms of shale reservoirs and proposes a method for identifying fracturing patterns based on the characteristic slopes of pressure buildup curves and curve combination patterns.The pressure buildup curve patterns are identified to represent three types of shale reservoirs in the Sichuan Basin,namely the moderately deep shale reservoirs with high pressure,deep shale reservoirs with ultra-high pressure,and moderately deep shale reservoirs with normal pressure.Based on this,the relationship between the typical pressure buildup curve patterns and the fracture network types are put forward.Fracturing effects of three types of shale gas reservoir are compared and analyzed.The results show that typical flow patterns of shale reservoirs include bilinear flow in primary and secondary fractures,linear flow in secondary fractures,bilinear flow in secondary fractures and matrix,and linear flow in matrix.The fracture network characteristics can be determined using the characteristic slopes of pressure buildup curves and curve combinations.The linear flow in early secondary fractures is increasingly distinct with an increase in primary fracture conductivity.Moreover,the bilinear flow in secondary fractures and matrix and the subsequent linear flow in the matrix occur as the propping and density of secondary fractures increase.The increase in the burial depth,in-situ stress,and stress difference corresponds to a decrease in the propping of primary fractures that expand along different directions in the shale gas wells in the Sichuan Basin.Four pressure buildup curve patterns exist in the Sichuan Basin and its periphery.The pattern of pressure buildup curves of shale reservoirs in the Yongchuan area can be described as 1/2/→1/4,indicating limited stimulated reservoir volume,poorly propped secondary fractures,and the forming of primary fractures that extend only to certain directions.The pressure buildup curves of shale reservoirs in the main block of the Fuling area show a pattern of 1/4/→1/2 or 1/2,indicating greater stimulated reservoir volume,well propped secondary fractures,and the forming of complex fracture networks.The pattern of pressure buildup curves of shale reservoirs in the Pingqiao area is 1/2/→1/4→/1/2,indicating a fracturing effect somewhere between that of the Fuling and Yongchuan areas.For reservoirs with normal pressure,it is difficult to determine fracture network characteristics from pressure buildup curves due to insufficient formation energy and limited liquid drainage.

Large Deformation Characterization of Porcine Thoracic Aortas: Inverse Modeling Fitting of Uniaxial and Biaxial Tests

The elastic behavior of arteries is nonlinear when subjected to large deformations. In order to measure their anisotropic behavior, planar biaxial tests are often used. Typically, hooks are attached along the borders of a square sample of arterial tissue. Cruciform samples clamped with grips can also be used. The current debate on the effect of different biaxial test boundary conditions revolves around the uniformity of the stress distribution in the center of the specimen. Uniaxial tests are also commonly used due to simplicity of data analysis, but their capability to fully describe the in vivo behavior of a tissue remains to be proven. In this study, we demonstrate the use of inverse modeling to fit the material properties by taking into account the non-uniform stress distribution, and discuss the differences between the three types of tests. Square and cruciform samples were dissected from pig aortas and tested equi-biaxially. Rectangular samples were used in uniaxial testing as well. On the square samples, forces were applied on each side of edge sample attached with hooks, and strains were measured in the center using optical tracking of ink dots. On the cruciform and rectangular samples, displacements were applied on grip clamps and forces were measured on the clamps. Each type of experiment was simulated with the finite element method. The parameters of the Mooney-Rivlin constitutive model were adjusted with an optimization algorithm so that the simulation predictions fitted the experimental results. Higher stretch ratios (>1.5) were reached in the cruciform and rectangular samples than in the square samples before failure. Therefore, the nonlinear behavior of the tissue in large deformations was better captured by the cruciform biaxial test and the uniaxial test, than by the square biaxial test. Advantages of cruciform samples over square samples include: 1) higher deformation range;2) simpler data acquisition and 3) easier attachment of sample. However, the nonuniform stress distribution in cruciform samples requires the use of inverse modeling adjustment of constitutive model parameters.

An Inverse Analysis of the Erichsen Test Applied for the Automatic Identification of Sheet Materials Behavior

Among the technological tests, the Erichsen drawing test gives a more appropriate material behavior, near the limit of the real manufactured process. In this paper an inverse finite element analysis of the Erichsen test is proposed. The new idea is to use a numerical simulation of the experimental test for the rheological identification of the constitutive equations available for sheet metals alloys. The inverse analysis is based on a robust optimization algorithm and uses simultaneously the experimental test data and the corresponding numerical one. A numerical inverse analysis software named OPTPAR was developed and improved for an automatically coupling with a commercial finite element code charged to simulate the experimental test. Results obtained for a virtual steel alloy will be analyzed numerically in order to validate the finite element model and the identification method. An application to an AA5182 aluminum alloy and a DC03 steel alloy will be presented.

A Reinterpretation of Historic Aquifer Tests of Two Hydraulically Fractured Wells by Application of Inverse Analysis, Derivative Analysis, and Diagnostic Plots

Aquifer test methods have greatly improved in recent years with the advent of inverse analysis, derivative analysis, and diagnostic plots. Updated analyses of past aquifer tests allow for improved interpretations of the data to enhance the knowledge and the predictive capabilities of the flow system. This work thoroughly reanalyzes a series of pre- and post-hydraulic fracturing, single-well aquifer tests conducted in two crystalline rock wells in New Hampshire as part of an early 1970’s study. Previous analyses of the data had relied on older manual type-curve methods for predicting the possible effects of hydraulic fracturing. This work applies inverse analysis, derivative analysis, and diagnostic plots to reanalyze the 1970’s aquifer test data. Our results demonstrate that the aquifer tests were affected by changes in flow regimes, dewatering of the aquifer and discrete fractures, and changes due to well development. Increases in transmissivities are related to well development prior to hydraulic fracturing, propagation of a single, vertical fracture hydraulically connecting the two wells after stimulation and expansion of troughs of depression. After hydraulic fracturing, the estimated total yield of the individual wells increased by 2.5 times due to the hydraulic fracturing. However, the wells may be receiving water from the same source, and well interference may affect any significant increase in their combined yield. Our analyses demonstrate the value in applying inverse analysis, derivative analysis, and diagnostic plots over the conventional method of manual type-curve analysis. In addition, our improvement in the aquifer test interpretation of the 1970’s test data has implications for more reliable estimates of sustained well yields.

Inverse Gaussian分布的假设检验

用一种新构造方法，实现了Inverse Gaussian分布的正态化，进而解决了一类长期遗留下来的Inverse Gaussian分布的假设检验问题。

Inverse analysis determining interfacial properties between metal film and ceramic substrate with an adhesive layer

In the present study, peel tests and inverse analysis were performed to determine the interracial mechanical parameters for the metal film/ceramic system with an epoxy interface layer between film and ceramic. Al films with a series of thicknesses between 20 and 250 μm and three peel angles of 90°, 135° and 180° were considered. A finite element model with the cohesive zone elements was used to simulate the peeling process. The finite element results were taken as the training data of a neural network in the inverse analysis. The interracial cohesive energy and the separation strength can be determined based on the inverse analysis and peel experimental result

An inverse problem in film/substrate indentation:extracting both the Young's modulus and thickness of films

In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,only the Young's modulus of the film was unknown.Once the effective Young's modulus of a film/substrate structure is determined at a given contact depth,the Young's modulus of the film can be uniquely determined,i.e.,there is a one-to-one relation between the Young's modulus of the film and the film/substrate effective Young's modulus.However,at times it is extremely challenging or even impossible to measure the film thickness.Furthermore,the precise definition of the layer/film thickness for a two-dimensional material can be problematic.In the current study,therefore,the thickness of the film and its Young's modulus are treated as two unknowns that must be determined.Unlike the case with one unknown,there are infinite combinations of film thickness and Young's modulus which can yield the same effective Young's modulus for the film/substrate.An inverse problem is formulated and solved to extract the Young's modulus and thickness of the film from the indentation depth-load curve.The accuracy and robustness of the inverse problem-solving method are also demonstrated.

Site-testing at Muztagh-ata site Ⅲ:temperature inversion in surface-layer atmosphere

In this article,we present detailed seasonal,monthly and daily statistics of temperature difference in the surface layer at the Muztagh-ata site based on the temperature measurements at two heights of 2 m and6 m.We find that temperature inversion occurs frequently at our site during nighttime,especially during the cold season.Strong temperature inversion always represents stable atmospheric turbulence,which is crucial for an optical observatory.By analyzing the behavior of temperature inversion and its correlation with wind and cloud amount,one conclusion can be made that radiation inversion is the main reason for the existence of temperature inversion in the surface-layer at the Muztagh-ata site.

Inversion of thicknesses of multi-layered structures from eddy current testing measurements

Luquire et al. ' s impedance change model of a rectangular cross section probe coil above a structure with an arbitrary number of parallel layers was used to study the principle of measuring thicknesses of multi-layered structures in terms of eddy current testing voltage measurements. An experimental system for multi-layered thickness measurement was developed and several fitting models to formulate the relationships between detected impedance/voltage measurements and thickness are put forward using least square method. The determination of multi-layered thicknesses was investigated after inversing the voltage outputs of the detecting system. The best fitting and inversion models are presented.

Evaluation of Treadmill Exercise Test Induced U-wave Inversion in Diagnosis of Coronary Artery Disease

U-wave changes during treadmill exercise test were compared between 34 patientswith CAD and 33 normal controls.All of them were confirmed by coronary angiographv.Tenpatients with CAD and 1 normal subject in this study showed U-wave inversion during exercisetest.Of the 10 patients,4 had significant stenosis in one-vessel and 6 in multi-vessel of the ma-jor coronary arteries.The sensitivity,specificity and predictive value of exercise-inducedU-wave inversion were 29％, 97％ and 91％,respectively.It is proposed that exercise inducedU-wave inversion is a reliable predictive index of CAD.

Experimental Studies on High-Frequency Performance of the Inverse Control Magneto-Rheological Damper

Severe vibration of underground structures may be induced under blast loads. According to the characteristics of the explosion-induced ground shock wave, a new-type damper, inverse control magneto-rheological（MR） damper was designed to control the vibration, The high-frequency performance test of the MR damper was carried out on the small shaking table. It is shown that the performance can be modeled by use of the modified Bouc-Wen model, and the Parameters of the model keep stable in the range of 15--50 Hz.

Conductivity Inversion of Unidirectional CFRP Laminates Using Eddy Current Testing

Due to the electrical anisotropy of carbon fiber reinforced polymer(CFRP),this paper presents a method to inverse the anisotropic conductivity of unidirectional CFRP laminate using eddy current testing(ECT). The relationship between the conductivity and probe signal of ECT is studied by means of numerical simulation. Finally,the accuracy of inversion result is improved by optimizing the initial conductivity by use of experimental data.

Inverse Length Biased Maxwell Distribution:Statistical Inference with an Application

In this paper,we suggested and studied the inverse length biased Maxell distribution(ILBMD)as a new continuous distribution of one parameter.The ILBMD is obtained by considering the inverse transformation technique of the Maxwell length biased distribution.Statistical characteristics of the ILBMD such as the moments,moment generating function,mode,quantile function,the coefficient of variation,coefficient of skewness,Moors and Bowley measures of kurtosis and skewness,stochastic ordering,stress-strength reliability,and mean deviations are obtained.In addition,the Bonferroni and Lorenz curves,Gini index,the reliability function,the hazard rate function,the reverse hazard rate function,the odds function,and the distributions of order statistics for the ILBMD,are presented.The ILBMD parameter is estimated using the maximum likelihood method,the method of moments,the maximum product of spacing technique,the ordinary and weight least square procedures,and the Cramer-Von-Mises methods.The Fishers information,as well as the Rényi and q-entropies,are derived.To investigate the usefulness of the proposed lifetime distribution and to illustrate the purpose of the study,a real dataset of the relief times of 20 patients receiving an analgesic is used.

An improved four-dimensional variation source term inversion model with observation error regularization

Aiming at the Four-Dimensional Variation source term inversion algorithm proposed earlier,the observation error regularization factor is introduced to improve the prediction accuracy of the diffusion model,and an improved Four-Dimensional Variation source term inversion algorithm with observation error regularization(OER-4DVAR STI model)is formed.Firstly,by constructing the inversion process and basic model of OER-4DVAR STI model,its basic principle and logical structure are studied.Secondly,the observation error regularization factor estimation method based on Bayesian optimization is proposed,and the error factor is separated and optimized by two parameters:error statistical time and deviation degree.Finally,the scientific,feasible and advanced nature of the OER-4DVAR STI model are verified by numerical simulation and tracer test data.The experimental results show that OER-4DVAR STI model can better reverse calculate the hazard source term information under the conditions of high atmospheric stability and flat underlying surface.Compared with the previous inversion algorithm,the source intensity estimation accuracy of OER-4DVAR STI model is improved by about 46.97%,and the source location estimation accuracy is improved by about 26.72%.

隧道不良地质识别:方法、现状及智能化发展方向

随着隧道施工对于不良地质识别精度要求的不断提高以及人工智能技术的发展,融合多源信息的不良地质智能化识别已成为发展趋势。本文首先阐述了常见的6种隧道不良地质类型及其地质成因,回顾分析了隧道主要的不良地质识别方法及现状,详细介绍了笔者在不良地质智能化识别方面的探索性研究:基于机器学习利用图像识别技术对隧道围岩岩性与裂隙特征进行智能识别;融合图像和光谱特征进行不良地质识别;将地化分析融入到传统的超前钻探中,融合随钻参数和地化信息进行不良地质随钻识别,既可以发挥超前钻探在感知岩体质量和地层信息变化方面的优势,又可以发挥地化分析在岩性和不良地质异常识别方面的优势;基于地质与物探联合反演进行不良地质识别,旨在实现掌子面前方不良地质体“形”(位置、形态、规模)和“性”(性质和类型)的精确识别。最后,对隧道不良地质智能化识别的发展趋势进行了展望。

基于压痕试验的岩石力学参数的卡尔曼滤波反演

压痕试验是间接评估岩体力学特性的重要方法,为了更好地利用压痕试验评估岩体力学参数,本研究采用一种基于扩展卡尔曼滤波算法的参数反演方法。基于室内压痕试验获得的荷载-压入深度曲线数据,通过弹塑性数值仿真模型构建岩体弹塑性参数与压痕曲线的响应关系,进而利用卡尔曼滤波算法实现对岩石弹塑性参数的反演和优化分析。结果表明,本研究方法具有较好的收敛性和较高的识别精度,反演得到的岩石弹性参数均值的相对误差在10%左右,塑性参数的反演结果相对误差小于20%。较传统单一经验公式方法,本研究方法具有更高的预测精度,对多种不同岩性岩体都有较好的适用性。

高水头抽水蓄能电站地应力综合测试及抗劈裂分析

抽水蓄能电站建设处于加速期,合理布置地应力试验和精准获取地应力数据对高水头抽水蓄能电站意义重大。以某高水头抽水蓄能电站为案例,针对高水头抽水蓄能电站地形特点及地应力特征,在不同勘察阶段进行综合地应力测量,同时联合地应力场反演,获取空间三维地应力场分布规律并对围岩进行抗劈裂分析。综合实测结果表明:二维地表深孔结果随埋深梯度有变化,中浅部有应力集中现象,符合高水头抽水蓄能电站的地应力分布特征。三孔交汇法及孔壁应变法的地应力量值及方向吻合,地下厂房和高压岔管部位岩体的最大实测主应力为20.9 MPa,最小主应力为7.0 MPa;岩体以水平应力为主。最大主应力方向集中为NEE向。有利用综合实测结果加强模型精细化进行初始地应力场反演,能较好地反映此类工程的地应力场特征,利于围岩的抗劈裂分析及后续衬砌方案的选择。综合地应力测试对于高水头抽水蓄能的地应力获取具有可行性。可供高水头抽水蓄能电站工程的地应力勘察参考。