# Monitoring of Ground Deformation due to Excessive Withdrawal of Natural Gas Using SBAS

Authors: Bo Hu
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The excessive withdrawal of natural gas can generate a wide range of ground deformation. The aim of this paper is to investigate ground deformation in Chiba Prefecture, Japan, using the Small Baseline Subset (SBAS) technique. In recent years, the local government has taken a series of measures such as artificial recharge and reduction of groundwater withdrawal to alleviate the ground settlement. In order to more fully investigate the spatial and temporal variation of deformation that occurred in Chiba Prefecture, we process a dataset of 34 ASAR images acquired from September 2006 to August 2010 based on SBAS technique and generate the mean deformation velocity map and time series deformation maps. The results show that (1) a relatively stable in Mao Yuan area, with an average subsidence velocity of −0.91 mm/y. (2) large uplift rates were highlighted in albino Ding and Kujukuri, which average deformation velocity were about 5.11 mm/y and 6.15 mm/y respectively.


---

## Body

## 1. Introduction

Various crustal movements and anthropogenic activities can result in ground deformation. Crustal movements can result from landslides, earthquakes, volcanic activity, glacial rebound, and subsidence. Anthropogenic causes of displacements include underground resources exploitation [1–5]. However, accurate detection of small deformation over large areas was very difficult in the past. In this background, Multitemporal InSAR has emerged as an effective technique to detect surface displacements with millimeter accuracy. Advanced InSAR technique known as Small Baseline Subset (SBAS) approach has been developed to analyze the temporal evolution of the deformation through the generation of time series of ground displacements [6–10].The formation of Kanto Plain around Tokyo Bay, Japan, is sediment in the oceanic environment through several hundred thousand years or even million years of geological activity. Under this geological activity, Ocean Plankton was transformed into high quality water soluble natural gas [11]. Since the 1930s, there were gas production activities near Tokyo. Large-scale production of gas in Chiba Prefecture area started from 1960. Production was conducted through withdrawal of groundwater from the depth of 500 to 2000 m. Land subsidence was found since the 1960s and large-scale land subsidence surveying was started in 1969. During the past 35 years, the maximum accumulated subsidence was 0.85 m and subsided area was 800 km2 [12]. The coastal line of Chiba Prefecture backed off 30–60 m from its position in the 1960s. In 2005, there is about 15 mm subsidence per year yet. From 2006 to 2010, the local government had taken recharge groundwater to alleviate the settlement. The previously subsidence area had been alleviated, or even uplift [13]. Therefore, accurate detecting and time-series analysis of subsidence change in water-soluble natural gas mining areas, together with the analysis of its driving force, are critically needed for the early hazard warning and sustainable mining.


## 2. Description of Study Area and Data Sets

Chiba Prefecture borders Ibaraki Prefecture to the north at the Tone River, Tokyo, and Saitama Prefecture to the west at the Edo River, the Pacific Ocean to the east and Tokyo Bay around its southern boundary. The geographic extent of Chiba Prefecture is marked by longitudes 140°13′E–114°26′E and latitudes 35°22′N–35°30′N with a total area of around 269.8 km2. Most of Chiba lies on the hilly Boso Peninsula, a rice farming region: the east coast, known as the Kujūkuri Plain, is an especially productive area. The most populous zone, in the northwest of the prefecture, is part of the Kantō region that extends into the urban agglomeration of Tokyo and Saitama [13], the study area shown in Figure 1.Figure 1
The yellow square overlaid on the ASAR intensity map outline of the study area.In this study, the ASAR dataset over the Chiba Prefecture area contains 34 descending acquisitions in IS2 acquisition mode and VV polarization. The temporal coverage is from 2006-09-04 till 2010-08-09 (see Table1).Table 1
Acquisition date of ASAR C-band ASAR data (descending, track: 489) used in this study.
Number
Acquisition date
1
2006/09/04
2
2006/10/09
3
2006/11/13
4
2006/12/18
5
2007/01/22
6
2007/05/07
7
2007/06/11
8
2007/07/16
9
2007/08/20
10
2007/09/24
11
2007/10/29
12
2007/12/03
13
2008/01/07
14
2008/02/11
15
2008/03/17
16
2008/04/21
17
2008/05/26
18
2008/06/30
19
2008/08/04
20
2008/09/08
21
2008/10/13
22
2008/11/17
23
2008/12/22
24
2009/03/02
25
2009/05/11
26
2009/07/20
27
2009/09/28
28
2009/12/07
29
2010/02/15
30
2010/03/22
31
2010/04/26
32
2010/05/31
33
2010/07/05
34
2010/08/09Besides the SAR data, an accurate digital elevation model (DEM), in an ellipsoidal reference system, is prepared before processing the SAR data. If a better DEM is not available, you can use free SRTM [14]. The ASAR scene acquired on Jue. 30, 2008, was used as the super master image. The super master is the reference image of the whole processing, and all the processed slant range pairs will be coregistered on this reference geometry.GivenN acquisitions, the maximum theoretical available interferograms are (N*(N-1))/2. From the whole set of interferograms formed by all possible pairs of images, only those with a perpendicular component of spatial baseline smaller than 350 m, a temporal baseline shorter than 180 days, and a relative Doppler centroid difference below 991 Hz were selected for the following interferometric analysis. By limiting the baselines, the most incoherent interferograms were rejected; the number of all generated interferograms is 213. The temporal baselines of these selected interferograms are shown in Figure 2. From the 213 selected interferograms, 108 interferograms having temporal baselines lower than 180 days and perpendicular baselines less than 100 m were considered as the most significant in terms of the description of the deformation and therefore were stacked (see Figure 3). Spatial phase unwrapping was carried out prior to the stacking analysis. It was accomplished in this work using the MCF phase unwrapping algorithm applied to a triangular irregular network.Figure 2
Time-baseline plots relevant to the ASAR images and interferograms used in this study. The green points represent the ASAR acquisition dates. The yellow point represents the super master image acquisition date.Figure 3
Subset of the whole interferograms that will be unwrapped in 3D. The blue points represent the ASAR acquisition dates. The yellow point represents the super master image acquisition date.


## 3. Principles and Process Flow of the SBAS

SBAS technique was proposed by Berardino et al. in 2002 [7]. This technique is based on the combination of differential interferograms, which is a postprocessing method to generate the mean deformation velocity maps and displacement time series. A detailed discussion on the SBAS algorithm can be referenced in [8].Assuming the existence ofN SAR images, the imaging time was t0,…tN-1. According to the vertical baseline threshold, the above images are grouped. Images within each grouped performs differential interferometric processing, the number of the resulting interferogram K met:
(1)N2≤K≤N(N-12).For a point(x,r) on the jth interferograms, the deformation phases of the two images tA and tB were acquired after removing the topographic phase:
(2)δϕj(x,r)=ϕ(tB,x,r)-ϕ(tA,x,r)≈4πλ[d(tB,x,r)-d(tB,x,r)],iiiiiiiiiiiiiiiiiiiiij=0,…,N-1,
where d(tB,x,r) and d(tA,x,r) are the radar line of sight projections of the cumulative deformation at tB and tA with respect to the first scene t0 assumed as a reference.Equation (2) allows us to define a system of M equations in the N unknowns ϕ(ti,x,r), which can be reorganized, by using matrix formalism, as follows:
(3)Aϕ=δϕ,
where A is an incidence-like matrix directly related to the set of interferograms generated from the available data. Accordingly, the dependence on the variables (x,r) has been neglected in (3) and this simplification will be maintained hereafter in the overall matrix representation.We may now manipulate the system of (3) in such a way to replace the present unknowns with the mean phase velocity between time adjacent acquisitions. Accordingly, the new unknowns become
(4)v=[v1=ϕ(t1,x,r)t1-t0,…,vN=ϕ(tN,x,r)-ϕ(tN-1,x,r)tN-tN-1].And by substituting (4) in (3), we finally get the new equation:
(5)Bv=δϕ,
where B is an M×N matrix corresponding to the unknown vector. In the end, the singular value decomposition (SVD) is applied; you can get the deformation velocity of each time period. At the same time, a filtering step is carried out to remove the atmospheric phase artifacts (APS) and orbital ramps [14–17]. In order to intuitively understand the SBAS, the basic flow of the SBAS is depicted in Figure 4.Figure 4
Basic flow of the SBAS.


## 4. Results and Interpretation

The following section presents results that have been achieved based on SBAS technique. The results are based on processing 34 ASAR data (see Table1). In the process of obtaining differential interferograms, we conducted 1 : 5 multilook in the distance and azimuth direction, respectively, so each pixel size is about 20 m × 20 m. The average surface displacement velocity is shown in Figure 5. The negative sign stands for subsidence, whereas the positive sign is uplift. It can be observed that several significant uplift areas are detected in the natural gas field area, which is shown in Figure 5 by superimposed on the SAR intensity image. It is clear that many ground subsidence funnels have been found around the gas field, in which the maximum subsidence velocity exceeds −11 mm/year in the LOS direction. In the meantime, some uplift areas have been found around some natural gas field after recharge, in which the maximum uplift velocity exceeds 12 mm/year in the LOS direction. According to the filed investigation, these regions are all in or close to the natural gas mining areas. The results showed a relatively stable pattern in Mao Yuan, with an average subsidence velocity of −0.91 mm/year. In contrast, relatively large uplift rates were highlighted in albino Ding and Kujukuri, in which average subsidence velocity was about 5.11 mm/year and 6.15 mm/year, respectively. Additionally, the local government has been taking recharge to alleviate the settlement since 2000. Therefore, we can infer that these land subsidence are induced by the natural gas mining in the underground and ground uplift are induced by water injection from 2006 to 2010.Figure 5
Mean velocity estimate using the SBAS technique. Black line is a hatch through gas field area, and A, B, C, and D are the feature points in the hatch. Red box indicates the area amplified in Figure7. Red triangles represent the location of groundwater wells; black pentagrams represent the location where the water is injected.In order to investigate the development of the ground displacements caused by gas extraction and water injection in Chiba Prefecture, the area located in the gas field, where the location is confined by the box in Figure5, is chosen to exhibit the time series of surface deformation resulting from the SBAS algorithm. As shown in Figure 6, each subplot represents the surface deformation during the period of Sep. 4, 2006, and the date plot in the upper left corner of the subplot. The black region indicates the mask area with low correlation. It is clear that a rebound center has been expanding in scale and aggravating in quantity during about 4 years. The largest rebound occurred in the center of the gas field after water injection for the value of more than 30 mm in the LOS direction since the year of 2006. Besides, a moderate ground rebound is found below the gas field. In addition, in order to intuitively analyze the deformation of the gas field area, the black hatch deformation is shown in Figure 7 and the time series deformation of the feature points is shown in Figure 8. From the time series deformation we can find the ground rebound with the increase of the amount of recharge after the year of 2006.Figure 6
The time series surface deformations at the 33 considered dates with respect to Sep. 4, 2006.Figure 7
The black hatch deformations at Aug. 9, 2010.Figure 8
The time series surface deformations at the 4 feature points (A, B, C, and D) along the hatch with respect to Sep. 4, 2006.


## 5. Conclusions

In this paper, we exploited SBAS technique to investigate the land deformation in LOS direction due to withdrawal of water-soluble natural gas in Chiba Prefecture (Japan) with 34 ASAR data between Sep. 4, 2006, and Aug. 9, 2010. The result in our paper reveals that the area of nature gas mining in Chiba Prefecture has experienced great ground rebound in nearly 4 years, where the maximum accumulation exceeds 30 mm. The presentation of result has much significance to the precaution of gas mining geologic hazard and environment disruption. Due to lack of the ground survey data, we cannot accurately compare the SBAS results currently, but next step is to collect detailed compared ground survey data, for example, leveling, GPS.



---

*Source: 674510-2014-08-27.xml*