Cross-matching with interpreted warping of 3D streamer and 3D ocean-bottom-cable data at Valhall for time-lapse assessment

Legacy streamer data and newer 3D ocean‐bottom‐cable data are cross‐matched and analysed for time‐lapse analysis of geomechanical changes due to production in the Valhall Field. The issues relating to time‐lapse analysis using two such distinctly different data sets are addressed to provide an optimal cross‐matching workflow that includes 3D warping. Additionally an assessment of the differences between the imaging using single‐azimuth streamer and multi‐azimuth ocean‐bottom‐cable data is provided. The 3D warping utilized in the cross‐matching procedure is sensitive to acquisition and processing differences but is also found to provide valuable insight into the geometrical changes that occur in the subsurface due to production. As such, this work also provides a demonstration of the use of high‐resolution 3D interpreted warping to resolve the 3D heterogeneity of the compaction and subsidence. This is an important tool for Valhall, and possibly other fields, where compaction and subsidence (and monitoring thereof) are key factors in the reservoir management since the predominant observed production‐induced changes are compaction of the soft, high‐porosity chalk reservoir, due to pore‐pressure reduction, and the resultant overburden subsidence. Such reservoir compaction could have significant implications for production by changing permeabilities and production rates. Furthermore the subsidence effects could impact upon subsea installations and well‐bore stability. Geomechanical studies that have previously been used to model such subsidence and compaction are only constrained by observed surface displacements and measured reservoir pressure changes, with the geological overburden being largely neglected. The approaches suggested herein provide the potential for monitoring and assessment in three dimensions, including the probable heterogeneity and shearing, that is needed for full understanding of reservoir compaction and the resultant effects on the overburden to, for example, mitigate well‐bore failures.     

用重力测量技术观测城市地表下沉的实验研究

本文用重力测量技术对城市地表下沉进行了实验研究,从2016年3月到2017年5月在武汉市内地表下沉较大的部分城区进行了7期流动重力观测实验,并用D-InSAR观测的垂向位移进行了验证.数值结果表明重力观测每期整网平差后点值平均精度都小于10×10~(-8)m·s~(-2),说明采用重力观测能在城市内获得高精度的区域重力变化.第7期相对于第1期的结果与D-InSAR在大致相同时间段内地表垂直位移结果比较表明,重力增加的大部分区域与D-InSAR观测到的地表下沉区域相一致,说明这些区域的重力增加主要是由地表下沉引起的.从第2到7期相对于第1期的重力变化说明在近12个月的时间内测区最大重力变化约40×10~(-8)m·s~(-2),且局部区域的重力值是逐渐增加的,说明地表下沉是持续进行的.本实验结果说明重力观测技术能为城市地表下沉提供重力观测约束和机制解释.     

Multiple-aquifer characterization from single borehole extensometer records

Measurement and analysis of aquifer-system compaction have been used to characterize aquifer and confining unit properties when other techniques such as flow modeling have been ineffective at adequately quantifying storage properties or matching historical water levels in environments experiencing land subsidence. In the southeastern coastal plain of Virginia, high-sensitivity borehole pipe extensometers were used to measure 24.2 mm of total compaction at Franklin from 1979 through 1995 (1.5 mm/year) and 50.2 mm of total compaction at Suffolk from 1982 through 1995 (3.7 mm/year). Analysis of the extensometer data reveals that the small rates of aquifer-system compaction appear to be correlated with withdrawals of water from confined aquifers. One-dimensional vertical compaction modeling indicates measured compaction is the result of nonrecoverable hydrodynamic consolidation of the fine-grained confining units and interbeds, as well as recoverable compaction and expansion of coarse-grained aquifer units. The calibrated modeling results indicate that nonrecoverable specific storage values decrease with depth and range from 1.5 x 10(-5)/m for aquifer units to 1.5 x 10(-4)/m for confining units and interbeds. The aquifer and Potomac system recoverable specific storage values were all estimated to be 4.5 x 10(-6)/m, while the confining units and interbeds had values of 6.0 x 10(-6)/m. The calibrated vertical hydraulic conductivity values of the confining units and interbeds ranged from 6.6 x 10(-4) m/year to 2.0 x 10(-3) m/year. These parameter values will be useful in future management and modeling of ground water in the Virginia Coastal Plain.     

影响四轨法D-InSAR形变测量精度误差的相关性分析

卫星成像的基线、视角、基线倾斜角、斜距、卫星距离地面径向距离以及地面分辨率等因素严重影响着合成孔径雷达差分干涉测量监测地面形变的能力和形变监测结果的精度。本文在分析四轨法D-InSAR基本原理和数据处理流程基础上,详细给出了相位测量误差对形变测量精度影响的定量关系式;分析讨论了基线长度误差、基线倾斜角误差、斜距误差、卫星高度误差和地形因素误差对形变测量精度的影响。从而在定量分析方面得出了这些误差对四轨法D-InSAR形变测量精度影响的结论。     

利用D—InSAR技术研究临汾水准巨幅形变异常

通过D-InSAR技术获取了临汾水准两次巨幅形变异常期间的区域垂直形变场,结果显示：2009年2月20日至2009年10月23日,临汾水准所在的龙祠区域无明显变形,临汾市区南部出现的10～20mm的下沉区域为周边密集分布的多个电厂大量开采地下水所致;2010年1月1日至2010年3月12日,大部分区域变形量基本为零,在临汾盆地内部存在一些小的、空间上不连续分布的沉降区域。幅度约为10mm左右,此类型的沉降区域多与城市发展过程中大规模抽取地下水有关。基于InSAR的结果,表明临汾水准BN,测点的巨幅形变异常在InSAR形变图上无法显示,属于局部地区小范围的变形。     

临汾水准巨幅形变的位错理论模拟与异常性质判定

在假设临汾台水准出现2次巨幅形变异常为断裂错动的情况下,笔者基于矩形断层位错模型模拟了罗云山断裂(土门-峪里段)错动所引起的垂直形变场分布,并通过D-InSAR技术对研究区域内的地面形变场进行了实测。分析结果表明:1从理论上讲,罗云山断裂(土门-峪里段)错动是可以产生长轴与断裂走向平行的椭圆状变形区域,其中,位于断裂上盘的区域中心变形量最大,变形量向外围逐渐衰减为零;2次错动导致的变形波及范围分别约为长轴18km和26km,短轴12km和17km;显著变形幅度分别约为1—3mm和4—14mm。2而同期D-InSAR实测形变场显示,临汾台水准出现2次巨幅异常期间,研究区域内未发现与断裂走向一致的连续变形区域,仅在盆地内部存在可能由于过量开采地下水所导致的地面沉降,其变形范围约为10—12mm和1—5mm。3实测形变场与理论形变场在变形区域和变形幅度上均不一致,说明断裂活动不是临汾台水准出现2次巨幅形变的主要原因,可能为断裂上盘的土层点局部变形所致。4通过断层位错模型的理论模拟与D-InSAR技术的实际监测相结合,可以有效地确定临汾台跨断层水准出现的2次巨幅形变异常的性质,可为重大水准异常的核实提供科学依据。     

Variation in the rate of land subsidence induced by groundwater extraction and its effect on the response pattern of soil microbial communities

Excessive extraction of groundwater leads to (irreversible) changes in the physical soil properties, causing land subsidence associated with soil compaction to occur. Using a combined image processing and field approach, we examined: (1) how variation in the land subsidence rate induces different soil compaction degrees; and (2) the response patterns of microbial communities to such variations. By using Sentinel Synthetic Aperture Radar image processing, we selected three locations that exhibited different land subsidence rates, including high (HSR), moderate (MSR), and low (LSR). Then, soil sampling was undertaken within these representative locations. Indicators of soil compaction, including total porosity, air-filled porosity, water-filled porosity, and bulk density, were measured. The soil microbial community was determined using qPCR and sequencing. The highest and lowest values for bulk density were observed in the HSR–MSR and LSR zones, respectively. The greatest values of total porosity and macropore volume were displayed in the LSR zone compared to other zones. Bacterial abundance in the LSR zone was significantly greater than that in the HSR and MSR zones. The relative abundances of bacterial taxa indirectly demonstrated that the anaerobic phyla were significantly increased (by 10–13%), and the aerobic phyla decreased (by 30–40%) in the HSR zone compared to the LSR zone. This result demonstrates that the aerobes declined as larger volumes of the soil became more anaerobic. Indeed, the increased abundance of anaerobes was not able to compensate for the larger decrease in the abundance of aerobes. Our work showed that at the increased rates of land subsidence, the abundance distribution of the microbial community critically declined. These findings highlight the critical impacts of increasing the land subsidence rate on the emergence of high soil compaction degrees, which can significantly affect the resilience thresholds of the microbial communities in dryland soils.     

星载合成孔径雷达差分干涉测量(D-InSAR)技术在形变监测中的应用概述

本文综述了地球表面形变的主要类型(包括开采沉陷、地表沉降、地壳运动、地震形变、火山运动、冰川运动及山体滑坡等)及其在我国的分布状况，结合合成孔径雷达干涉测量(包括InSAR及D-InSAR，统称InSAR技术)的技术原理及特点，介绍了国内外InSAR技术近年来在形变监测领域的应用与发展。通过与传统形变监测及GPS监测技术的对比后指出，由于InSAR特有的技术特点，使其在各类形变监测应用中具有传统方法无可比拟的技术优势，必将对形变监测的发展起到极大的推动作用。     

拉分盆地形成机制三维数值模拟——以海原断裂带老龙湾盆地为例

拉分盆地是走滑断层系中受拉伸作用形成的断陷盆地.一般在两条平行断层控制下发育.盆地形似菱形,几何形态主要受两条主控走滑断层错距和叠接长度影响.本文以青藏高原东北缘海原断裂带老龙湾拉分盆地第四纪所处的构造环境为基础,参考盆地周围断层几何分布,建立了三维有限元数值模型,模拟该拉分盆地的演化过程;进一步分析了断层力学性质、地壳分层结构等各因素对盆地形成和演化的影响.模拟结果显示,盆地地表沉降伴随有下地壳物质的上涌,此上涌对盆地地表沉降存在阻碍作用.各因素的影响具体表现为:(1)断层力学性质(弹性模量和黏滞系数)越弱,其对构造应力较低的传递效率导致盆地两端差异性运动越明显,从而形成较大的盆地地表沉降和明显的上地壳减薄.(2)平行主控断层的叠接长度反映盆地形成的拉伸作用范围,叠接长度越大,相同的差异性运动在单位面积形成的拉伸应力越小,盆地地表沉降较小.(3)下地壳流变性影响其物质的上涌量,下地壳黏滞系数越小,其对上部拉伸作用的响应越明显,上涌量越大,此上涌对上地壳沉降形成的阻碍作用也越明显.根据老龙湾拉分盆地所处的构造格局,将平行断层的叠接长度取20km,当断层黏滞系数取值为周围基岩的1/10,参考该盆地第四纪构造演化历史,模拟得到的盆地第四纪下沉量与盆地内第四系沉积层厚度在规模上近似,下地壳黏滞系数取值在(2.5~5.0)×1021 Pa·s范围内时,盆地下沉量模拟结果与老龙湾拉分盆地第四系地层厚度吻合较好.     

Subsidence of the Atlantic-type continental margin off New York

The continental margin off New York consists of a thick sequence of at least 4.8 km of mainly shallow-water sediments. Such large thicknesses of shallow-water sediments cannot be produced by the effects of sediment loading alone. We have used biostratigraphic data from the COST B-2 well to examine the origin of the subsidence of this margin. The contribution of sediment loading to the subsidence has been evaluated and removed. Corrections for effects of compaction, water depth and changes of sea-level were also included. The remaining subsidence has been interpreted in terms of a simple thermal model for the cooling lithosphere. Based on this model a thickness of the thermal lithosphere of 113–139 km is estimated. The total subsidence and crustal thinning beneath the COST B-2 well are also estimated and used to place constraints on models for the origin of the margin.     

D-InSAR技术在矿山沉陷和地面沉降监测中的应用

差分雷达干涉(D-InSAR)测量技术,是合成孔径雷达(SAR)卫星应用的一个拓展。雷达图像的差分干涉图可用于监测厘米级或更微小的地面形变,具有全天候、大面积监测地面沉降和矿山沉陷的优势,本文以武安矿山沉陷和沧州市地面沉降监测为例,介绍了这一新技术在灾害监测领域中的实际应用。     

利用InSAR短基线技术估计洛杉矶地区的地表时序形变和含水层参数

美国南加州洛杉矶地区是自然和人为活动引起的地质构造活跃、石油及地下水抽取和回灌频繁的区域.本文利用19景ENVISAT ASAR降轨影像生成了71幅垂直基线小于300 m、时间间隔小于3年的解缠差分干涉图,并基于短基线集技术(SBAS),GPS和地下水水位数据估计了该区域2003年9月~2009年8月的地表时序形变及含水层贮水系数等物理参数.研究结果表明:(1)在InSAR干涉图中可以清楚的识别多处沉降明显的区域.例如,主要由于含水层地下水的抽取与回灌引起地表沉降的Pasadena盆地(~-2.5 cm/a)、San Gabriel流域(~-2 cm/a)、San Bernardino盆地(~-2.5 cm/a)、Pomona-Ontario盆地(~-4 cm/a)和Santa Ana盆地(~-2.5 cm/a),以及由石油抽取引起地面形变的Santa Fe Springs区域(~-1 cm/a)和Wilmington区域(~-1 cm/a)等;(2)InSAR时间序列形变与GPS投影在雷达视线方向上的形变结果具有较高的一致性,平均形变速率差异的均方差为0.39 cm/a;(3)InSAR时间序列形变与含水层地下水位的变化基本一致,并基于相关理论计算出了含水层的弹性贮水系数和非弹性贮水系数,分析了含水层的形变机理.     

L和C波段雷达干涉数据矿区地面沉降监测能力分析

雷达成像的波长、入射角、地面分辨率等参数严重影响着SAR差分干涉测量地面沉降的监测能力和精度,论文通过理论推导和矿区实际沉降差分干涉相位模拟,从监测到的最大沉降梯度和沉降量、保相能力、对微小沉降的敏感程度等方面对L和 C波段雷达干涉数据的矿区地面沉降监测能力进行分析;精化双轨D-InSAR数据处理的流程、方法和相应参数,使用ALOS PALSAR和ENVISAT ASAR数据获取济宁某矿区2009年12月到2010年02月期间更为精确的矿区地面沉降结果,并对沉降结果进行详细比较和系统分析.理论推导、相位模拟和真实数据实验都表明,相对于C波段的雷达干涉数据而言,L波段雷达干涉数据具有较强的保相能力,能够更好地降低失相干和相位不连续性的影响,更容易监测到沉降梯度和沉降量较大的矿区地面沉降,但对微小矿区地面沉降的敏感程度较低.     

2020年1月19日新疆伽师MS6.4地震InSAR同震形变场特征及发震构造初步探讨

基于Sentinel-1 SAR升、降影像,利用D-InSAR技术获取新疆伽师M S6.4地震的同震形变场,结果表明,本次地震引起的同震形变场整体呈近椭圆状分布,形变区东西长约66 km,南北宽约40 km,整个形变场由南部隆升区和北部沉降区组成,南部最大隆升量约7 cm,北部最大沉降量约3 cm。本次地震发生在块体俯冲界面处的低倾角逆冲推覆构造带上,隆升和沉降两个中心均位于逆冲推覆体的上盘,形变主要以隆升形变为主,符合低倾角逆断层中强震的变形特征。在沉降区与隆升区之间干涉条纹连续分布,未出现表征地表破裂位置的空间失相关带,表明地震未引起明显的地表破裂。结合震源机制、余震精定位及区域构造特征,初步推断认为伽师地震的发震构造可能为柯坪塔格推覆构造前缘的N倾的柯坪断裂。     

Land Subsidence Near Oil and Gas Fields,Houston, Texasa

Subsidence profiles across 29 oil and gas fields in the 12,200-km² Houston, Texas, regional subsidence area, which is caused by decline of ground-water level, suggest that the contribution of petroleum withdrawal to local land subsidence is small. Despite large volumes of petroleum production, subsidence at most fields was not increased by oil and gas withdrawal. Local increases of subsidence were detected at only six fields—Alco-Mag, Chocolate Bayou, Goose Creek, Hastings, Mykawa, and South Houston. With the exception of the 1-m subsidence from 1917 to 1925 at Goose Creek, differential subsidence across oil and gas fields was smaller by a factor of two or more than subsidence caused by aquifer compaction. At four fields—Barbers Hill, Cedar Bayou, Humble, and Pierce Junction—subsidence was substantially less than in the surrounding area. Except for Cedar Bayou, these fields are associated with shallow salt domes that partly occupy the aquifer system; for the three fields, subsidence during the periods of record came to less than half the subsidence in the surrounding area. In addition to land subsidence, faults with an aggregate length of more than 240 km (150 mi) have offset the land surface in historical time. Natural geologic deformation, ground-water pumping, and petroleum withdrawal have all been considered as potential causes of the historical offset across these faults. The minor amount of localized land subsidence associated with oil and gas fields, however, suggests that petroleum withdrawal is not a major cause of the historical faulting, at least by a differential compaction mechanism.     

Estimation of the Change in Storage Capacity above Mined Longwall Panels

Accurate estimation of the change in groundwater storage capacity (S) above mined longwall panels is vital for analysis of postmining void water level recovery in coal mines, and assessment of water quality impacts. At present, there is no generalized representation of the spatial distribution of changes in S around a panel. Current estimates are generally bulk averages with high uncertainty, precluding calculation of groundwater velocities in various parts of the subsurface. In this work, a recently published hydrogeological conceptual model of longwall caving is used in conjunction with observations from borehole extensometers, goaf height measurements, and pumping/drawdown records for mine pools to develop a subsurface spatial distribution of changes in S following longwall caving, with reduced uncertainty in their magnitudes. The assumption of saturation in the disturbed zone proved critical for obtaining accurate results and in reconciling widely varying published estimates of S. Results indicate that the goaf and collapsed zones each absorb over 30% of the mined volume, and about 20% is absorbed by the surface subsidence trough. The increase in S in the collapsed zone is inversely proportional to the amount of surface subsidence. The conceptual model is updated with these results to present the spatial distribution of S after caving. The results allow calculation of water velocities in various zones, and may provide greater accuracy in estimation of water level rebound and water quality processes. Most of the S participating in groundwater flows is provided by defects rather than the matrix.     

Determining the dilation factor in 4D monitoring of compacting reservoirs by rock-physics models

Hydrocarbon depletion and fluid injection cause compaction and stretching of the reservoir and overburden layers. 4D prestack seismic data can be used to detect these changes because compaction/stretching causes changes in traveltimes and seismic velocities. We show that, by using two different petro‐elastic models at varying effective pressures, a good approximation is to assume that the fractional changes in layer thickness, ΔL/L, and seismic velocity, Δv/v, are related by a linear function of ΔL/L. The slope of this function (the dilation factor, α= (Δv/v)/(ΔL/L) ) is negative and its absolute value generally decreases (shale, low porosity) or increases (sandstone, high porosity) with increasing layer thickness and decreasing effective pressure. The analysis is mainly performed for isotropic deformations. The dilation factor for uniaxial deformations is smaller in absolute value. The dilation factor, which can be calculated from time‐lapse data, can be used to predict reservoir compaction/stretching as a function of depth and surface subsidence.     

Regional subsidence of Mexico City and its effects on seismic response

Changes expected in Mexico City's seismic response due to land subsidence caused by groundwater extraction are examined. Settlements of the land surface observed in the recent past are compiled and processed to predict regional subsidence in the near future. This information is used in combination with the iso-depth and iso-period maps stipulated by the building code to estimate the compaction of surface sediments and the shortening of soil periods. The time-progressive evolution of the official seismic zoning maps and the associated site-specific response spectra is predicted. The same is done with spectral response contours expressed in terms of the site and structure periods. The subsidence effects are considered through an impact factor, which represents the increase or decrease in the spectral ordinates with respect to the current values specified for design purposes. The detrimental or beneficial effects depend on the relative values of the structure and site periods.     

Subsidence of Surtsey volcano, 1967–1991

The Surtsey marine volcano was built on the southern insular shelf of Iceland, along the seaward extension of the east volcanic zone, during episodic explosive and effusive activity from 1963 to 1967. A 1600-m-long, east-west line of 42 bench marks was established across the island shortly after volcanic activity stopped. From 1967 to 1991 a series of leveling surveys measured the relative elevation of the original bench marks, as well as additional bench marks installed in 1979, 1982 and 1985. Concurrent measurements were made of water levels in a pit dug on the north coast, in a drill hole, and along the coastline exposed to the open ocean. These surveys indicate that the dominant vertical movement of Surtsey is a general subsidence of about 1.1±0.3 m during the 24-year period of observations. The rate of subsidence decreased from 15–20 cm/year for 1967–1968 to 1–2 cm/year in 1991. Greatest subsidence is centered about the eastern vent area. Through 1970, subsidence was locally greatest where the lava plain is thinnest, adjacent to the flanks of the eastern tephra cone. From 1982 onward, the region closest to the hydrothermal zone, which is best developed in the vicinity of the eastern vent, began showing less subsidence relative to the rest of the surveyed bench marks. The general subsidence of the island probably results from compaction of the volcanic material comprising Surtsey, compaction of the sea-floor sediments underlying the island, and possibly downwarping of the lithosphere due to the laod of Surtsey. The more localized early downwarping near the eastern tephra cone is apparently due to greater compaction of tephra relative to lava. The later diminished local subsidence near the hydrothermal zone is probably due to a minor volume increase caused by hydrous alteration of glassy tephra. However, this volume increase is concentrated at depth beneath the bottom of the 176-m-deep cased drillhole.     

Three-Dimensional Numerical Investigation of Pore Water Pressure and Deformation of Pumped Aquifer Systems

Excessive groundwater withdrawal has caused severe land subsidence worldwide. The pore water pressure and the deformation of pumped hydrostratigraphic units are complex. A fully coupled three-dimensional numerical simulation was carried out for different pumping plans in this paper. When groundwater is pumped from a confined aquifer, the great compaction occurs in the pumped aquifer and its upper and lower adjacent aquitard units. Land subsidence is smaller and the area affected by land subsidence is greater when groundwater is pumped from the deeper confined aquifer. The pore water pressure in the pumped confined aquifer changes immediately with pumpage. In the adjacent aquitard units, however, the pore water pressure increases in the early pumping time and decreases in the early recharging time. The decrease in the pore water pressure vertically spreads from the interface between aquitard and pumped aquifer to the other surface of the aquitard. The pumped aquifer compacts and rebounds immediately with pumping and non-pumping or recharging actions, while the compaction and rebounding of the aquitard units clearly lag behind. The compaction of the adjacent aquitard unit first occurs near the interface between aquitard and pumped aquifer units, and the compaction zone spreads outward as the pumping goes on. The aquitards may expand vertically within some zones. Due to the inelastic deformation of soil skeleton, different pumping plans result in different land subsidence. For the same net pumpage, maximal land subsidence and horizontal displacement are the smallest for constant discharge and the greatest for recharge-discharge cycle.