共查询到20条相似文献，搜索用时 15 毫秒
1.
A 3D relocation technique permits precise locations of induced earthquakes. Geostatistical processing using the data of 87 boreholes provides the basis of a precise 3D structure, with a dome geometry. Conventional laboratory mechanical tests performed on deep rock samples (1000 m to 5000 m) define the rock properties at depths similar to those of the seismic events (1<M
L<4.2) that range from 1 to 7 km.In the studied period, most (85%), of the events were located above the gas reservoir, with very few located in the reservoir itself. Because the production parameters (50 MPa depletion of the gas pressure reservoir) are homogeneous throughout the gas field, the lateral inhomogeneity of the seismic rupture locations are a consequence of variations in the rheological response of the dome to the deformation induced by gas production.Here a ratio of two is found between the elastic modulus of the seismic rock matrix and the elastic modulus of the aseismic rock matrix. The contrast in strength is at least as great, if not greater. Repeated measured surface deformations involve the whole structure. Spatial and temporal deformations indicate that aseismic deformation is quantitatively the main process of this structural deformation. The heterogeneous stress pattern inferred fromP-axes of induced earthquakes disagrees with the tectonic regional stress field. The radial distribution ofP-axes towards the gas reservoir probably reflects the production induced deformation. The inferred deformation of the dome occurs in response to weak induced stresses. 相似文献
2.
Yipeng Zhang Mark Person John Rupp Kevin Ellett Michael A. Celia Carl W. Gable Brenda Bowen James Evans Karl Bandilla Peter Mozley Thomas Dewers Thomas Elliot 《Ground water》2013,51(4):525-538
A series of Mb 3.8–5.5 induced seismic events in the midcontinent region, United States, resulted from injection of fluid either into a basal sedimentary reservoir with no underlying confining unit or directly into the underlying crystalline basement complex. The earthquakes probably occurred along faults that were likely critically stressed within the crystalline basement. These faults were located at a considerable distance (up to 10 km) from the injection wells and head increases at the hypocenters were likely relatively small (～70–150 m). We present a suite of simulations that use a simple hydrogeologic‐geomechanical model to assess what hydrogeologic conditions promote or deter induced seismic events within the crystalline basement across the midcontinent. The presence of a confining unit beneath the injection reservoir horizon had the single largest effect in preventing induced seismicity within the underlying crystalline basement. For a crystalline basement having a permeability of 2 × 10?17 m2 and specific storage coefficient of 10?7/m, injection at a rate of 5455 m3/d into the basal aquifer with no underlying basal seal over 10 years resulted in probable brittle failure to depths of about 0.6 km below the injection reservoir. Including a permeable (kz = 10?13 m2) Precambrian normal fault, located 20 m from the injection well, increased the depth of the failure region below the reservoir to 3 km. For a large permeability contrast between a Precambrian thrust fault (10?12 m2) and the surrounding crystalline basement (10?18 m2), the failure region can extend laterally 10 km away from the injection well. 相似文献
3.
D.A. Angus J.‐M. Kendall Q.J. Fisher J.M. Segura S. Skachkov A.J.L. Crook M. Dutko 《Geophysical Prospecting》2010,58(5):901-914
In this paper, we investigate production induced microseismicity based on modelling material failure from coupled fluid‐flow and geomechanical simulation. The model is a graben style reservoir characterized by two normal faults subdividing a sandstone reservoir into three compartments. The results are analysed in terms of spatial and temporal variations in distribution of material failure. We observe that material failure and hence potentially microseismicity is sensitive to not only fault movement but also fluid movement across faults. For sealing faults, failure is confined to the volume in and around the well compartment, with shear failure localized along the boundaries of the compartment and shear‐enhanced compaction failure widespread throughout the reservoir compartment. For non‐sealing faults, failure is observed within and surrounding all three reservoir compartments as well as a significant distribution located near the surface of the overburden. All shear‐enhanced compaction failures are localized within the reservoir compartments. Fault movement leads to an increase in shear‐enhanced compaction events within the reservoir as well as shear events located within the side‐burden adjacent to the fault. We also evaluate the associated moment tensor mechanisms to estimate the pseudo scalar seismic moment of failure based on the assumption that failure is not aseismic. The shear‐enhanced compaction events display a relatively normal and tight pseudo scalar seismic moment distribution centred about 106 Pa, whereas the shear events have pseudo scalar seismic moments that vary over three orders of magnitude. Overall, the results from the study indicate that it may be possible to identify compartment boundaries based on the results of microseismic monitoring. 相似文献
4.
随着二氧化碳地质封存、深部地热开采、地下储气库建设、页岩气开发、二次驱油/驱气等工业应用的快速发展,与地下流体注入有关的诱发地震活动呈现一定的增加趋势.利用声发射实验观测油气田典型岩石在三轴压缩条件下变形破坏过程与声发射活动特征,对研究注水诱发地震过程有着重要意义.本文利用四川盆地现场采集的震旦系白云岩及页岩,采用实验室声发射技术观测研究岩石三轴压缩变形破坏过程中地震波速度等物性参数及声发射事件时空分布特征.实验结果表明:震旦系白云岩及页岩在变形破坏过程中均有一定的声发射活动.根据声发射定位结果,声发射主要集中在破坏前后的较短时间内,页岩的层理面为结构弱面,控制最终破坏面的形态及声发射特征.根据应力-应变结果,白云岩在压缩的后期阶段有一定的扩容现象,但页岩在整个压缩阶段均没有明显扩容现象.研究结果表明四川盆地较古老的白云岩及页岩具有脆性破坏特征,地下流体注入容易诱发微震活动,形成裂缝,有利于页岩气压裂开采.微震活动有利于监测裂缝的发生发展,但同时在页岩气开发及二氧化碳地质封存时应采取相应的预防控制措施进行安全合理的储盖层管理,避免灾害性诱发地震的发生. 相似文献
5.
水库诱发地震与构造地震震源参数特征差异性研究——以龙滩水库为例 总被引：8，自引：0，他引：8
首先讨论龙滩库区水库蓄水与地震活动之间的关系,发现龙滩水库诱发地震特征十分明显,地震共分5丛呈丛集分布.利用库区架设的24个固定和流动台站记录的数字记录资料,在研究得到龙滩库区非弹性衰减和台站场地响应的基础上,精确测定得到了该地区总共1616个ML≥0.1级地震的震源参数,比较了水库诱发地震与构造地震震源参数特征的差异,得到了以下主要结论:1)龙滩水库地震活动与水库蓄水关系密切,不同蓄水阶段5丛的地震活动状态不同,局部断裂构造发育以及岩石透水性能影响着地震活动对蓄水过程的响应.2)龙滩水库诱发地震的地震矩M0随震级ML的增大而增大,两者之间存在较好的线性关系,统计关系为LogM0=1.07 ML+10.17.应力降与地震大小之间的关系和Nuttli的板内地震为增加应力降(ISD)模型的结果比较吻合,统计关系为LogΔσ=0.71 ML-2.89.3)龙滩水库地区地震辐射能量和地震视应力均随震级的增大而增大,后者意味着大地震是比小地震更高效率的地震能量辐射体.4)总体上不同丛地震应力降水平存在差异.地震应力降空间分布上与库水深度有较好的一致性,即库水深的区域应力降水平高.5)与同震级的构造地震相比,水库诱发地震的应力降值比前者明显偏低,大约小10倍.这可能是由于水库蓄水造成地下介质孔隙压力增大或者水的润滑作用,从而导致在一个比较低的构造应力情况下发生水库诱发地震. 相似文献
6.
7.
《地震工程与结构动力学》2018,47(6):1544-1565
This study evaluates capabilities of the VecTor4 computer program to assess the seismic performance of a concrete containment structure subjected to both design‐basis and beyond design‐basis earthquakes. Based on a detailed nonlinear finite element model of the nuclear power plant containment structure, the seismic performance in terms of stress and strain distributions, cracks, yielding of reinforcement bars and tendons, as well as overall failure mechanism, is thoroughly reviewed. In addition, the post‐seismic performance of the containment structure subjected to a subsequent pressure increase is also investigated. Consideration is also given to the time‐dependent parameters of creep, shrinkage, and relaxation of prestressing tendons. It is found that the time‐dependent parameters and earthquake have a non‐negligible impact on pressure‐induced structural performance but that this impact can be predicted for a given earthquake. The proposed method is useful for assessing the progressive failure behaviour of containment structure under chain events. 相似文献
8.
糯扎渡水库地区地震活动和震源参数研究 总被引：1，自引：0，他引：1
本文利用糯扎渡水库台网和景洪水库台网记录的地震波形反演震源谱, 计算了2011年1月至2014年3月糯扎渡水库附近143个1.0级以上地震的震源参数。 分析研究发现: ① 糯扎渡水库蓄水后, 库区内地震活动明显增强, 尤其是在水位上升3~4个月后; ② 水库开始蓄水后, 库区内地震震源深度变浅, 一段时间以后则与库区外的差别减小; ③ 地震拐角频率随地震矩增大而减小, 且蓄水以后, 库区内地震拐角频率的对数与地震矩的对数线性关系更明显; ④ 应力降、 视应力与地震矩存在正相关关系, 且相同地震矩时库区内的应力降与视应力低于蓄水前和库区外的地震应力降与视应力值; ⑤ 蓄水对库区内地震的影响可延伸至库底10 km深度处, 且在3~6 km深度内影响最大、 库区内外地震频度与视应力均值差异最显著。 相似文献
9.
10.
From August 1976 to June 1979, more than 220 seismic events were recorded by a 7-station network operating over an active gas field in France. 120 earthquakes with magnitudes greater than 1.0 were selected for a detailed study in order to understand the seismic response to a decrease in the gas reservoir pressure of 50 MPa. Composite focal mechanisms were determined for these events. We found that where preexisting fractures were known from well-logs, faults with the most favorable orientation were reactivated. In the part of the gas field where no preexisting faults were known, orientation of the fracture planes was close to the orientations of stratigraphic interfaces.Our kinematic model, constrained by the distribution of the earthquakes and their focal mechanisms, consists of the downward dropping of a block located close to the top of the gas field. The main faults, as defined by the spatial density of the seismic activity, define the limits of this block. The displacements observed along these faults can be either normal or reverse mechanism, dependent on the preexisting orientation of each fault. Another important result is the evidence that the depletion of the gas reservoir exerts the major influence on local seismic displacements, obscuring the effect of regional setting. 相似文献
11.
During the 2000 activity of Miyake-jima volcano, Japan, we detected long period seismic signals with initial pulse widths of 1-2 s, accompanied by infrasonic pulses with almost the same pulse widths. The seismic signals were observed from 13 July 2000, a day before the second summit eruption. The occurrences of the seismic signals were intermittent with a gradual increase in their magnitudes and numbers building toward a significant explosive eruption on 18 August. After the eruption, the seismic and infrasonic events ceased. The results of a waveform inversion show that the initial motions were excited by an isotropic inflation source beneath the south edge of the caldera at a depth of 1.4 km. On the other hand, the sources of the infrasonic pulses were located in the summit caldera area. The times at which the infrasonic pulses were emitted at the surface were delayed by about 3 s from the origin times of the seismic events. It is suggested that small isotropic inflations excited seismic waves in the crust and simultaneously caused acoustic waves that traveled in the conduit and produced infrasonic pulses at the crater bottom. Considering the observed time differences and gas temperatures emitted from the vent, the conduit should have been filled with vapor mixed with SO2 gas and volcanic ash. The change of the time differences between the seismic and infrasonic signals suggests that the seismic source became shallower within half a day before the August 18 explosive eruption. We interpret the source process as a fragmentation process of magma in which gas bubbles burst and quickly released part of the pressure that had been sustained by the tensional strength of magma. 相似文献
12.
Kusnowidjaja Megawati Tso-Chien Pan Kazuki Koketsu 《Soil Dynamics and Earthquake Engineering》2005,25(1):11-25
Singapore and Kuala Lumpur, the capital of Malaysia, may well represent the classic examples of area with low seismic hazard but with high consequence. Both cities are located in a low-seismicity region of Southeast Asia, where active seismic sources are located more than 300 km away. Seismic designs have not been implemented in this seemingly low-hazard region though distant earthquakes in Sumatra had frequently shaken high-rise structures in the two cities. Several studies have been conducted to systematically assess the seismic hazards of Singapore and the Malay Peninsula. The present research particularly addresses issues in deriving a new set of attenuation relationships of peak ground acceleration (PGA), peak ground velocity (PGV) and response spectral acceleration (RSA) for the Sumatran-subduction earthquakes. To be relevant for the seismic hazard assessment of the remote metropolises, the derived attenuation relationships cover a long distance range from 150 to 1500 km. The attenuation relationships are derived using synthetic seismograms that account for source and path effects. The uncertainties in rupture parameters, such as stress drop, strike, dip and rake angles, have been defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. The seismic potential of the Sumatran subduction zone are high in the region from 2°N to 5°S as there has been no recurrence of great thrust events since 1861. A large event with Mw greater than 7.8 in this particular subduction zone may be capable of generating destructive ground motions in Singapore and Kuala Lumpur, even at a distance of 700 km. 相似文献
13.
Real-time envelope cross-correlation detector: application to induced seismicity in the Insheim and Landau deep geothermal reservoirs 总被引：1，自引：0，他引：1
Margarete Vasterling Ulrich Wegler Jan Becker Andrea Brüstle Monika Bischoff 《Journal of Seismology》2017,21(1):193-208
We develop and test a real-time envelope cross-correlation detector for use in seismic response plans to mitigate hazard of induced seismicity. The incoming seismological data are cross-correlated in real-time with a set of previously recorded master events. For robustness against small changes in the earthquake source locations or in the focal mechanisms we cross-correlate the envelopes of the seismograms rather than the seismograms themselves. Two sequenced detection conditions are implemented: After passing a single trace cross-correlation condition, a network cross-correlation is calculated taking amplitude ratios between stations into account. Besides detecting the earthquake and assigning it to the respective reservoir, real-time magnitudes are important for seismic response plans. We estimate the magnitudes of induced microseismicity using the relative amplitudes between master event and detected event. The real-time detector is implemented as a SeisComP3 module. We carry out offline and online performance tests using seismic monitoring data of the Insheim and Landau geothermal power plants (Upper Rhine Graben, Germany), also including blasts from a nearby quarry. The comparison of the automatic real-time catalogue with a manually processed catalogue shows, that with the implemented parameters events are always correctly assigned to the respective reservoir (4 km distance between reservoirs) or the quarry (8 km and 10 km distance, respectively, from the reservoirs). The real-time catalogue achieves a magnitude of completeness around 0.0. Four per cent of the events assigned to the Insheim reservoir and zero per cent of the Landau events are misdetections. All wrong detections are local tectonic events, whereas none are caused by seismic noise. 相似文献
14.
为了解水库蓄水过程中,水库诱发地震活动的动态响应机制,本文建立了针对水库诱发地震(Reservoir-induced seismicity,RIS)定量化研究的数理模型,并以紫坪铺水库为例,对库区地质构造及水文地质结构条件、水库蓄水后库区小震活动时空演化特征进行了详细的研究.在此基础上,利用有限元方法计算了水库蓄水过程中弹性附加应力场、有效附加应力场、孔隙压力和断层稳定性的动态变化,讨论了RIS时空演化与库水加卸载及渗透过程的动态响应关系.结果表明:(1) RIS诱发机制的定量化模型可分为2个层次:一是以孔隙介质为载体的流体渗流对岩体变形和稳定性的影响,由流-固耦合形式的岩体变形与孔隙渗流模型进行描述;二是对断层相关的RIS定量研究可将水库附加水头压力沿断层面(区)的扩散与断层库仑应力变化联系起来.两种形式模型方法的结合能为RIS定量研究提供一个相对宏观的力学框架;(2) 断裂渗透结构对孔隙压力变化下断裂的力学响应具有重要的影响,研究区主干断裂可能属于一种上盘破碎带导水、下盘地层及断层核阻水的"下阻上导型"的渗透结构类型,不同程度的具有使地表水体向深部渗流的通道性.库区深部岩体渗透稳定性的差异在很大程度上导致了诱发地震活动对岩性条件的依赖.(3) 紫坪铺水库蓄水后,小震活动在空间分布上呈现出条带状分布、丛集分布和地震迁移的特点,小震震源深度优势分布在地下4~10 km范围内,在通济场断裂与安县-灌县断裂的深部汇聚区域震源分布最为密集.同时,小震活动主要集中发生在脆性程度高、渗透稳定性低的碳酸盐岩地层中,而在岩性较软弱、渗透稳定性高的三叠系须家河组砂泥岩和煤系地层中很少有地震发生.在水库蓄水后地震活动的时间响应特征上,水库西南侧和东北侧两个丛集区的小震活动可能属于"快速响应型"RIS,而都江堰小震群活动可能属于"滞后响应型"诱发地震活动;(4) RIS的发生与库水加卸载及渗透过程中库底岩体有效应力的变化密切相关.在以挤压为主的构造应力环境中,库体荷载作用的结果一般会使库底断层更趋向稳定,而水库附加水头压力扩散的效应则是促使断层趋向失稳,正是这个矛盾双方相互制约与平衡的动态过程,控制了断层库仑应力变化的取向,从而决定了RIS时空演化的规律. 相似文献
15.
Qian Wang Xinxin Yin Changsheng Jiang Cong Jiang Yan Zhang Hongyu Zhai Yanbao Zhang Guijuan Lai Fengling Yin 《地震科学（英文版）》2021,34(3):286-298
In the new types of industrial activities including unconventional energy extraction associated with shale gas and hot dry rock, gas reservoir operations, CO2 geological storage, undergoing research on induced earthquake forecasting has become one of the forward positions of current seismology. As for the intense actual demand, the immature research on induced earthquake forecasting has already been applied in pre-assessment of site safety and seismic hazard and risk management. This work will review systematically recent advances in earthquake forecasting induced by hydraulic fracturing during industrial production from four aspects: earthquake occurrence probability, maximum expected magnitude forecasting, seismic risk analysis for engineering and social applications and key scientific problems. In terms of earthquake occurrence probability, we introduce statistical forecasting models such as an improved ETAS and non-stationary ETAS and physical forecasting models such as Seismogenic Index (SI) and hydro-mechanism nucleation. Research on maximum expected magnitude forecasting has experienced four stages of linear relationship with net injection volume of fluid, power exponential relationship and physical forecasting regarding fault parameters. For seismic risk analysis, we focus on probabilistic seismic hazard assessment and quantitative geological susceptibility model. Furthermore, this review is extended to key scientific problems that contain obtaining accurate fault scale and environmental stress state of reservoir, critical physical process of runaway rupture, complex mechanism of fault activation as well as physical mechanism and modeling of trailing effect. This work in understanding induced earthquake forecasting may contribute to unconventional energy development and production, seismic hazard mitigation, emergency management and scientific research as a reference. 相似文献
16.
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. 相似文献
17.
P-wave velocity structure beneath reservoirs and surrounding areas in the lower Jinsha River 
下载免费PDF全文
下载免费PDF全文 Changzai Wang Jianping Wu Lihua Fang Yaning Liu Jing Liu Yan Cai Poren Li 《地震科学（英文版）》2023,36(1):64-75
The lower reaches of the Jinsha River are rich in hydropower resources because of the high mountains, deep valleys, and swift currents in this area. This region also features complex tectonic structures and frequent earthquakes. After the impoundment of the reservoirs, seismic activity increased significantly. Therefore, it is necessary to study the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds, thus providing seismological support for subsequent earthquake prevention and disaster reduction work in reservoir areas. In this study, we selected the data of 7,670 seismic events recorded by the seismic networks in Sichuan, Yunnan, and Chongqing and the temporary seismic arrays deployed nearby. We then applied the double-difference tomography method to this data, to obtain the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds. The results showed that the Jinsha River basin has a complex lateral P-wave velocity structure. Seismic events are mainly distributed in the transition zones between high- and low-velocity anomalies, and seismic events are particularly intense in the Xiluodu and Baihetan reservoir areas. Vertical cross-sections through the Xiangjiaba and Xiluodu reservoir areas revealed an apparent high-velocity anomaly at approximately 6 km depth; this high-velocity anomaly plays a role in stress accumulation, with few earthquakes distributed inside the high-velocity body. After the impoundment of the Baihetan reservoir, the number of earthquakes in the reservoir area increased significantly. The seismic events in the reservoir area north of 27° N were related to the enhanced activity of nearby faults after impoundment; the earthquakes in the reservoir area south of 27° N were probably induced by additional loads (or regional stress changes), and the multiple microseismic events may have been caused by rock rupture near the main faults under high pore pressure. 相似文献
18.
Sawasdee Yordkayhun Alexandra Ivanova Rüdiger Giese Christopher Juhlin Calin Cosma 《Geophysical Prospecting》2009,57(1):125-139
In 2004 three seismic surface sources (VIBSIST, accelerated weight drop and MiniVib) were tested in a pilot study at the Ketzin test site, Germany, a study site for geological storage of CO2 (EU project CO2SINK). The main objectives of this pilot study were to 1) evaluate the response of the Ketzin site to reflection seismics, especially at the planned injection depth, 2) test different acquisition parameters and 3) use the results to guide the planning of the 3D survey. As part of these objectives, we emphasize the source performance comparison in this study. The sources were tested along two perpendicular lines of 2.4 km length each. Data were acquired by shooting at all stations (source and receiver spacing of 20 m) on both lines, allowing common‐midpoint stacked sections to be produced. The sources' signal characteristics based on signal‐to‐noise ratio, signal penetration and frequency content of raw shot records were analysed and stacked sections were compared. The results show that all three surface sources are suitable for reflection seismic studies down to a depth of about 1 km and provide enough bandwidth for resolving the geological targets at the site, i.e., the Weser and Stuttgart Formations. Near surface conditions, especially a thick weathering layer present in this particular area, strongly influence the data quality, as indicated by the difference in reflectivity and signal‐to‐noise ratio of the two common‐midpoint lines. The stacked sections of the MiniVib source show the highest frequency signals down to about 500 ms traveltime (approximately 500 m depth) but also the shallowest signal penetration depth. The VIBSIST source generates signals with the highest signal‐to‐noise ratio and greatest signal penetration depth of the tested sources. In particular, reflections below 900 ms (approximately 1 km depth) are best imaged by the VIBSIST source. The weight drop performance lies in between these two sources and might be recommended as an appropriate source for a 3D survey at this site because of the shorter production time compared to the VIBSIST and MiniVib sources. 相似文献
19.
本文系统地介绍了地震区三角网平差和利用平差结果反演断层参数的随机尝试————动态平差方法.通过新丰江水库5.3级地震区三角测量资料的处理, 进一步阐述如何应用这种方法, 并给出反演5.3级地震断层参数结果:断层走向北东62, 倾向南东, 倾角59, 断层长4.7公里, 上界面深0.45公里, 下界面深4.5公里, 倾向滑距73.4毫米, 走向滑距25.5毫米, 为正断右旋走滑断层, 地震矩 m0=2.01024达因厘米, 应力降=12巴, 释放的应变能下限 E=1.21019尔格.文中依据所得结果, 结合地震机制和构造资料讨论了新丰江水库地震的力源问题, 指出区域应力场是主要力源, 水压应力场是主要的诱发力源。 相似文献
20.
The present work describes the results of a modeling study addressing the geological sequestration of carbon dioxide (CO2) in an offshore multi-compartment reservoir located in Italy. The study is part of a large scale project aimed at implementing carbon capture and storage (CCS) technology in a power plant in Italy within the framework of the European Energy Programme for Recovery (EEPR). The processes modeled include multiphase flow and geomechanical effects occurring in the storage formation and the sealing layers, along with near wellbore effects, fault/thrust reactivation and land surface stability, for a CO2 injection rate of 1 × 106 ton/a. Based on an accurate reproduction of the three-dimensional geological setting of the selected structure, two scenarios are discussed depending on a different distribution of the petrophysical properties of the formation used for injection, namely porosity and permeability. The numerical results help clarify the importance of: (i) facies models at the reservoir scale, properly conditioned on wellbore logs, in assessing the CO2 storage capacity; (ii) coupled wellbore-reservoir flow in allocating injection fluxes among permeable levels; and (iii) geomechanical processes, especially shear failure, in constraining the sustainable pressure buildup of a faulted reservoir. 相似文献