Land surface subsidence due to excessive groundwater pumping is an increasing concern in California, USA. Interferometric Synthetic Aperture Radar (InSAR) is a remote sensing technique for measuring centimeter-to-millimeter surface deformation at 10–100 m spatial resolution. Here, a data-driven approach that attributes deformation to individual depth intervals within an aquifer system by integrating head data acquired from each of three screened intervals in a monitoring well with InSAR surface deformation measurements was developed. The study area was the Colusa Basin in northern Central Valley. To reconstruct the surface deformation history over the study area, 13 ALOS-PALSAR scenes acquired between 2006 and 2010 were processed. Up to ~3-cm year?1 long-term subsidence and up to ~6 cm seasonal subsidence were observed using the InSAR technique. The technique developed in this paper integrates the InSAR-observed seasonal deformation rate and the co-located head measurements in multiple depth intervals to estimate the elastic skeletal storage coefficient, the time delay between the head change and the observed deformation, and subsequently the deformation of each depth interval. This technique can be implemented when hydraulic head measurements within each depth interval are not correlated with each other. Using this approach, the depth interval that contributed the most to the total subsidence, as well as storage parameters for all intervals, are estimated. The technique can be used for identification of the depth interval within the aquifer system responsible for deformation.
Seismic anisotropy has an important influence on seismic data processing and interpretation. Although the frequency-domain seismic wavefield simulation has a problem of solving the large scale linear sparse matrix due to the computational limitations, it has some advantages over the time-domain seismic wavefield simulation including efficient inversion using only a limited number of frequency components and easy implementation of multiple sources. To accurately simulate seismic wave propagation in the frequency domain, we also need to choose the absorbing boundary conditions to absorb artificial reflections from edges of the model as we do in the time domain. Compared with the classical boundary conditions including the perfectly matched layer and complex frequency-shifted perfectly matched layer, the complex frequency-shifted multi-axial perfectly matched layer has been proven to effectively suppress the unwanted reflections at grazing incidence and solve the instability problem in the time-domain seismic numerical modelling in anisotropic elastic media. In this paper, we propose to extend the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition to the frequency-domain seismic wavefield simulation in anisotropic elastic media. To test the validity of our proposed algorithm, we compare the results (snapshots and seismograms) of the frequency-domain seismic wavefield simulation with those of the time-domain modelling. The model studies indicate that the complex frequency-shifted multi-axial perfectly matched layer absorbing boundary condition is stable in the frequency-domain seismic wavefield simulation in anisotropic media, and provides better absorbing performance than the complex frequency-shifted perfectly matched layer boundary condition. 相似文献
The identification of fractures is of great importance in gravity and magnetic data processing and interpretation.In this study,four fracture identification methods are applied,and widely used in processing and analysis of the gravity anomaly,including vertical second derivative method,tilt derivative method,theta map method and normalized differential method,for gravity data acquired in a given area in Heilongjiang.By comparing the distribution of the zero contour or maximum contour,we summarize the application effects,and both advantages and disadvantages of each method.It is found that tilt derivative method and normalized differential method provide better effects than other two methods:the narrower anomaly gradient belt and higher identification precision of fracture or geological boundary.The inferred fractures and geological boundaries have a great match with the results obtained from geologic map and remote sensing data interpretation.Those study results have definitely provided the theoretical foundation for identifying faults and the geological boundaries. 相似文献