A computing technique is described for routine processing of geophysical field readings A( k) observed in equidistant points along traverses across two-dimensional geological formations. For every point of observation x( k) a number of adjacent readings, N at each side, are included in the calculation. The arithmetical mean of these readings is taken as the temporary regional anomaly A ̄ (k) . The temporary local anomaly A 1 is defined as A(k + i) − A ̄ (k) for ¦ i ¦ ≤ N , and as zero for ¦ i ¦ > N. The local anomaly is approximated by an interpolation polynomial of infinite degree F, being a sum of unit functions E, each multiplied by its corresponding value A 1: F(x)= ∑ i=n N A 1(k+i)·E(λ,ξ,0) ξ= {x−x(k+i)} D E(λ, ξ, 0) = sin ( πξ λ ) πξ in which λ is the degree of smoothing and D the point spacing. If it is assumed that the observed potential field component has constant direction along a straight line any desired component or derivative of the field anywhere in a homogeneous semi-space may be computed with the aid of the potential of the harmonic unit function E( λ, ξ, gz), in which ζ is the normalized vertical variable. Routine calculations are performed by the formula : Q{x(k),z)= ∑ i=−N N C i·A(k+i) where Q is the desired geophysical quantity, and c i are coefficients which are generated by a special subprogram included in a general program for two-dimensional problems. Some applications on model data are demonstrated.
Read full abstract