Iterative deblending with multiple constraints based on shaping regularization |

where denotes the deblended data after th iteration. The blending noise section after th iteration can be expressed as . During the iterations, the blending noise section will contain coherent signal. I propose to iteratively orthogonalize the deblended data and its corresponding noise section, using local signal-and-noise orthogonalization Chen and Fomel (2015). The orthogonalization problem can be expressed as a regularized inversion problem:

where is the local orthogonalization weight (LOW), , which denotes Hadamard (or Schur) product, and denotes the diagonal matrix composed of an input vector. denotes a smoothing regularization operator. When using equation 7, I assume that the signal and blending noise do not correlate with each other. The weighting vector can be solved using shaping regularization Fomel (2007) with a local-smoothness constraint:

where is a triangle smoothing operator, , and is a scaling parameter. The orthogonalized deblended data can be expressed as:

Combining equations 3 and 9 I obtain the following iterative framework using the proposed iterative orthogonalization and seislet thresholding:

The iterative orthogonalization and seislet thresholding framework corresponds to an iterative shaping framework with two projections: thresholding in the seislet transform domain, and the orthogonalization between the deblended data and blending noise. It should be mentioned that the proposed approach can work robust in the case that inappropriate threshold parameters are selected because of an extra denoising compensation during each iteration that minimize the useful energy damages.

h1-1,hs-1
Simulated synthetic example. (a) Unblended data. (b) Blended data.
Figure 2. |
---|

hdbor1-1,hdbst1-1,hdbft1-1,hdbfx1-1
Simulated synthetic example. (a) Deblended data using the proposed approach. (b) Deblended data using seislet domain thresholding. (c) Deblended data using
domain thresholding. (d) Deblended data using
deconvolution.
Figure 3. |
---|

heor1,hest1,heft1,hefx1
Simulated synthetic example. (a) Estimation error section using the proposed approach. (b) Estimation error section using seist domain thresholding. (c) Estimation error section using
domain thresholding. (d) Estimation error section using
deconvolution.
Figure 4. |
---|

h1-z,hs-z,hdbor1-z-0,hdbst1-z-0,hdbft1-z,hdbfx1-z
Zoomed section comparisons of the simulated synthetic example. (a) Unblended data. (b) Blended data. (c) Deblended data using the proposed approach. (d) Deblended data using seislet domain thresholding. (e) Deblended data using
domain thresholding. (f) Deblended data using
deconvolution.
Figure 5. |
---|

hsnrsa,htc
(a) Comparison of signal-to-noise ratio of the simulated synthetic example. "@" refers to the proposed approach. "+" refers to seislet thresholding. "*" refers to
deconvolution. "o" refers to
thresholding. (b) Comparison of the amplitude between 1.35s and 1.38s of the 25th trace in the simulated synthetic data example. Black solid line denotes the unblended trace (true trace). Blue double dot line corresponds to the proposed approach. Red dot dash line corresponds to seislet thresholding. Green dash line corresponds to
thresholding. Yellow long dash line corresponds to
deconvolution.
Figure 6. |
---|

Iterative deblending with multiple constraints based on shaping regularization |

2015-09-15