NUMERICAL METHOD OF PUMPING TEST ANALYSIS USING MICROCOMPUTERS

Abstract. This paper describes how a numerical method of pumping test analysis, which has proved to be useful in many practical situations, can be run on microcomputers. Full details of a program in BASIC and a test problem are provided. The need to perform all the calculations to a sufficient accuracy is stressed, and the choice of suitable mesh spacings and time steps is discussed.     

用参考台技术计算大地电磁张量阻抗

Gaubau等人于1978年提出的参考台技术,从理论上完全消除了计算张量阻抗时所具有的偏倚误差。但在实际测量中,参考台的距离多大才合适,尚无定论。本文所给的实例,采取了固定主测台、移动参考台的方法。结果表明,参考台距离从20 km变到40 km不会给测量结果带来明显的变化。     

利用地球重力位模型计算重力和重力梯度

高阶高精度地球重力场模型具有广泛的用途。本文利用地球重力位模型计算重力和重力梯度在应用中很有实用阶值,同时也是计算重力场其它量的关键。利用伪局部笛卡尔坐标与球坐标的关系计算了重力与重力梯度在伪局部笛卡尔坐标系下的分量;利用张量变换的原理给出了已知重力与重力梯度在某一坐标系下的分量求它们在另一坐标系下分量的方法,并具体给出了重力与重力梯度在局部笛卡儿坐标系下的分量计算公式,同时还给出计算重力场五参量与垂线偏差的计算公式,本研究推进了地球重力场的可视化进程。     

用参考台技术计算大地电磁张量阻抗

Gaubau等人于1978年提出的参考台技术,从理论上完全消除了计算张量阻抗时所具有的偏倚误差。但在实际测量中,参考台的距离多大才合适,尚无定论。本文所给的实例,采取了固定主测台、移动参考台的方法。结果表明,参考台距离从20 km变到40 km不会给测量结果带来明显的变化。     

PHENOMENA AND CONTROLLING FACTORS IN DRAWDOWN FLUSHING PROCESSES

IINTRODUCTIONByutilizingthenatUralpowerofflow,hydraulicflushingisaneconomicaltechniquetoremovethedepositedsedimentineitherasettlingbasinorareservoir.Inmanyreservoirsaroundtheworld,hydraulicflushinghasbeenprovenasaneffectivemethodtosustainthereservoirstoragecapacity.Hydraulicflushingprocessesmayinvolveboththeprocessesofscouringpreviousdepositsandpassingincomingsediment-ladenflow,suchasgravitycurrent,throughareservoirtobereleased.Inthisarticle,onlytheformercaseisdiscussed.Thekeythinginthef…     

一种利用单台强震记录计算任意场地三分量地面运动的方法

本文利用近场理论格林函数反褶积的方法计算土层的地震反应.解答中包含了产生地面运动的各种反射、透射和转换波震相,也计入了土层与基底介质的相互作用.将计算的土层响应与已获得的近场加速度记录褶积,可算出不同土层场地上垂直向和水平向三分量的地面运动.文中利用1982年卢龙地震的近场加速度记录,检验了方法的可行性.     

一种利用单台强震记录计算任意场地三分量地面运动的方法

本文利用近场理论格林函数反褶积的方法计算土层的地震反应.解答中包含了产生地面运动的各种反射、透射和转换波震相,也计入了土层与基底介质的相互作用.将计算的土层响应与已获得的近场加速度记录褶积,可算出不同土层场地上垂直向和水平向三分量的地面运动.文中利用1982年卢龙地震的近场加速度记录,检验了方法的可行性.     

一种利用单台强震记录计算任意场地三分量地面运动的方法

本文利用近场理论格林函数反褶积的方法计算土层的地震反应.解答中包含了产生地面运动的各种反射、透射和转换波震相,也计入了土层与基底介质的相互作用.将计算的土层响应与已获得的近场加速度记录褶积,可算出不同土层场地上垂直向和水平向三分量的地面运动.文中利用1982年卢龙地震的近场加速度记录,检验了方法的可行性.     

COMPUTING WELL HYDRAULICS SOLUTIONS

Abstract. Type curves, and drawdown analysis curves of other sorts, have a well-established position in the ground-water profession. Although a plethora of solutions to well hydraulics problems exists, many of these arc not available in useful type curve formats. lpurthcr, if one wants to perform parameter estimation using many of these solutions, a computer-based form of the solution is needed. Those who have programmed analytical solutions know the difficulties that will be encountered in coding even the most innocuous-looking result. In this paper, a computer-based method of type curve generation is prcscnted based on Stehfest's method of inverting the Laplace transform solution of ground-water flow problems. The latter are usually easier to obtain than a complete closed-form solution. The method is very robust for a large class of well hydraulics problems. A companion diskette with FORTRAN source and executable example program is available from the author.     

A NOTE ON POSSIBLE USE OF MICROCOMPUTERS FOR AQUIFER STUDIES

Abstract. Microcomputers can be a major tool for the hydrologist. Small and portable, microcomputers can be useful for the acquisition and modeling of an aquifer while on site. This would result in very low cost for computer usage and fast processing time. Their low cost makes the availability almost universal. Tests were run on the TRS-80 system using the Prickett finite difference model and the Theis equation. The results indicate that with the exception of memory storage no problems should arise. With newer generations of microcomputers core memory storage will no longer be a problem.     

岩性地震电算模型技术

根据我国的条件,我们实现了岩性地震电算模型技术。这一技术的理论是建立在二维波动理论的基础上,根据波动方程的克希霍夫积分解,能够研究与介质有关的地震反射振幅特征。其数学公式为     

岩性地震电算模型技术

Here in this paper are presented the results of application of seismic modeling system built on home-made mini-computers of the type 719.     

TWO-DIMENSIONAL MASS DISTRIBUTIONS FROM GRAVITY ANOMALIES*

The oscillatory character of the solution obtained by Qureshi and Mula's method for the direct determination of two-dimensional mass distributions from gravity anomalies is examined.     

超高阶扰动场元的计算方法

空间重力测量技术的飞速发展为超高阶全球重力场模型的建立提供了条件和更迫切的应用需求.本文从完全规格化缔合勒让德函数及其一、二阶导数的标准前向列递推公式出发,研究了超高阶扰动场元球谐展开式中完全规格化缔合勒让德函数及其一、二阶导数的数值特征,改进了上述标准前向列递推关系式,提出了在微机上计算超高阶扰动场元的实用方法,避免了数字溢出.     

GEOELECTRICAL MODEL CALCULATIONS FOR TWO-DIMENSIONAL RESISTIVITY DISTRIBUTIONS*

For the calculation of geoelectrical model curves for a two-dimensional resistivity distribution, the potential equation is transformed by means of a Fourier cosine transform into a two-dimensional Helmholtz equation containing the separation parameter λ. The numerical solution of this equation for different values of λ for an irregular grid is obtained using the method of finite differences combined with the method of overrelaxation. The method by which derivatives are replaced by finite differences turned out to be very important, especially for high resistivity contrasts. After testing several methods designed to deal with any type of resistivity distribution, a method of discretization similar to that used by Brewitt—Taylor and Weaver (1976) for magnetotelluric modeling for H polarization was found the best. Examples are given of model curves for Schlumberger soundings over a vertical fault covered by overburden. The incorrect use of horizontal-layer models leads to erroneous interpretations that are more complex than the real subsurface situations.     

THREE-DIMENSIONAL SEISMIC MODELING FOR ARBITRARY VELOCITY DISTRIBUTIONS *

The growing awareness of the significance of three-dimensional interpretation of seismic reflection data carries with it a need for better understanding of the role of the velocity configuration in the production of reflection patterns. Modeling the reflection responses of a simple dipping plane reflector through velocity models of various degrees of complexity demonstrates the importance of overlying velocities in determining the reflection pattern. Modeling is accomplished using a raytracing technique which determines total travel time of the normal incidence raypath through an arbitrary iso-velocity layer model.