共查询到20条相似文献，搜索用时 46 毫秒
1.
Hydraulic fractures generated by fluid injection in rock formations are often mapped by seismic monitoring. In many cases, the microseismicity is asymmetric relative to the injection well, which has been interpreted by stress gradient along the direction of the hydraulic fracture. We present a mathematical model of asymmetric hydrofracture growth based on relations between the solid‐phase stress and the fracture hydraulics. For single fracture and single injection point, the model has three parameters, hydraulic conductivities of the fracture wings, and normalised stress gradient and predicts the positions of the fracture tips as functions of time. The model is applied to a set of microseismic event locations that occurred during and after an injection process. Two different methods are suggested that make it possible to delineate the fracture tips from the set of microseismic events. This makes it possible to determine the model parameters and to check the agreement between the model prediction and the measured data. The comparison of the measured and modelled growth of fracture wings supports both the assumption of the non‐zero stress gradient and the existence of the post‐injection unilateral growth. 相似文献
2.
Numerical evaluation of the effect of fracture network connectivity in naturally fractured shale based on FSD model 总被引：1，自引：0，他引：1
In this paper, the effect of pre-existing discrete fracture network (DFN) connectivity on hydraulic fracturing is numerically investigated in a rock mass subjected to in-situ stress. The simulation results show that DFN connectivity has a significant influence on the hydraulic fracture (HF) & DFN interaction and hydraulic fracturing effectiveness, which can be characterized by the total interaction area, stimulated DFN length, stimulated HF length, leak-off ratio, and stimulated total length. In addition, even at the same fluid injection rate, simulation models exhibit different responses that are strongly affected by the DFN connectivity. At a low injection rate, total interaction area decreases with increasing DFN connectivity; at a high injection rate, total interaction area increases with the increase of DFN connectivity. However, for any injection rate, the stimulated DFN length increases and stimulated HF length decreases with the increase of connectivity. Generally, this work shows that the DFN connectivity plays a crucial role in the interaction between hydraulic fractures, the pre-existing natural fractures and hydraulic fracturing effectiveness; in return, these three factors affect treating pressure, created microseismicity and corresponding stimulated volume. This work strongly relates to the production technology and the evaluation of hydraulic fracturing effectiveness. It is helpful for the optimization of hydraulic fracturing simulations in naturally fractured formations. 相似文献
3.
In unconventional reservoirs, small faults allow the flow of oil and gas as well as act as obstacles to exploration; for, (1) fracturing facilitates fluid migration, (2) reservoir flooding, and (3) triggering of small earthquakes. These small faults are not generally detected because of the low seismic resolution. However, such small faults are very active and release sufficient energy to initiate a large number of microseismic events (MEs) during hydraulic fracturing. In this study, we identified microfractures (MF) from hydraulic fracturing and natural small faults based on microseismicity characteristics, such as the time–space distribution, source mechanism, magnitude, amplitude, and frequency. First, I identified the mechanism of small faults and MF by reservoir stress analysis and calibrated the ME based on the microseismic magnitude. The dynamic characteristics (frequency and amplitude) of MEs triggered by natural faults and MF were analyzed; moreover, the geometry and activity types of natural fault and MF were grouped according to the source mechanism. Finally, the differences among time–space distribution, magnitude, source mechanism, amplitude, and frequency were used to differentiate natural faults and manmade fractures. 相似文献
4.
《中国科学:地球科学(英文版)》2021,(3)
Reservoir reconstructions implemented in unconventional oil and gas exploration usually adopt hydraulic fracturing techniques to inject high-pressure fluid into the reservoir and change its pore-fracture connection structure to enhance production. Hydraulic fracturing changes the reservoir stress and causes the rocks to crack, thus generating microseismic events.One important component of microseismic research is the source mechanism inversion. Through the research on the microseismic focal mechanism, information on the source mechanisms and in-situ stress status variations can be quantitatively revealed to effectively optimize the reservoir reconstruction design for increasing production. This paper reviews the recent progress in hydraulic fracturing induced microseismic focal mechanism research. We summarize their main principles and provide a detailed introduction of the research advances in source modeling, microseismic data synthesis, and focal mechanism inversion. We also discuss the challenges and limitations in the current microseismic focal mechanism research and propose prospects for future research ideas and directions. 相似文献
5.
水力压裂是页岩气开发过程中的核心增产技术,微地震则广泛用于压裂分析、水驱前缘监测和储层描述.微地震反演过程中,用于反演的速度模型往往基于测井、地震或标定炮资料构建,忽略了压裂过程中裂缝及孔隙流体压力变化对地层速度的影响.本文首先基于物质守恒、渗流理论和断裂力学模拟三维水力压裂过程,得到地下裂缝发育特征和孔隙压力分布.继而根据Coates-Schoenberg方法和裂缝柔量参数计算裂缝和孔隙压力对速度场的影响,得到压裂过程中的实时速度模型.最后利用三维射线追踪方法正演微地震走时和方位信息,并采用常规微地震定位方法反演震源位置及进行误差分析.数值模拟结果表明,检波器空间分布影响定位精度,常规方法的定位误差随射线路径在压裂带中传播距离增加而变大,且不同压裂阶段的多点反演法与单点极化法精度相当. 相似文献
6.
In hydraulic fracturing treatments, locating not only hydraulic fractures but also any pre‐existing natural fractures and faults in a subsurface reservoir is very important. Hydraulic fractures can be tracked by locating microseismic events, but to identify the locations of natural fractures, an additional technique is required. In this paper, we present a method to image pre‐existing fractures and faults near a borehole with virtual reverse vertical seismic profiling data or virtual single‐well profiling data (limited to seismic reflection data) created from microseismic monitoring using seismic interferometry. The virtual source data contain reflections from natural fractures and faults, and these features can be imaged by applying migration to the virtual source data. However, the imaging zone of fractures in the proposed method is strongly dependent on the geographic extent of the microseismic events and the location and direction of the fracture. To verify our method, we produced virtual reverse vertical seismic profiling and single‐well profiling data from synthetic microseismic data and compared them with data from real sources in the same relative position as the virtual sources. The results show that the reflection travel times from the fractures in the virtual source data agree well with travel times in the real‐source data. By applying pre‐stack depth migration to the virtual source data, images of the natural fractures were obtained with accurate locations. However, the migrated section of the single‐well profiling data with both real and virtual sources contained spurious fracture images on the opposite side of the borehole. In the case of virtual single‐well profiling data, we could produce correct migration images of fractures by adopting directional redatuming for which the occurrence region of microseismic events is divided into several subdivisions, and fractures located only on the opposite side of the borehole are imaged for each subdivision. 相似文献
7.
在水力压裂施工中,如何有效获取压裂过程中产生的裂缝形态以及裂缝的动态扩展过程一直是困扰学术界和工业界的问题。目前,常规利用微地震事件定位结果进行分析的方法存在需要人工干预、散点信息表示能力不足等问题;采用数值模拟分析的方法往往因复杂的地下介质情况而引入计算偏差。本文基于非监督学习算法,通过提取微地震事件的空间和时间信息,实现对裂缝平面的识别以及裂缝网络拓展路径的分析;并通过引入水力压裂岩石物理实验,利用实际监测获得的声发射数据以及对应的真实破裂情况的CT扫描数据,检验方法的可行性。最终结果表明,本文所提方法对主断裂有较好的识别效果,识别结果与CT扫描的真实结果吻合性较好。 相似文献
8.
Yang Zhang Leo Eisner William Barker Michael C. Mueller Kevin L. Smith 《Geophysical Prospecting》2013,61(5):919-930
We develop a methodology to obtain a consistent velocity model from calibration shots or microseismicity observed on a buried array. Using a layered 1D isotropic model derived from checkshots as an initial velocity model, we invert P‐wave arrival times to obtain effective anisotropic parameters with a vertical axis of symmetry (VTI). The nonlinear inversion uses iteration between linearized inversion for anisotropic parameters and origin times or depths, which is specific to microseismic monitoring. We apply this technique to multiple microseismic events from several treatments within a buried array. The joint inversion of selected events shows a largely reduced RMS error indicating that we can obtain robust estimates of anisotropic parameters, however we do not show improved source locations. For joint inversion of multiple microseismic events we obtained Thomsen anisotropic parameters ε of 0.15 and δ of 0.05, which are consistent with values observed in active seismic surveys. These values allow us to locate microseismic events from multiple hydraulic fracture treatments separated across thousands of metres with a single velocity model. As a result, we invert the effective anisotropy for the buried array region and are able to provide a more consistent microseismicity mapping for past and future hydraulic fracture stimulations. 相似文献
9.
Borehole fluid injections are accompanied by microseismic activity not only during but also after termination of the fluid injection. Previously, this phenomenon has been analysed, assuming that the main triggering mechanism is governed by a linear pressure diffusion in a hydraulically isotropic medium. In this context the so‐called back front of seismicity has been introduced, which allows to characterize the hydraulic transport from the spatiotemporal distribution of post‐injection induced events. However, rocks are generally anisotropic, and in addition, fluid injections can strongly enhance permeability. In this case, permeability becomes a function of pressure. For such situations, we carry out a comprehensive study about the behaviour and parametrization of the back front. Based on a model of a factorized anisotropic pressure dependence of permeability, we present an approach to reconstruct the principal components of the diffusivity tensor. We apply this approach to real microseismic data and show that the back front characterizes the least hydraulic transport. To investigate the back front of non‐linear pore‐fluid pressure diffusion, we numerically consider a power‐law and an exponential‐dependent diffusivity. To account for a post‐injection enhanced hydraulic state of the rock, we introduce a model of a frozen (i.e., nearly unchanged after the stimulation) medium diffusivity and generate synthetic seismicity. We find that, for a weak non‐linearity and 3D exponential diffusion, the linear diffusion back front is still applicable. This finding is in agreement with microseismic data from Ogachi and Fenton Hill. However, for a strong non‐linear fluid–rock interaction such as hydraulic fracturing, the back front can significantly deviate from a time dependence of a linear diffusion back front. This is demonstrated for a data set from the Horn River Basin. Hence, the behaviour of the back front is a strong indicator of a non‐linear fluid–rock interaction. 相似文献
10.
3D anisotropic waveform inversion could provide high-resolution velocity models and improved event locations for microseismic surveys. Here we extend our previously developed 2D inversion methodology for microseismic borehole data to 3D transversely isotropic media with a vertical symmetry axis. This extension allows us to invert multicomponent data recorded in multiple boreholes and properly account for vertical and lateral heterogeneity. Synthetic examples illustrate the performance of the algorithm for layer-cake and ‘hydraulically fractured’ (i.e. containing anomalies that simulate hydraulic fractures) models. In both cases, waveform inversion is able to reconstruct the areas which are sufficiently illuminated for the employed source-receiver geometry. In addition, we evaluate the sensitivity of the algorithm to errors in the source locations and to band-limited noise in the input displacements. We also present initial inversion results for a microseismic data set acquired during hydraulic fracturing in a shale reservoir. 相似文献
11.
WANG Xiaoqiong XU Jianguo ZHAO Chenxu LIU Tongyuan GE Hongkui SHEN Yinghao WU Shan YU Jiayao HUANG Rongyan 《中国地震研究》2020,34(2):167-186
The premise of hydraulic fracturing is to have an accurate and detailed understanding of the rock mechanical properties and fracture propagation law of shale reservoirs. In this paper, a comprehensive evaluation of the mechanical properties of the shale oil reservoir in the south of Songliao Basin is carried out. Based on the experiments and the in-situ stress analysis, the fracture propagation law of three types of shale reservoirs is obtained, and the suggestions for fracturing are put forward. The results have shown that the fracture propagation of pure shale and low mature reservoir is easy to open along the bedding plane under compression loading, which is greatly influenced by the bedding. Sand-bearing shale is slightly better, the fractures of which are not easy to open along the bedding plane. The mechanical experimental results show that all the samples have the characteristics of low compressive strength, low Young''s modulus and strong anisotropy, indicating that the shale oil reservoir is certain plastic, which is related to its high clay mineral content and controlled by the bedding development. Compared with pure shale and low mature shale, the sand-bearing shale has less clay content and less developed bedding, which maybe the main reason for its slightly better brittleness. Overall, the expansion of hydraulic fracture is controlled by in-situ stress and bedding. Because of the development of bedding, it is easy to form horizontal fractures. Thus it is not suitable for horizontal well fracturing. Because of the high content of clay minerals, the applicability of conventional slick hydraulic fracturing fluid is poor. It is suggested to use vertical well or directional well to carry out volume fracturing. In this way, the effect of bedding can be effectively used to open and connect the bedding and form a larger fracture network. 相似文献
12.
Hydraulic stimulation and geothermal reservoir operation may compromise the rock mechanical stability and trigger microseismic events. The mechanisms leading to this induced seismicity are still not completely understood. It is clear that injection causes an overpressure that reduces the effective stress, bringing the system closer to failure conditions. However, rock instability may not result only from hydraulic effects, but also from thermal effects. In fact, hydro-mechanical (i.e., isothermal) models often fail to reproduce field observations because the injection of cold water into a hot reservoir induces thermal stresses due to rock contraction. Thus, rock instability is likely to result from the superposition of hydraulic and thermal effects. Here, we perform coupled thermo-hydro-mechanical and hydro-mechanical simulations to investigate the effects of cold water injection in a fracture zone-intact rock system. Results show that thermal effects induce a significant perturbation on the stress in the intact rock affected by the temperature drop. This perturbation is likely to trigger induced seismicity in the surroundings of critically oriented fractures near the injection well. Hydro-mechanical simulations show that the behavior depends on the orientation of the faults and on the initial stress tensor. In the direction of the fractures, where the strains are more constrained, total stress increases with increasing pressure; thus, deviatoric stress increases or decreases depending on the initial stress state. The comparison between hydraulic and thermal effects shows that, when the largest confining stress acts perpendicular to the fractures, thermoelastic effects dominate and could trigger induced seismicity. 相似文献
13.
The hydrocarbon industry is moving increasingly towards tight sandstone and shale gas resources – reservoirs that require fractures to be produced economically. Therefore, techniques that can identify sets of aligned fractures are becoming more important. Fracture identification is also important in the areas of coal bed methane production, carbon capture and storage (CCS), geothermal energy, nuclear waste storage and mining. In all these settings, stress and pore pressure changes induced by engineering activity can generate or reactivate faults and fractures. P‐ and S‐waves are emitted by such microseismic events, which can be recorded on downhole geophones. The presence of aligned fracture sets generates seismic anisotropy, which can be identified by measuring the splitting of the S‐waves emitted by microseismic events. The raypaths of the S‐waves will have an arbitrary orientation, controlled by the event and geophone locations, meaning that the anisotropy system may only be partly illuminated by the available arrivals. Therefore to reliably interpret such splitting measurements it is necessary to construct models that compare splitting observations with modelled values, allowing the best fitting rock physics parameters to be determined. Commonly, splitting measurements are inverted for one fracture set and rock fabrics with a vertical axis of symmetry. In this paper we address the challenge of identifying multiple aligned fracture sets using splitting measured on microseismic events. We analyse data from the Weyburn CCS‐EOR reservoir, which is known to have multiple fracture sets, and from a hydraulic fracture stimulation, where it is believed that only one set is present. We make splitting measurements on microseismic data recorded on downhole geophone arrays. Our inversion technique successfully discriminates between the single and multiple fracture cases and in all cases accurately identifies the strikes of fracture sets previously imaged using independent methods (borehole image logs, core samples, microseismic event locations). We also generate a synthetic example to highlight the pitfalls that can be encountered if it is assumed that only one fracture set is present when splitting data are interpreted, when in fact more than one fracture set is contributing to the anisotropy. 相似文献
14.
In June/July 2000, a hydraulic stimulation experiment took place at the geothermal EGS site of Soultz-sous-Forêts (Alsace,
France) in order to enhance the permeability of the fractured granitic massif at 5 km depth. As it is well known that fluid
injections tend to induce microseismic events, a downhole and a surface seismological network have been installed to monitor
the seismic activity during the stimulation test. 23400 m3 of fluid have been injected in the rock volume through the open-hole section (4400 m–5000 m) of the well GPK2 at increasing
rates of 30 l.s−1, 40 l.s−1 and then 50 l.s−1. More than 7200 microseismic events in the magnitude range –0.9 to 2.6 have been precisely located through a simultaneous
inversion of the seismic velocity structure and location parameters. The analysis of the behavior of the seismicity relative
to the hydraulic parameters gives important information about the geothermal reservoir. It appears that the evolution of the
seismicity strongly depends on the variations of the injection rate: An increase or a decrease leads to changes of the velocity
structure, the number and magnitude of microseismic events. This involves different hydro-mechanical processes between the
fluid flow and the fracture planes, which will control the final shape of the microseismic cloud. Moreover, the study of the
variations of the b-value with time suggests that the stimulation experiment produces a large proportion of small earthquakes, but records of
events of magnitude higher than 2 indicate that fluid injection could reactivate structures whose dimensions allow the generation
of such earthquakes. 相似文献
15.
Rui-Zhao Yang Zheng-Guang Zhao Wei-Jun Peng Yu-Bo Gu Zhan-Gang Wang Xi-Qin Zhuang 《应用地球物理》2013,10(2):157-169
The development of unconventional resources, such as shale gas and tight sand gas, requires the integration of multi-disciplinary knowledge to resolve many engineering problems in order to achieve economic production levels. The reservoir heterogeneity revealed by different data sets, such as 3D seismic and microseismic data, can more fully reflect the reservoir properties and is helpful to optimize the drilling and completion programs. First, we predict the local stress direction and open or close status of the natural fractures in tight sand reservoirs based on seismic curvature, an attribute that reveals reservoir heterogeneity and geomechanical properties. Meanwhile, the reservoir fracture network is predicted using an ant-tracking cube and the potential fracture barriers which can affect hydraulic fracture propagation are predicted by integrating the seismic curvature attribute and ant-tracking cube. Second, we use this information, derived from 3D seismic data, to assist in designing the fracture program and adjusting stimulation parameters. Finally, we interpret the reason why sand plugs will occur during the stimulation process by the integration of 3D seismic interpretation and microseismic imaging results, which further explain the hydraulic fracture propagation controlling factors and open or closed state of natural fractures in tight sand reservoirs. 相似文献
16.
微地震事件的空间分布可以用来监测水力压裂过程中裂缝的发育情况。因此,震源定位是微震监测中重要的环节。震源定位依赖准确的速度模型,而震源位置和速度模型的耦合易导致线性迭代的同时反演方法陷入局部极小值。邻近算法作为一种非线性全局优化算法,能够最大程度地避免陷入局部最优解。本文将邻近算法应用于单井监测的微震定位和一维速度模型同时反演,首先利用邻近算法搜索一维速度模型,再使用网格搜索方法进行震源定位,并根据定位的走时残差产生新的速度模型,最后通过若干次迭代使其收敛到最优解。理论和实际数据结果均表明该方法能够避免局部最优解,得到较为可靠的震源位置和一维速度模型。 相似文献
17.
Microseismicity variations in space and time due to pore pressure changes in an inhomogeneous permeable sample are considered. The analysis is based on laboratory experiments performed to study the relation between acoustic emission (which corresponds to microseismic events in real scale) and pore pressure changes due to water injection into, or release from, a porous sample under load. The study showed not only some characteristics of microseismicity variations during non-stationary fluid flow but also the possibility to solve the inverse problem of estimating local permeability from variations of microseismic activity in a particular volume of the porous medium. 相似文献
18.
During seismic monitoring of hydraulic fracturing treatment, it is very common to ignore the deviations of the monitoring or treatment wells from their assumed positions. For example, a well is assumed to be perfectly vertical, but in fact, it deviates from verticality. This can lead to significant errors in the observed azimuth and other parameters of the monitored fracture‐system geometry derived from microseismic event locations. For common hydraulic fracturing geometries, a 2° deviation uncertainty on the positions of the monitoring or treatment well survey can cause a more than 20° uncertainty of the inverted fracture azimuths. Furthermore, if the positions of both the injection point and the receiver array are not known accurately and the velocity model is adjusted to locate perforations on the assumed positions, several‐millisecond discrepancies between measured and modeled SH‐P traveltime differences may appear along the receiver array. These traveltime discrepancies may then be misinterpreted as an effect of anisotropy, and the use of such anisotropic model may lead to the mislocation of the detected fracture system. The uncertainty of the relative positions between the monitoring and treatment wells can have a cumulative, nonlinear effect on inverted fracture parameters. We show that incorporation of borehole deviation surveys allows reasonably accurate positioning of the microseismic events. In this study, we concentrate on the effects of horizontal uncertainties of receiver and perforation positions. Understanding them is sufficient for treatment of vertical wells, and also necessary for horizontal wells. 相似文献
19.
Large sand intrusions often exhibit conical morphologies analogous to magmatic intrusions such as saucer-shaped or cup-shaped sills. Whereas some physical processes may be similar, we show with scaled experiments that the formation of conical sand intrusions may be favoured by the pore-pressure gradients prevailing in the host rock before sand injection. Our experiments involve injecting air into a permeable and cohesive analogue material to produce hydraulic fractures while controlling the pore pressure field. We control the state of overpressure in the overburden by applying homogeneous basal pore pressure, and then adding a second local pore pressure field by injecting air via a central injector to initiate hydraulic fractures near the injection point. In experiments involving small vertical effective stresses (small overburden, or high pore fluid overpressure), the fracturing pressure (λfract) is supralithostatic and two dipping fractures are initiated at the injection point forming a conical structure. From theoretical considerations, we predict that high values of λfract are due to strong cohesion or high pore fluid overpressure distributed in the overburden. Such conditions are favoured by the pore pressure/stress coupling induced by both pore pressure fields. The dips of cones can be accounted for elastic-stress rotation occurring around the source. Contrary to magmatic chamber models, the aqueous fluid overpressure developed in a parent sandbody (and prevailing before the formation of injectites) may diffuse into the surrounding overburden, thus favouring stress rotation and the formation of inclined sheets far from the parent source. For experiments involving higher vertical effective stresses (thick overburden or low pore fluid overpressure), the fracturing pressure is lower than the lithostatic stress, and a single fracture is opened in mode I which then grows vertically. At a critical depth, the fracture separates into two dilatant branches forming a flat cone. We make use of a P.I.V. (Particle Imaging Velocimetry) technique to analyse plastic deformation, showing that these inclined fractures are opened in mixed modes. Close to the surface, they change into steep shear bands where fluids can infiltrate. The final morphology of the fracture network is very similar to the common tripartite architecture of various injection complexes, indicating that different mechanisms may be involved in the formation of dykes. Feeder dykes under the sill zones may open as tensile fractures, while overlying dykes may be guided by the deformation induced by the growth of sills. These deformation conditions may also favour the formation of fluid escape structures and pockmarks. 相似文献
20.
G. Michael Hoversten Michael Commer Eldad Haber Christoph Schwarzbach 《Geophysical Prospecting》2015,63(6):1508-1526
As motivation for considering new electromagnetic techniques for hydraulic fracture monitoring, we develop a simple financial model for the net present value offered by geophysical characterization to reduce the error in stimulated reservoir volume calculations. This model shows that even a 5% improvement in stimulated reservoir volume for a 1 billion barrel (bbl.) field results in over 1 billion U.S. dollars (US$) in net present value over 24 years for US$100/bbl. oil and US$0.5 billion for US$50/bbl. oil. The application of conductivity upscaling, often used in electromagnetic modeling to reduce mesh size and thus simulation runtimes, is shown to be inaccurate for the high electrical contrasts needed to represent steel‐cased wells in the earth. Fine‐scale finite‐difference modeling with 12.22‐mm cells to capture the steel casing and fractures shows that the steel casing provides a direct current pathway to a created fracture that significantly enhances the response compared with neglecting the steel casing. We consider conductively enhanced proppant, such as coke‐breeze‐coated sand, and a highly saline brine solution to produce electrically conductive fractures. For a relatively small frac job at a depth of 3 km, involving 5,000 bbl. of slurry and a source midpoint to receiver separation of 50 m, the models show that the conductively enhanced proppant produces a 15% increase in the electric field strength (in‐line with the transmitter) in a 10‐Ωm background. In a 100‐Ωm background, the response due to the proppant increases to 213%. Replacing the conductive proppant by brine with a concentration of 100,000‐ppm NaCl, the field strength is increased by 23% in the 100‐Ωm background and by 2.3% in the 10‐Ωm background. All but the 100,000‐ppm NaCl brine in a 10‐Ωm background produce calculated fracture‐induced electric field increases that are significantly above 2%, a value that has been demonstrated to be observable in field measurements. 相似文献