共查询到20条相似文献，搜索用时 437 毫秒
1.
Multi-source seismic technology is an efficient seismic acquisition method that requires a group of blended seismic data to be separated into single-source seismic data for subsequent processing. The separation of blended seismic data is a linear inverse problem. According to the relationship between the shooting number and the simultaneous source number of the acquisition system, this separation of blended seismic data is divided into an easily determined or overdetermined linear inverse problem and an underdetermined linear inverse problem that is difficult to solve. For the latter, this paper presents an optimization method that imposes the sparsity constraint on wavefields to construct the object function of inversion, and the problem is solved by using the iterative thresholding method. For the most extremely underdetermined separation problem with single-shooting and multiple sources, this paper presents a method of pseudo-deblending with random noise filtering. In this method, approximate common shot gathers are received through the pseudo-deblending process, and the random noises that appear when the approximate common shot gathers are sorted into common receiver gathers are eliminated through filtering methods. The separation methods proposed in this paper are applied to three types of numerical simulation data, including pure data without noise, data with random noise, and data with linear regular noise to obtain satisfactory results. The noise suppression effects of these methods are sufficient, particularly with single-shooting blended seismic data, which verifies the effectiveness of the proposed methods. 相似文献
2.
In seismic migration, it is important to sample a range of dips around the local structural dip at each image point. Meaningful images are obtained only where this condition holds. For cross-hole seismic reflection surveys, the distribution of dips sampled at each image point is controlled principally by the survey geometry, including source and receiver array lengths and their element spacings. Using a real data set as an example, we show how survey geometry can limit imaging capability close to the boreholes and even in the middle of the section between the boreholes. At the processing stage, effective removal of direct waves and accurate estimation of the velocity field are essential for optimizing image quality. For migration, we propose a generalized Berryhill (GB) scheme which is based on the Kirchhoff integral and takes into account both the near-field and far-field terms. This should improve the ability to image close to source and receiver arrays, provided that the element spacing in the nearby array is small enough. 相似文献
3.
Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays 
下载免费PDF全文
下载免费PDF全文 In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result. 相似文献
4.
Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration 下载免费PDF全文
下载免费PDF全文 In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice. 相似文献
5.
Ambient seismic noise or microtremor observations used in spatial auto-correlation (SPAC) array methods consist of a wide frequency range of surface waves from the frequency of about 0.1 Hz to several tens of Hz. The wavelengths (and hence depth sensitivity of such surface waves) allow determination of the site S-wave velocity model from a depth of 1 or 2 m down to a maximum of several kilometres; it is a passive seismic method using only ambient noise as the energy source. Application usually uses a 2D seismic array with a small number of seismometers (generally between 2 and 15) to estimate the phase velocity dispersion curve and hence the S-wave velocity depth profile for the site. A large number of methods have been proposed and used to estimate the dispersion curve; SPAC is the one of the oldest and the most commonly used methods due to its versatility and minimal instrumentation requirements. We show that direct fitting of observed and model SPAC spectra generally gives a superior bandwidth of useable data than does the more common approach of inversion after the intermediate step of constructing an observed dispersion curve. Current case histories demonstrate the method with a range of array types including two-station arrays, L-shaped multi-station arrays, triangular and circular arrays. Array sizes from a few metres to several-km in diameter have been successfully deployed in sites ranging from downtown urban settings to rural and remote desert sites. A fundamental requirement of the method is the ability to average wave propagation over a range of azimuths; this can be achieved with either or both of the wave sources being widely distributed in azimuth, and the use of a 2D array sampling the wave field over a range of azimuths. Several variants of the method extend its applicability to under-sampled data from sparse arrays, the complexity of multiple-mode propagation of energy, and the problem of precise estimation where array geometry departs from an ideal regular array. We find that sparse nested triangular arrays are generally sufficient, and the use of high-density circular arrays is unlikely to be cost-effective in routine applications. We recommend that passive seismic arrays should be the method of first choice when characterizing average S-wave velocity to a depth of 30 m (Vs30) and deeper, with active seismic methods such as multichannel analysis of surface waves (MASW) being a complementary method for use if and when conditions so require. The use of computer inversion methodology allows estimation of not only the S-wave velocity profile but also parameter uncertainties in terms of layer thickness and velocity. The coupling of SPAC methods with horizontal/vertical particle motion spectral ratio analysis generally allows use of lower frequency data, with consequent resolution of deeper layers than is possible with SPAC alone. Considering its non-invasive methodology, logistical flexibility, simplicity, applicability, and stability, the SPAC method and its various modified extensions will play an increasingly important role in site effect evaluation. The paper summarizes the fundamental theory of the SPAC method, reviews recent developments, and offers recommendations for future blind studies. 相似文献
6.
B. URSIN 《Geophysical Prospecting》1978,26(4):722-749
The use of arrays to separate primary reflections from unwanted coherent seismic events is common practice in land seismic surveys. Very long source and receiver arrays have been used recently to reduce the effects of waterbottom multiples on marine seismic data. The source array consists of five uniformly spaced identical subarrays, each with five different airguns, where the distance between the subarrays may vary from 20 m56 m. The volume of each subarray is 10.3 1 (630 cu.in.) which gives a total volume of the array of 51.5 1 (3150 cu.in.) operated at a pressure of 14 MPa (2000 psi). In order to have a flexible receiver system it was decided to implement the extended receiver array in data processing by computing a weighted sum of two to five traces. The hydrophone cable consists of fifty-four channels with a group length of 50 m. Data shot with the superlong airgun array are processed by a combination of standard techniques and special procedures. In particular, the quality of the stack section is improved by using a weighted stack. The stack weights are computed by a program which takes into account the primary-to-multiple ratio. Comparisons with conventional data show significant improvements in data quality obtained by using the superlong airgun array. Examples show that the waterbottom multiples have been strongly attenuated and the deep seismic events have been enhanced. The combined array response function for dipping events is given in an appendix. 相似文献
7.
A new method for the 2D inversion of induced polarization (IP) data in the time domain has been developed. The entire IP transients were observed and inverted into 2D Cole-Cole earth models, including resistivity, chargeability, relaxation time and the frequency constant. Firstly, a modified 1D time-domain electromagnetic algorithm was used to calculate the response of a layered polarizable ground. The transient signals were then inverted using the Marquardt method to derive the Cole-Cole parameters of each layer. However, model calculations showed that the EM effects could be neglected for the time range (>1 ms) and for the transmitter–receiver distances (<50 m) used in this study. Therefore, the induction effects were not considered for the solution of the 2D inverse problem and a DC solution was applied. An approximative forward algorithm was introduced in order to calculate the IP transients directly in the time domain and in order to speed up the inverse procedure. The approximation is highly accurate, and this is demonstrated by comparing the approximations with their exact solutions up to 3D. The inverse algorithm presented consists of two steps. The transient voltages of an array data set were inverted separately into a two-dimensional resistivity model for each time channel. The time-dependent resistivity of each cell was then interpreted as the response of a homogeneous half-space. In the 2D inversion algorithm, a 3D DC algorithm was used as a forward operator. The method only requires a standard 2D DC inversion and a homogenous half-space Cole-Cole inversion. The developed algorithm has been successfully applied to synthetic data sets and to a field data set obtained from a waste site situated close to Düren in Germany. 相似文献
8.
Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L1‐norm) inversion and smoothness‐constrained least‐squares (L2‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L2‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations. 相似文献
9.
基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27. 相似文献
10.
激发极化法在金属矿、硫化矿等资源勘探方面应用较广.随着勘探设备与计算机硬件的发展,所采集到的数据具有观测数据量大、布极方式多样等特点.针对实测数据的特点,我们研究完成了基于并行技术的激发极化法对数反演算法,该算法具有如下特点：（1）通过压缩存储技术和并行技术的集成实现了可以处理大数据特征的反演算法;（2）利用对数反演来约束每次计算得到新模型的电导率恒为正,充电率在0到1之间变化,从而保证反演的稳定性和可靠性.本文首先设计了两组对比模型进行试验,通过对不同区块数据采用不同加权的方式来减弱噪声对反演结果影响的效果;其次,采用并行技术提高了反演的计算速度,并利用理论模型分析了不同电极装置对反演分辨率的影响.最后,在甘肃某金属矿开发前景区利用激发极化法开展了中梯装置的采集工作,利用加权后的观测数据反演推断出了测区金属矿开发靶区的大致位置及分布特征. 相似文献
11.
Emil Blias 《Geophysical Prospecting》2012,60(6):1054-1067
This paper presents a new explicit method for the estimation of layered vertical transverse isotropic (VTI) anisotropic parameters from walkaway VSP data. This method is based on Dix‐type normal moveout (NMO) inversion. To estimate interval anisotropic parameters above a receiver array, the method uses time arrivals of surface‐related double‐reflected downgoing waves. A three‐term NMO approximation function is used to estimate NMO velocity and a non‐hyperbolic parameter. Assuming the vertical velocity is known from zero‐offset VSP data, Dix‐type inversion is applied to estimate the layered Thomsen anisotropic parameters ?, δ above the receivers array. Model results show reasonable accuracy for estimates through Dix‐type inversion. Results also show that in many cases we can neglect the influence of the velocity gradient on anisotropy estimates. First breaks are used to estimate anisotropic parameters within the walkaway receiver interval. Analytical uncertainty analysis is performed to NMO parameter estimates. Its conclusions are confirmed by modelling. 相似文献
12.
A new approximate method to calculate the space-time acoustic wave motion generated by an impulsive point source in a horizontally layered configuration is presented. The configuration consists of a stack of fluid layers between two acoustic half-spaces where the source and the receiver are located in the upper half-space. A distorted-wave Born approximation is introduced; the important feature of the method is the assumption of a background medium with vertical varying root-mean-square acoustic wave speed. A closed-form expression for the scattered field in space and time as a function of the contrast parameters is deduced. The result agrees closely with rigorously calculated synthetic seismograms. In the inverse scheme the wave speed and mass density can be reconstructed within a single trace. Results of the inversion scheme applied to synthetic data are shown. 相似文献
13.
由于地表电阻率法受到浅地表局部异常体的干扰,反演精度受到影响,井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰,但是由于钻井位置的局限性,数据量得不到保障,也会导致反演精度降低.为此,本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先,利用有限差分法实现正演模拟,采用非线性共轭梯度法（NLCG）恢复电阻率结构;其次,调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算,当开辟12个进程时,反演程序获得最大加速比4.51;最后,通过合成数据和实测数据算例证明该反演算法的有效性,也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果. 相似文献
14.
A. Hering R. Misiek A. Gyulai T. Ormos M. Dobroka L. Dresen 《Geophysical Prospecting》1995,43(2):135-156
For the exploration of near-surface structures, seismic and geoelectric methods are often applied. Usually, these two types of method give, independently of each other, a sufficiently exact model of the geological structure. However, sometimes the inversion of the seismic or geoelectric data fails. These failures can be avoided by combining various methods in one joint inversion which feads to much better parameter estimations of the model than the independent inversions. A suitable seismic method for exploring near-surface structures is the use of dispersive surface waves: the dispersive characteristics of Rayleigh and Love surface waves depend strongly on the structural and petrophysical (seismic velocities) features of the near-surface Underground. Geoelectric exploration of the structure Underground may be carried out with the well-known methods of DC resistivity sounding, such as the Schlumberger, the radial-dipole and the two-electrode arrays. The joint inversion algorithm is tested by means of synthetic data. It is demonstrated that the geoelectric joint inversion of Schlumberger, radial-dipole and two-electrode sounding data yields more reliable results than the single inversion of a single set of these data. The same holds for the seismic joint inversion of Love and Rayleigh group slowness data. The best inversion result is achieved by performing a joint inversion of both geoelectric and surface-wave data. The effect of noise on the accuracy of the solution for both Gaussian and non-Gaussian (sparsely distributed large) errors is analysed. After a comparison between least-square (LSQ) and least absolute deviation (LAD) inversion results, the LAD joint inversion is found to be an accurate and robust method. 相似文献
15.
Lei Li 《Journal of Seismology》2011,15(1):165-170
A so-called teleseismic fluctuation wavefield method (TFWM) was adopted by Shen and Ritter to study small-scale heterogeneities
beneath the Lanzhou array, a small-aperture seismic array locating in the northeastern margin of the Tibetan plateau. Through
this method, they speculated a correlation length of 2.4–8.2 km and a root-mean-squared (RMS) velocity perturbation of 1–3%
for the crust under the array. However, for such a tectonically active and heavily deformed region where great earthquakes
up to M7 occurred, such a low RMS velocity perturbation proposed seems to be unconvincing and conflicts with other studies.
To confirm the suspicion, the acoustic finite difference method is used to model scattering wavefields in the random medium
created through parameters suggested by Shen and Ritter. TFWM is applied to the synthetic data from both the small- and large-aperture
arrays. For synthetic data from arrays with similar aperture as the Lanzhou array, it fails to reproduce the parabolic fit
for real data, implying that the statistical parameters proposed by Shen and Ritter are incorrect. Another test presuming
a correlation length of 2.4 km and a RMS velocity perturbation of 5% can reproduce the fit for real data. The numerical experiments
prove that it is improper to apply TFWM to small-aperture arrays like the Lanzhou array, or else the heterogeneity will be
heavily underestimated. The aperture of an array must be much larger than the correlation length and wavelength to get unbiased
measurements. 相似文献
16.
无线电波透视法是常用的工作面地质构造探测方法之一,目前普遍使用的SIRT方法层析分辨率不高。本文采用约束正则化方法,推导Tikhonov正则化和全变差正则化的最小化问题表达式,讨论影响层析结果的主要因素,对典型理论模型进行了层析成像实验。结果表明:正则化方法具有比SIRT方法更好的分辨率;射线条数越多、噪声水平越低,层析分辨率越高;Tikhonov正则化在正则参数增大时层析结果更光滑,减小时则更贴近异常,全变差正则化与其相反。最后对实际坑透数据进行层析,识别出的异常构造基本吻合已知疑似构造位置,从而说明正则化方法在无线电波透视应用中的可行性。 相似文献
17.
面波成像是研究地壳上地幔横波速度结构的一种重要方法．通常,面波相速度或群速度成像假设面波沿大圆路径传播．但是,在地下介质速度结构变化较大时,面波会偏离大圆路径传播,从而导致基于大圆路径假设下的面波成像结果存在一定的误差. 我们采用基于射线追踪的面波成像方法,研究了面波的偏离大圆路径传播对四川西部地区面波相速度成像结果的影响．使用快速行进法(fast marching method)进行面波传播路径的射线追踪,采用子空间反演法(subspace inversion)进行迭代反演,对理论模型合成数据和川西台阵的短周期背景噪声相速度频散数据进行成像分析,并与使用大圆路径传播的成像结果进行对比．对理论模型的测试结果表明,当速度结构变化较大时,基于偏离大圆路径传播的面波成像能够更好地恢复模型异常．对川西台阵的真实数据反演结果显示:在短周期为6 s时,基于偏离大圆路径传播的反演方法较基于大圆路径传播的反演方法所获得的相速度异常的幅度更大些,在四川盆地区域两者的差异接近0.2 km/s;在周期为10 s时,两种反演方法的差异显著减小,基本都在0.1 km/s以内．这主要是因为6 s周期的面波相速度对复杂的上地壳浅层结构更为敏感,从而使得面波的偏离大圆路径传播效应对反演结果的影响更为显著．本文结果表明,当某一周期不同路径的面波相速度测量值变化较大,例如相对于平均相速度的异常超过10％时,则需考虑采用基于偏离大圆路径传播的面波成像方法,否则速度异常较大区域的反演结果可能会造成较大的偏差． 相似文献
18.
为了对比天然源面波勘探不同台阵布局的探测效果, 筛选出探测成果可靠、 效率高和便于野外施工的天然源面波勘探台阵阵形, 在天水市黄土覆盖区的同一场地分别用4种常见的阵形进行数据采集试验, 并对各种阵形数据使用空间自相关法或扩展空间自相关法提取相应的频散曲线, 通过反演得到了试验点地下的浅层速度结构模型. 分析对比试验结果表明: 4种台阵提取的频散曲线数值很相近; 频散谱能量集中度较高的是嵌套式等边三角形和圆形台阵, L形和直线形台阵相对分散; L形台阵低频段(4—8 Hz)比直线形台阵差, 其高频段(8—40 Hz)比直线形台阵好. 针对直线形台阵在高频段信噪比较低的情况, 在确保探测成果可靠性的前提下, 为了提高探测效率, 提出了在同一直线形台阵开展天然源与人工源面波联合勘探的数据采集方法. 实验结果证实, 这种联合勘探方法不仅可弥补直线形台阵高频段的不足, 确保探测精度和结果的可靠性, 而且还能实现“高低”频兼顾, 即“深浅”兼顾. 相似文献
19.
M.M. ROKSANDIý 《Geophysical Prospecting》1986,34(8):1154-1166
Two factors are responsible for the fact that an extended marine source array performs better than a point source: 1. a higher degree of transmission of the radiated seismic energy through the water-sediment interface providing a better penetration; 2. filtering effects. The higher degree of transmission is due to: (a) the directivity of extended sources, (b) the lower reflection coefficient at the water-sediment interface for seismic waves radiated from an extended source array than for spherical seismic waves radiated from a point source, (c) the lower amplitude decay of the pulses from an extended source than from a point source. In addition, signature characteristic of an extended source array and Fresnel zone of waves generated by such a source differ from those corresponding to a point source. The propagating wavelet radiated from a point source array may not be, in a sedimentary sequence below the sea-floor, the linear combination of wavelets emitted from point sources. In such cases, there is a noticeable difference between the performance of a field-implemented source array and that of the corresponding simulated source array. The performances of the field-implemented and simulated extended receiver arrays can be identical if the recording system is adequate and the processing technique appropriate. 相似文献
20.
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 Kmax/Kmin 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. 相似文献