TRANSFORMATION OF LINE-SOURCE RESISTIVITY DATA TO POINT-SOURCE DATA AND VICE VERSA USING THE MATRIX METHOD1

1991 ◽  
Vol 39 (1) ◽  
pp. 95-104
Author(s):  
SRI NIWAS
Geophysics ◽  
1975 ◽  
Vol 40 (4) ◽  
pp. 689-693 ◽  
Author(s):  
Sri Niwas ◽  
S. K. Upadhyay

Investigations of apparent resistivity due to a point source over an inclined contact have been reported by Aldredge (1937), Unz (1953), Maeda (1955), and Chastenet de Gery and Kunetz (1956). In these investigations either the image or the harmonic method has been utilized. In this note, we propose to solve the same problem as follows: (1) Transform point‐source potential data into line‐source apparent resistivity data. (2) Interpret transformed apparent resistivities by the master curves provided.


Geophysics ◽  
1978 ◽  
Vol 43 (5) ◽  
pp. 930-942 ◽  
Author(s):  
Irshad R. Mufti

Highly efficient finite‐difference resistivity modeling algorithms which yield accurate results are put forward. The given medium is discretized and divided into rectangular blocks by using a very coarse system of vertical and horizontal grid lines, whose distance from the source(s) increases logarithmically. Expressions are derived to compute the longitudinal conductance and transverse resistance associated with each of these blocks for a parallel‐layer medium followed by a generalized treatment to accommodate arbitrarily shaped structures. Since the values of Dar Zarrouk parameters are derived from the exact resistivity distribution of the given medium, fine features such as a thin but anomalously resistive bed which ordinarily would be missed entirely in coarse discretization can be taken into account. Further reduction in the size of the model is achieved by making use of a symmetry wherever possible. In most cases the computation of the potential field which involves the inversion of a small sparse matrix requires about 0.5 sec of computer time. Moreover, changes in geology affect neither the size nor the zero structure of the matrix. Therefore, when more than one model is to be computed, the factorization of the matrix can be done symbolically only once for all models, followed by numeric factorization for each individual model. The coarse grid algorithm was applied to a number of horizontally layered models involving a point source. The results obtained for each model were in excellent agreement with the corresponding analytical data. Finite‐difference investigation of the potential field for two‐dimensional structures and a line source dipole indicates that as long as one is interested only in the evaluation of the Schlumberger‐type apparent resistivity curves, the line‐source results may be a much better approximation to the corresponding point‐source data than is commonly believed.


Geophysics ◽  
1993 ◽  
Vol 58 (9) ◽  
pp. 1335-1348 ◽  
Author(s):  
Lasse Amundsen

In seismic processing, plane‐wave decomposition has played a fundamental role, serving as a basis for developing sophisticated processing techniques valid for depth‐dependent models. By comparing analytical expressions for the decomposed wavefields, we review several processing algorithms of interest for the geophysicist. The algorithms may be applied to marine point‐source data acquired over a horizontally layered viscoelastic and anisotropic medium. The plane‐wave decomposition is based on the Fourier transform integral for line‐source data and the Hankel transform integral for point‐source data. The involved wavenumber integrals of the cosine or Bessel‐function type are difficult to evaluate accurately and efficiently. However, a number of the processing techniques can easily be run as a filtering operation in the spatial domain without transforming to the wavenumber domain. The mathematical expressions for the spatial filters are derived using plane wave analysis. With numerical examples, we demonstrate the separation of upgoing and downgoing waves from the pressure, the removal of the source ghost from the pressure, and the transformation of point‐source pressure data to the corresponding line‐source data. The filters for these three processes work satisfactorily. Limited spatial aperture is discussed both for point‐source and line‐source data. The resolution kernels relating finite‐aperture decomposed data to infinite‐aperture decomposed data are given. The kernels are approximately equal in the asymptotic limit when the minimum offset is zero.


Geophysics ◽  
2021 ◽  
pp. 1-45
Author(s):  
Qingjie Yang ◽  
Bing Zhou ◽  
Mohamed Kamel Riahi ◽  
Mohammad Al-Khaleel

We present a simple yet effective transform function to convert 3D point-source seismic data to equivalent 2D line-source data, which is required when applying efficient 2D migration and full-waveform inversion to field data collected along a line. By numerically comparing the 3D and corresponding 2D Green’s tensors in various media, the phase shift around 45° and the offset amplitude compensation factor, as well as small fluctuations of the amplitude ratios are observed in all nonzero components of the wave-equation solutions. Based on these observations, we derive a transform function comprised of (1) a simple filter for compensating amplitude and phase shift, and (2) stretching scalars for scaling amplitude differences for different components. We employ the 3D and 2D analytical wave solutions in various homogeneous media to demonstrate the accuracy of the proposed transform function, and then apply it to a heterogeneous, viscoelastic, anisotropic model and a modified Marmousi model. All of these results indicate that the proposed transform function is applicable for the conversion of point-source data to equivalent line-source data for imaging 2D subsurface structure.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Jianbin Tao ◽  
XiangBing Kong

AbstractA gridded social-economic data is essential for geoscience analysis and multidisciplinary application. Spatial allocation of carbon dioxide statistics data is an important issue in the context of global climate change, which involves the carbon emissions accounting and decomposition of responsibility for carbon emission reductions. In this research a new spatial allocation method for non-point source anthropogenic carbon dioxide emissions (ACDE) fusing multi-source data using Bayesian Network (BN) was introduced. In addition to common-used DMSP (Defense Meteorological Satellite Program), PD (population density) and GDP (Gross Domestic Production) data, the land cover and vegetation data was imported into the model as prior knowledge to optimize the model fitting. The prior knowledge here was based on the understanding that ACDE was dominated by human activities and has strong correlations with land cover and vegetation conditions. A 1 km gridded ACDE map integrated emissions form point-source and non-point source was generated and validated. The model predicts ACDE with high accuracies and great improvement can be observed when fusing land cover and vegetation as prior knowledge. The model can achieve successful statistics data downscaling on national scale provided adequate sample data are available, offering a novel method for ACDE accounting in China.


2020 ◽  
Vol 142 (5) ◽  
Author(s):  
Yanan Ding ◽  
Xiaoyan Meng ◽  
Daoyong Yang

Abstract A robust and pragmatic method has been developed and validated to analytically determine dynamic dispersion coefficients for particles flowing in a parallel-plate fracture, in which gravity settling has been considered due to its significant impact on particle flowing behavior. More specifically, a two-dimensional (2D) advection–diffusion equation together with the initial and boundary conditions has been formulated to describe the flow behavior of finite-sized particles on the basis of coupling the Poiseuille flow with vertical settling. Meanwhile, three types of instantaneous source conditions (i.e., point source, uniform line source, and volumetric line source) have been considered. Explicit expressions, which can directly and time-efficiently calculate dynamic dispersion coefficient, have been derived through the moment analysis and the Green’s function method. By performing the simulation based on the random walk particle tracking (RWPT) algorithm, the newly developed model has been verified to determine particle dispersion coefficients agreeing well with those obtained from the RWPT simulations. It is found that the point source is the most sensitive to gravity effect among different source conditions, while the volumetric line source is affected more than the uniform line source. For particle size larger than its critical value, an increased particle size leads to a decreased asymptotical dispersion coefficient for all the source conditions due to the significant gravity effect, while gravity positively affects the dispersion coefficient at early times for the point source condition. In addition, average flow velocity positively affects the dispersion coefficient for all the source conditions, while the associated gravity effect is influenced only at early times for the point source condition.


2012 ◽  
Vol 616-618 ◽  
pp. 804-811
Author(s):  
Quan Tang Fang ◽  
Wei Chen ◽  
Rong Wang

The transient flowing model of slotted liner completion was established by superposition principle based on the geometric model of slotted liners, with the point source function and the single slotting equal to line source, and then the optimized model of slotted liner completion parameter was established with the skin factor of slotted liners completion as evaluation index. After analyzing the parameter sensitivity with cases, the slot density is confirmed as the main reason leading to flow convergence and additional flow resistance. Furthermore, the optimization principles of slotted liners completion of horizontal well are determined. These results are significant in optimizing the slot distribution pattern and parameter allocation.


Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1370-1386 ◽  
Author(s):  
Henry Brysk ◽  
Douglas W. McCowan

The proper implementation of the τ-p method for surface data excited by a point source requires a cylindrical slant stack. Usually the common (Cartesian) slant stack is computed instead as an approximation to the geometrically correct procedure. Here we describe a formulation of the cylindrical slant stack as a weighted sum of Cartesian slant stacks; our cylindrical slant stack is computationally efficient to perform. We show how, although the usefulness of the slant stack is most easily seen with Cartesian coordinates, it can also be used with Fourier‐Bessel transforms. An example of the method shows results computed from data recorded on the West Florida Shelf. Severe edge‐effect noise which overwhelms the Cartesian slant stack is attenuated by the cylindrical slant‐stacking. Applications of the cylindrical slant stack to other seismological calculations, such as Lamb’s problem, are also discussed. In particular, we prove that the plane‐wave reflection coefficients apply exactly in the τ-p domain; hence an amplitude‐versus‐offset analysis is unambiguous in that domain.


Sign in / Sign up

Export Citation Format

Share Document