ISRN Applied MathematicsVolume 2012 (2012), Article ID 539359, 10 pagesdoi:10.5402/2012/539359

Research Article

Green's Theorem for Sign Data

Louis M. Houston

The University of Louisiana at Lafayette, Lafayette, LA 70504-4210, USA

Received 14 March 2012; Accepted 19 April 2012

Academic Editors: A. Bellouquid and M.-H. Hsu

Copyright © 2012 Louis M. Houston. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Sign data are the signs of signal added to noise. It is well known that a constant signal can be recovered from sign data. In this paper, we show that an integral over variant signal can be recovered from an integral over sign data based on the variant signal. We refer to this as a generalized sign data average. We use this result to derive a Green's theorem for sign data. Green's theorem is important to various seismic processing methods, including seismic migration. Results in this paper generalize reported results for 2.5D data volumes in which Green's theorem applies to sign data based only on traditional sign data recovery.

1. Introduction

In certain cases, data, consisting of coherent signal and random noise, is summed repeatedly in order to enhance the signal and reduce the noise. It has been found that when the signal-to-noise ratio is between 0.1 and 1, if only the signs of the data are retained prior to summation, then the signal can be recovered [1]. We refer to this type of data as sign data. Sign data has advantages in significantly reducing the amount of space needed to record the data. The reduction in space can amount to a ratio of 20 to 1 bits of data storage [1].

In the seismic industry, it has been discovered that most of the data processing methods used on regular data are also effective on sign data [2]. Many of these data processing methods are dependent on Green’s theorem. In particular, Kirchhoff migration is dependent on Green’s theorem. Migration is a seismic process that corrects data coordinates which are distorted by attributes of the seismic experiment. In accordance with Green’s theorem, we show that an integral of an operation on sign data over a volume is equal to an integral of an operation of the original data over a surface.

This result generalizes the work of Houston and Richard [3] in which, based only on traditional sign data recovery, it was shown that Green’s theorem is satisfied when the sign data encompasses a 2.5D volume. 2.5D data is data which is three-dimensional but only has variations along two dimensions. In this paper, we find that Green’s theorem is satisfied for sign data when the data volume is largely symmetric. For cases in which the data volume is of arbitrary shape we have derived a variant of Green’s theorem.

2. A Generalized Sign Data Average

It was shown by O’Brien et al. [1] that there is a recovery of amplitude from sign data in the presence of uniform random noise. If the random noise has amplitude , the signal is designated by the function , and the number of iterations is , then we can write where

Figure 1 illustrates the data recovery process including sign data.

Figure 1: A synthetic signal () and the computer-generated uniform random noise used to examine sign-bit amplitude recovery. This test has noise of unit magnitude () and a signal-to-noise ratio of one. Shown from bottom to top is the signal, the signal plus noise, the sign of the signal plus noise, the average over 200 iterations of signal plus noise, and the average over 200 iterations of the sign of signal plus noise.

Multiply both sides of (2.1) by the function : Now sum both sides of (2.3) over the index: If we allow , then (2.4) becomes or It is clear that if and , then (2.6) becomes which is essentially (2.1). A continuous version of (2.6) might be written as An argument for the consistency of (2.8) is as follows. Let . Then we have Now integrate over all values of : Let be a uniform probability density, .

That implies and (2.10) becomes which was shown by Houston et al. [4].

3. A Generalized Average for Sign Data Derivatives

Consider the following th order forward finite difference [5]: with .

If we make the variable, discrete by choosing a small real interval, and writing where is an integer index, then (3.1) becomes Because the finite difference is a linear operator, we can use (2.6) to derive Equation (3.4) takes into account the fact that where is an arbitrary index.

Since (3.4) suggests the following equation:

Let and .

Then from (3.7) we have Once again, integrate over all , let , and (3.8) becomes which was shown by Houston et al. [4].

4. The Application to Green’s Theorem

Equation (3.7) implies the special case: Employing three variables in (4.1) yields which can be simplified if , , and : Summation over the variables yields Consequently, (4.4) can be written as If we make the mapping then (2.8) becomes Summation over the variables yields which can be written as Subtracting (4.9) from (4.5) yields Green’s theorem is Consequently, we can write a variant of Green’s theorem for sign data as If the spatial intervals are uniform, that is, , , and , then (4.12) becomes In the case of (4.13), we see that Green’s theorem is satisfied by replacing the function with sign data and dividing the surface integral by the noise amplitude. Therefore, all processing of real data based on Green’s theorem will be effective for sign data with an associated variance. We also note that the effectiveness of Green’s theorem on sign data is enhanced by operating over a symmetric volume. When the volume is not symmetric, we can apply the variant of Green’s theorem given by (4.12).

5. The 2.5D Case

Now let us consider the special case for which the functions and have only two-dimensional variation. That is, Without using a generalized sign data average, this implies or Equation (5.3) is Green’s theorem for sign data when the data encompasses a 2.5D volume and is consistent with results reported in Houston and Richard [3].

6. Computational Tests

Equation (2.6) derives from a weighted average of the function . We can thus write the expectation value of as We find that The variance, , can be written as This can be reduced to a simpler form by using the fact that Thus, we have We can demonstrate (2.6), computationally. Let Consequently, we want to demonstrate that We can compute the average percentage error as We should see a correlation between and . The results of this comparison are shown in Table 1.

Table 1: Twenty trials which compare the values of the indicated expressions for iterations per trial. The trials are divided into ten trials for two different noise amplitudes. Recall that . Let and .

3540.594404
3529.425407
(1)

3569.469208
(2)

3460.5978
(3)

3540.394914
(4)

3542.955134
(5)

3644.993387
(6)

3441.978126
(7)

3510.495952
(8)

3544.098984
(9)

3645.64803
(10)

%

3540.594404
3578.540737
(1)

3468.32837
(2)

3063.26568
(3)

3629.793227
(4)

3562.692845
(5)

3605.376332
(6)

3284.005063
(7)

3375.67384
(8)

3356.467161
(9)

3500.608956
(10)

%

7. Conclusions

Using the results of sign data signal recovery leads to a derivation of a generalized sign data average. Extending these results to incorporate derivatives leads to a variant of Green’s theorem for sign data. We find that Green’s theorem directly applies to sign data when the data volume is symmetric and the surface integral is divided by the noise amplitude. A specific application of this result is that Green’s theorem applies to sign data when the data volume is 2.5D and a generalized sign data average is not required.

Acknowledgment

Discussions with Maxwell Lueckenhoff are appreciated.

References

J. T. O'Brien, W. P. Kamp, and G. M. Hoover, “Sign-bit amplitude recovery with applications to seismic data,” Geophysics, vol. 47, no. 11, pp. 1527–1539, 1982.

N. A. Anstey, Seismic Prospecting Instruments, Gebruder Borntraeger, Berlin, Germany, 2nd edition, 1981.

L. M. Houston and B. A. Richard, “The Helmholtz-Kirchoff 2.5D integral theorem for sign-bit data,” Journal of Geophysics and Engineering, vol. 1, no. 1, pp. 84–87, 2004.

L. M. Houston, G. A. Glass, and A. D. Dymnikov, “Sign data derivative recovery,” ISRN Applied Mathematics, vol. 2012, Article ID 630702, 7 pages, 2012.

W. G. Kelley and A. C. Peterson, Difference Equations: an Introduction with Applications, Academic Press, San Diego, Calif, USA, 1st edition, 1991.