Microearthquake analysis for hydraulic frac－turing proces

Ｍｉｃｒｏｅａｒｔｈｑｕａｋｅａｎａｌｙｓｉｓｆｏｒｈｙｄｒａｕｌｉｃｆｒａｃ┐ｔｕｒｉｎｇｐｒｏｃｅｓＹＩＮＧ－ＰＩＮＧＬＩ（李应平）ＤｅｐａｒｔｍｅｎｔｏｆＥａｒｔｈ，ＡｔｍｏｓｐｈｅｒｉｃａｎｄＰｌａｎｅｔａｒｙＳｃｉｅｎｃｅｓ，Ｍａｓａｃｈ...     

Recognition of a Locked State in Plate Subduction from Microearthquake Seismicity

—A tectonic state of a locked subduction is considered to be a possible source of a future interplate earthquake. Discriminating an actually locked state to verify its extent is therefore essential in constructing an accurate prospect against the forthcoming earthquake. Micorearthquake seismicity is an effective tool for such an analysis because it is considered to be a faithful indicator of the stress state, and is expected to exhibit a characteristic pattern in the area where the locked state in the subduction appears with a certain stress concentration. Focusing on the microearthquake seismicity around the Tokai district in central Japan, where a large interplate earthquake is feared to occur, we tried to identify such an area of locked subduction on the Philippine Sea plate, possibly related to the future earthquake. We investigated the microearthquake seismicity from various perspectives. First, the hypocenter distribution was analyzed to identify the extent of the locked area. The characteristic profile of the distribution was presumed to represent a stress concentrated area induced from the mechanical contact between both plates. The second approach is to interpret stress patterns reflected in focal mechanisms. The locked state was recognized and verified by a comparison of the P-axis distribution pattern with that expected from a model imaging a partially locked subduction. The third approach is to monitor the temporal change of the seismic wave spectrum. Analyzing predominant frequencies of P and S waves and monitoring their changes for a period of 10 years, we found a trend of gradual increase common to both waves. This means an increase of stress drop in microfracturings, and in its turn implies accumulation of stress around the focus area. The rate of the stress change converted from the frequency change was compared with the result derived from a numerical simulation. The simulation, performed on the basis of a constitutive friction law for a stick sliding on the plate interface, computed a changing rate of the maximum shear stress around the locked zone and showed its spatial variation along the subduction axis. Thus the simulated result indicated a certain compatibility with the observed one. Although ambiguities and uncertainties still exist in the study, all the results derived here seem to indicate an identical conclusion that the plate subduction is actually locked in this region at present.     

Microearthquake observations around Thien Dam in northwest Himalayas

A microearthquake survey has been conducted around Thien Dam (32°26.5′N, 75°43.5′E) during March–April 1983 by operating seven high gain stations. The data recorded at these stations have been supplemented by other seismological stations in the region under the Beas and Salal Projects. Surficial velocity structure has been obtained from waves generated from blasts in the dam area and hypocentres were determined using the program HYPO 71.Microearthquake distribution shows four main zones of active seismicity. Deeper microearthquakes with focal depths greater than 10 km below the dam were attributed to the Satlitta fault. Composite focal mechanism solution shows that the Satlitta fault at this depth is of normal type. Shallow microearthquakes below the dam are due to the Suruin-Mastgarh anticline which is of the thrust type. Microearthquakes located northeast of the dam are apparently associated with a ‘back’ thrust.     

Microearthquake seismicity and active tectonics of northwestern Greece

We carried out a microearthquake survey lasting for six weeks in northwest Greece using 18 portable seismograph stations to examine a region in which normal and thrust faulting have been reported in close proximity to one another. With this array we located 148 events and determined fault plane solutions for eight events using only rays radiated upwards. The seismicity of the region is diffuse with events extending to depths of nearly 30 km, and there is a minimum in activity near a depth of 15 km. The fault plane solutions exhibit a wide spectrum of fault types and orientations and are not consistent with simple zones of shortening or extension. Neither tractions applied to the edges or bottom of the region nor deviatoric stresses that compensate for lateral variations in crustal thickness can account for the variety of fault plane solutions. We think that the complicated behavior is a manifestation of inhomogeneous deformation due, at least in part, to a pre-existing complicated juxtaposition of structures and formations.     

Microearthquake studies in Egypt carried out by the geological survey of Egypt

Extensive microearthquake studies have been conducted in Egypt as a joint project between scientists from the Egyptian Geological Survey and Mining Authority (EGSMA) and U.S. scientists. At this stage, a great part of the data has been analyzed and two intensively active areas have been located: one in the Abu Dabbab area of the Eastern Desert, the second at the mouth of the Gulf of Suez near Gubal Island (Daggett et al., 1980). Both sites have been reported to be the epicenters of large earthquakes in 1955 and 1969, respectively. A few scattered earthquakes have also been located in the northern part of the Red Sea, some of which lie along its median axis (Daggett et al., 1986) adding to evidence for the medial opening of the northern Red Sea. After the occurrence of an earthquake (M = 5.5) in the Aswan region on 14 November 1981, continuous recording of the many aftershocks was carried out by EGSMA for about seven months from December 1981 to July 1982, when the temporary network was replaced by a network of telemetered seismographs installed and operated by Helwan Institute of Astronomy and Geophysics in cooperation with scientists from Lamont and Doherty Geological Observatory (LDGO). The majority of epicenters are concentrated in the vicinity of G. Marawa about 65 km upstream of Aswan Dam, along the E-W Kalabsha fault. The observed focal mechanism is consistent with a right-lateral strike-slip motion on the Kalabsha fault. Analysis of Aswan microearthquakes has been done by EGSMA in cooperation with scientists from California Division of Mines and Geology (CDMG).     

High-resolution modelling of a large-scale river plume

Evolution of a large-scale river plume is studied numerically using the Massachusetts Institute of Technology general circulation model. The model parameters were set close to those observed in the area of the Columbia River mouth. The fine-resolution grid along with the non-hydrostatic dispersion included in the model allowed for the reproduction of detailed inner plume structure, as well as a system of internal waves radiated from the plume’s boundary. It was found that not only first-mode but also second- and third-mode internal waves are radiated from the plume at the latest stages of its relaxation when the velocity of the front propagation drops below an appropriate wave phase speed of internal baroclinic mode. The model output shows that the amplitude of these high-mode waves is of the same order as the leading first-mode waves, which in combination with the specific vertical structure (location of the maximum structure function beyond the pycnocline layer) creates favourable conditions for the generation of shear instabilities. High-resolution model output also reveals evidence of a fine internal structure of the plume characterized by the presence of secondary fronts inside the plume and secondary internal wave systems propagated radially from the lift-off area to the outer boundary. These structures intensify the mixing processes within the propagating plume with predominance of the entrainment mechanism developing on the lower boundary between the plume’s body and underlying waters. The scheme of horizontal circulation in the plume was reproduced by the methodology of Lagrange drifters released near the mouth at different depths.

     

High-resolution modelling of a large-scale river plume

Evolution of a large-scale river plume is studied numerically using the Massachusetts Institute of Technology general circulation model. The model parameters were set close to those observed in the area of the Columbia River mouth. The fine-resolution grid along with the non-hydrostatic dispersion included in the model allowed for the reproduction of detailed inner plume structure, as well as a system of internal waves radiated from the plume’s boundary. It was found that not only first-mode but also second- and third-mode internal waves are radiated from the plume at the latest stages of its relaxation when the velocity of the front propagation drops below an appropriate wave phase speed of internal baroclinic mode. The model output shows that the amplitude of these high-mode waves is of the same order as the leading first-mode waves, which in combination with the specific vertical structure (location of the maximum structure function beyond the pycnocline layer) creates favourable conditions for the generation of shear instabilities. High-resolution model output also reveals evidence of a fine internal structure of the plume characterized by the presence of secondary fronts inside the plume and secondary internal wave systems propagated radially from the lift-off area to the outer boundary. These structures intensify the mixing processes within the propagating plume with predominance of the entrainment mechanism developing on the lower boundary between the plume’s body and underlying waters. The scheme of horizontal circulation in the plume was reproduced by the methodology of Lagrange drifters released near the mouth at different depths.     

Induced Microearthquake Patterns in Hydrocarbon and Geothermal Reservoirs: Six Case Studies

—?The injection or production of fluids can induce microseismic events in hydrocarbon and geothermal reservoirs. By deploying sensors downhole, data sets have been collected that consist of a few hundred to well over 10,000 induced events. We find that most induced events cluster into well-defined geometrical patterns. In many cases, we must apply high-precision, relative location techniques to observe these patterns. At three sedimentary sites, thin horizontal strands of activity are commonly found within the location patterns. We believe this reflects fracture containment between stratigraphic layers of differing mechanical properties or states of stress. At a massive carbonate and two crystalline sites, combinations of linear and planar features indicate networks of intersecting fractures and allow us to infer positions of aseismic fractures through their influence on the location patterns. In addition, the fine-scale seismicity patterns often evolve systematically with time. At sedimentary sites, migration of seismicity toward the injection point has been observed and may result from slip-induced stress along fractures that initially have little resolved shear. In such cases, triggering events may be critical to generate high levels of seismic activity. At one crystalline site, the early occurrence of linear features that traverse planes of activity indicate permeable zones and possible flow paths within fractures. We hope the continued development of microseismic techniques and refinement of conceptual models will further increase our understanding of fluid behavior and lead to improved resource management in fractured reservoirs.     

甘东南流动台阵微震监测结果

利用高密度宽频带流动台阵资料,应用网格搜索地震定位方法对甘东南2009年11月-2011年11月期间的微震进行了重新定位.在研究区内有3个以上台站记录到的微震共3 014个,比常规台网记录多1 693个.重新定位结果显示:在各断裂带附近地震分布更加密集,沿断裂呈带状或簇集状分布;震源深度分布优势范围为5～20 km,其平均深度为11 km;地震深度分布与断裂的关系密切,地震条带与断裂构造可以较好的相对应.     

Active faulting of the southeastern-most Zagros (Iran): Microearthquake seismicity and crustal structure

The microseismicity of the southeastern-most Zagros is examined by high-resolution data recorded by a temporary dense local seismic network. The seismicity defines a diffuse pattern, mostly located beneath folds in the southern part of the High Zagros Fault (HZF). Seismicity dips gently northward in the depth range 6–25 km, implying slip on a major intracrustal thrust fault extending to the north of the Main Zagros Reverse Fault (MZRF) which seems to connect to the Mountain Frontal Fault (MFF). Furthermore, observed focal mechanisms suggest transpressive motion on the HZF located west of the Zendan-Minab-Palami (ZMP) fault system and striking obliquely to the convergent motion. These observations suggest that the transition zone between the Zagros continental collision zone and the Makran oceanic subduction zone is not confined to the east of the ZMP and some part of the this diffuse transition is accommodated north of the Hormuz Strait in the west by partitioning between strike-slip and shortening components. The Zagros reverse domain is terminated by a transpressive tectonic regime. Moho depth beneath the MZRF, deduced from receiver functions, is almost 45 km thinner than is observed in the central and northern parts of the Zagros. These observations support a model of active underthrusting of the Arabian plate beneath central Iran in the southeastern-most Zagros.     

High-resolution slug testing

The hydraulic conductivity (K) variation has important ramifications for ground water flow and the transport of contaminants in ground water. The delineation of the nature of that variation can be critical to complete characterization of a site and the planning of effective and efficient remedial measures. Site-specific features (such as high-conductivity zones) need to be quantified. Our alluvial field site in the Kansas River valley exhibits spatial variability, very high conductivities, and nonlinear behavior for slug tests in the sand and gravel aquifer. High-resolution, multilevel slug tests have been performed in a number of wells that are fully screened. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been used to analyze the data, employing an automated processing system that runs within the Excel spreadsheet program. It is concluded that slug tests can provide the necessary data to identify the nature of both horizontal and vertical K variation in an aquifer and that improved delineation or higher resolution of K structure is possible with shorter test intervals. The gradation into zones of higher conductivity is sharper than seen previously, and the maximum conductivity observed is greater than previously measured. However, data from this project indicate that well development, the presence of fines, and the antecedent history of the well are important interrelated factors in regard to slug-test response and can prevent obtaining consistent results in some cases.     

Microearthquake survey in the hot spring area in Parbati Valley (Himachal Pradesh, India)

A number of hot springs occur in the Parbati Valley in Himachal Pradesh in India. Temperatures range from 21 to 96°C, the boiling point of water at that altitude. Geological conditions, temperature variations and chemical composition of spring water in the Parbati Valley hydrogeological unit indicate that the deep thermal fluids are of meteoric origin. The maximum temperature acquired by water during its circulation is estimated to come close to 200°C. In order to assess the possibility of extracting geothermal energy, a seismic survey was arranged to locate the hypocentres of microearthquakes associated with the thermal source. A total of eight microearthquake units was set up at an interstation spacing of about 10 km and two months recording were obtained. During this period an average of 2–3 events per day was recorded with S—P interval less than 5 seconds. The data have been analysed with the help of Hypo 71, a Fortran IV computer program designed to determine the hypocentral parameters of earthquakes from seismic data. The results indicate faulting but there is no apparent spatial relationship to surface manifestations of geothermal energy.     

地震序列中“转折时间t0”的力学意义及作为前兆性特征参量的探讨

本文提出了一种利用地震序列“转折时间ｔ０”来判断震源区附近介质粘性系数η相对大小的方法。η的相对大小一方面反应了地震序列持续时间的长短，另一方面又隐含着地震类型的信息，于是，本文探讨了将ｔ０作为分析地震序列持续时间和估计震源区附近发生后续强震的前兆性特征参可能性。     

High-resolution Electrical Resistivity Tomography monitoring of a tracer test in a confined aquifer

A permanent geoelectrical subsurface imaging system has been installed at a contaminated land site to monitor changes in groundwater quality after the completion of a remediation programme. Since the resistivities of earth materials are sensitive to the presence of contaminants and their break-down products, 4-dimensional resistivity imaging can act as a surrogate monitoring technology for tracking and visualising changes in contaminant concentrations at much higher spatial and temporal resolution than manual intrusive investigations. The test site, a municipal car park built on a former gasworks, had been polluted by a range of polycyclic aromatic hydrocarbons and dissolved phase contaminants. It was designated statutory contaminated land under Part IIA of the UK Environmental Protection Act due to the risk of polluting an underlying minor aquifer. Resistivity monitoring zones were established on the boundaries of the site by installing vertical electrode arrays in purpose-drilled boreholes. After a year of monitoring data had been collected, a tracer test was performed to investigate groundwater flow velocity and to demonstrate rapid volumetric monitoring of natural attenuation processes. A saline tracer was injected into the confined aquifer, and its motion and evolution were visualised directly in high-resolution tomographic images in near real-time. Breakthrough curves were calculated from independent resistivity measurements, and the estimated seepage velocities from the monitoring images and the breakthrough curves were found to be in good agreement with each other and with estimates based on the piezometric gradient and assumed material parameters.     

2020年古冶MS5.1地震序列微震检测及活动性研究

2020年7月12日河北省古冶(39.78°N, 118.44°E)发生M_S5.1地震, 震源深度10 km.本文通过对其周边500 km范围内的84个固定台站连续观测地震数据, 采用基于图形处理器的模板匹配定位识别技术(GPU-M&L)对地震前后(2020年6月1日至2020年8月31日)共三个月的连续数据进行了微震检测和地震活动性分析.我们首先从台网目录的362个古冶地震序列事件中挑选了信噪比较高的118个地震作为模板, 通过模板扫描和检测, 获得共1017个地震, 约为台网目录的2.8倍, 然后运用双差定位方法(hypoDD)对669个地震事件进行了重定位.最后我们对古冶地震的事件序列进行时空分布、震源机制解以及b值空间分布等分析.结果显示古冶M_S5.1地震前, 震中区周边地震活动性显著, 地震前后地震序列分布优势方向为NE向, 震源深度方向近直立分布; 发震断裂为唐山—古冶断裂或其延伸部分, 与1976年唐山M_S7.8地震同为唐山断裂带, 发震断层走向基本一致.据此, 我们推测古冶M_S5.1地震是1976年唐山M_S7.8地震后余震区应力调整所触发.

     

一种高分辨率温度传感器

介绍一种基于Pt1000惠斯通平衡电桥的高分辨率温度传感器。 实验测试表明, 该温度传感器具有精确、 稳定、 分辨率高的特点, 测试温度的分辨率可以达到1×10^-5℃。 其温度测试范围可调, 从而可以满足不同测试环境温度测试的需要, 具有较广泛的应用范围, 并可同通用数据采集器连接构成精密温度观测系统。     

Mechanics and performance of a tied-back wall under seismic loads

In this paper the stability of a tied-back wall subjected to seismic loads is analysed for a predetermined mode of failure (rotation about the top of the wall) and the analysis is compared with data from tests on this type of wall using the seismic simulator at the State University of New York at Buffalo. We carried out a pseudo-static analysis of the problem using the Mononobe-Okabe earth pressure coefficients, wherein the dynamic effects due to the seismic loading are converted into equivalent static loads. The acceleration ratio at which the wall fails by rotation about the top was obtained by considering the moments due to the various lateral earth pressure resultants and the inertial forces induced in the soil due to the seismic loading. We found that the presence of wall friction on the passive side significantly enhances the stability of the flexible retaining wall under seismic loads. Thus, flexible retaining walls supporting dry cohesionless soil can be very efficient during earthquakes. Under moderate earthquakes, an increase in the depth of embedment increases the dynamic factor of safety significantly. However, beyond a certain acceleration ratio for a soil with a particular value of ø, any increase in the depth of emdedment has no effect in impeding failure, irrespective of any change in the geometry of the system. Seismic design charts are presented to evaluate the stability of, and to design, flexible retaining walls embedded in dry cohesionless soils under seismic loading.     

Mechanics of Sediment Transport

After12yearssinceProfDr.NingChienpassedawayin1986,theEnglishversionofhismasterpiece"MechanicsofSedimentTransport"wasfinallypublishedbyASCEpressin1998.Thebookisimmeadialelyfeaturedunder"ASCEBestSellingTit-tles"andwelcomedbyhydraulicengineersandscientists.ThereforethemarketingassistantofASCE,Mr.BestyShepardwroteMrs.ChienforC0ngratulation.Dr.Chienspent30yearsforthebook.Intheprefacehewrote:"ThewritingofthisbookwasbeganwhileIwasstillatanageintheprimeoflifeandendedattheagewithgreyingtem…     

High-resolution computed microtomography for the characterization of a diffusion tensor imaging phantom

This paper addresses the issue of the quantitative characterization of the structure of the calibration model (phantom) for b-matrix spatial distribution diffusion tensor imaging (BSD-DTI) scanners. The aim of this study was to verify manufacturing assumptions of the structure of materials, since phantoms are used for BSD-DTI calibration directly after manufacturing. Visualization of the phantoms’ structure was achieved through optical microscopy and high-resolution computed microtomography (µCT). Using µCT images, a numerical model of the materials structure was developed for further quantitative analysis. 3D image characterization was performed to determine crucial structural parameters of the phantom: porosity, uniformity and distribution of equivalent diameter of capillary bundles. Additionally calculations of hypothetical flow streamlines were also performed based on the numerical model that was developed. The results obtained in this study can be used in the calibration of DTI-BST measurements. However, it was found that the structure of the phantom exhibits flaws and discrepancies from the assumed geometry which might affect BSD-DTI calibration.