共查询到20条相似文献，搜索用时 46 毫秒
1.
2.
High-resolution seismic surveys were performed on an underground gas storage of Gaz de France at Gournay-sur-Aronde (Oise) in order to obtain information about gas bubble boundaries within the reservoir. For that purpose, a light seismic source Soursile and spreads with small spacing between geophone groups (10 m) and small geophone group extension (10 m or 20 m) were used. In December 1982, a line recorded on the top of the anticline structure of the gas reservoir provided an estimation of the lateral extension of the bubble based on an amplitude anomaly (dim spot) associated with reservoir seismic horizons. In order to estimate the displacement of the gas bubble within the reservoir, three lines were recorded on the southeastern part of the anticline at different time periods (April 1981 and January 1982) with different volumes of gas in the reservoir. It is assumed that a variation of gas volume introduces a seismic velocity variation and as a result a travel-time variation for the seismic waves that pass through the reservoir. A method based on the residual time shifts observed on reflectors below the reservoir from one survey to the other was developed and implemented. This method permitted a detection of a movement gas bubble boundary between April 1981 and January 1982 which agrees with drill hole data (water level measurements). 相似文献
3.
The Sebei gasfield is the largest biogas accumulation found in China and many reservoirs and seal rocks superposed on a syndepositional anticline in Quaternary. The biogas charging and dissipating process and its distribution have been a research focus for many years. The authors suggest a diffusing and accumulating model for the biogas, as they find that the shallower the gas producer, the more methane in the biogas, and the lighter stable carbon isotope composition of methane. Based on the diffusing model, diffused biogas is quantitatively estimated for each potential sandy reservoir in the gasfield, and the gas charging quantity for the sandy reservoir is also calculated by the diffused gas quantity plus gas reserve in-place. A ratio of diffusing quantity to charging quantity is postulated to describe biogas accumulating state in a sandy reservoir, if the ratio is less than 0.6, the reservoir forms a good gas-pool and high-production layer in the gasfield, which often occurs in the reservoirs deeper than 900 m; if the ratio is greater than 0.6, a few gas accumulated in the reservoir, which frequently exists in the reservoirs shallower than 900 m. Therefore, a biogas accumulation model is built up as lateral direct charging from gas source for the sands deeper than 900 m and indirect charging from lower gas-bearing sands by diffusion at depth shallower than 900 m. With this charging and diffusion quantitative model, the authors conducted re-evaluation on each wildcat in the central area of the Qaidam Basin, and found many commercial biogas layers. 相似文献
4.
《中国科学D辑(英文版)》2008,(Z1)
The Sebei gasfield is the largest biogas accumulation found in China and many reservoirs and seal rocks superposed on a syndepositional anticline in Quaternary.The biogas charging and dissipating process and its distribution have been a research focus for many years.The authors suggest a diffusing and accumulating model for the biogas,as they find that the shallower the gas producer,the more methane in the biogas,and the lighter stable carbon isotope composition of methane.Based on the diffusing model,diffused biogas is quantitatively estimated for each potential sandy reservoir in the gasfield,and the gas charging quantity for the sandy reservoir is also calculated by the diffused gas quantity plus gas reserve in-place.A ratio of diffusing quantity to charging quantity is postulated to describe biogas accumulating state in a sandy reservoir,if the ratio is less than 0.6,the reservoir forms a good gas-pool and high-production layer in the gasfield,which often occurs in the reservoirs deeper than 900 m;if the ratio is greater than 0.6,a few gas accumulated in the reservoir,which frequently exists in the reservoirs shallower than 900 m.Therefore,a biogas accumulation model is built up as lateral direct charging from gas source for the sands deeper than 900 m and indirect charging from lower gas-bearing sands by diffusion at depth shallower than 900 m.With this charging and diffusion quantitative model,the authors conducted re-evaluation on each wildcat in the central area of the Qaidam Basin,and found many commercial biogas layers. 相似文献
5.
The Sebei gasfield is the largest biogas accumulation found in China and many reservoirs and seal rocks superposed on a syndepositional anticline in Quaternary. The biogas charging and dissipating process and its distribution have been a research focus for many years. The authors suggest a diffusing and accumulating model for the biogas, as they find that the shallower the gas producer, the more methane in the biogas, and the lighter stable carbon isotope composition of methane. Based on the diffusing model, diffused biogas is quantitatively estimated for each potential sandy reservoir in the gasfield, and the gas charging quantity for the sandy reservoir is also calculated by the diffused gas quantity plus gas reserve in-place. A ratio of diffusing quantity to charging quantity is postulated to describe biogas accumulating state in a sandy reservoir, if the ratio is less than 0.6, the reservoir forms a good gas-pool and high-production layer in the gasfield, which often occurs in the reservoirs deeper than 900 m; if the ratio is greater than 0.6, a few gas accumulated in the reservoir, which frequently exists in the reservoirs shallower than 900 m. Therefore, a biogas accumulation model is built up as lateral direct charging from gas source for the sands deeper than 900 m and indirect charging from lower gas-bearing sands by diffusion at depth shallower than 900 m. With this charging and diffusion quantitative model, the authors conducted re-evaluation on each wildcat in the central area of the Qaidam Basin, and found many commercial biogas layers.
相似文献6.
The South Pars, the largest gas field in the world, is located in the Persian Gulf. Structurally, the field is part of the Qatar–South Pars arch which is a regional anticline considered as a basement-cored structure with long lasting passive folding induced by salt withdrawal. The gas-bearing reservoir belongs to Kangan and Dalan formations dominated by carbonate rocks. The fracture role is still unknown in gas accumulation and distribution in this reservoir. In this paper, the Scattering Index (SI) and the semblance methods based on scattered waves and diffraction signal studies, respectively, were used to delineate the fracture locations. To find the relation between fractures and gas distribution, desired facies containing the gas, were defined and predicted using a method based on Bayesian facies estimation. The analysis and combination of these results suggest that preference of fractures and/or fractured zones are negligible (about 1% of the total volume studied in this paper) and, therefore, it is hard to conceive that they play an important role in this reservoir. Moreover, fractures have no considerable role in gas distribution (less than 30%). It can be concluded from this study that sedimentary processes such as digenetic, primary porosities and secondary porosities are responsible for the gas accumulation and distribution in this reservoir. 相似文献
7.
P wave velocity structure of upper and middle crust beneath the Three Gorges reservoir dam and adjacent regions 总被引：3，自引：0，他引：3
The three-dimensional P wave velocity structures of upper and middle crust beneath the Three Gorges reservoir dam and adjacent
regions were reconstructed by means of the seismic tomography, using the absolute P-arrivals from the events observed by the
Three Gorges seismic network. Our preliminary study of the deep structure shows: (1) the lithology of the Huangling anticline,
integral and highly solidified, exhibits a uniformity. Down to a depth of 20 km, the Huangling anticline still shows high-velocity
anomalies. Note that the characteristics of the velocity images with the different depths are not the same; (2) the depth
of the sedimentary rocks beneath the entire Zigui basin is not constant, which is less than 6 km; (3) the top boundary of
the crystalline basement varies from 0 to 7 km, whereas the bottom one between 14 and 16 km; (4) the partially molten substances
of the upper mantle are located in the Badong-Zigui and its adjacent regions as well as nearby Peishi area. The development
and evolution of the current tectonics are associated with the cooling and solidifying process of the molten substances within
the upper mantle; (5) the different velocity structures of the deep crust in both east and west side of the Jiuwanxi fault
are to be investigated in order to understand the origin of gravitational gradient zone in the Three Gorge areas; and (6)
the velocities within the upper crust would be affected by reservoir water permeability.
Supported by National Natural Science Foundation of China (Grant No. 40574039), National Key Technology R& D Program (Grant
No. 2008BAC38B02) 相似文献
8.
Magnetotelluric surveys were conducted to investigate the structure in the El Hito Anticline in the southern Loranca Basin, Spain. The structure of this anticline is partly known from nearby data (geological cross sections, geological maps, well-logs) from vintage geophysical and geological basin-scale surveys. Unfortunately, these surveys do not have the appropriate resolution to determine certain characteristics of the anticline accurately, such as the thickness or geometry of geological units. To address this deficiency, magnetotelluric data were acquired at 51 sites along three profiles to image the electrical resistivity of the anticline. To identify a geologically reasonable resistivity model, a conceptual model based on previous geological and geophysical information (cross sections and well-log data) was generated. Several inversions were performed using the conceptual model, which played a key role in the interpretation of the magnetotelluric data, and the construction of a suitable initial model was essential in producing geologically meaningful models. Using these results, we obtained 2-D resistivity models that provide information on the main geological units and structures. We also performed sensitivity tests to understand the morphology and structure of the resistive basement better. The results indicate the presence of basement highs that show different structural styles for the basement and its cover, and an increase in the basement depth towards the south. The final 2-D resistivity models provide new information regarding the structure of the anticline and demonstrate a clear correlation between the main geological units and resistivity, as well as correlations with mapped surface faults. 相似文献
9.
Statistical analyses of the Karun III reservoir seismicity reveal a remarkable correlation between seismicity rate and water-level
harmonic changes. It seems that seismicity in this dam depends on rapid water-level changes. The three biggest earthquakes
of Karun III, measuring 4, 4.1, and 4.3 on the Richter scale (ML), occurred after two stages of rapid filling of the dam on
March 22, 2005 and May 12, 2006. These earthquakes happened when the water reached the maximum operational level. Since the
beginning of filling the reservoir on November 8, 2004 until March, 2006, most reservoir-induced seismicity has been localized
in three main clusters. The majority of the earthquakes occurred in the frontal anticline of Keyf Malek; the second and third
clusters happened near Karun Blind Fault (KBF) and Mountain Front Fault (MFF), respectively. Filling Karun III reservoir immediately
led to an increase in the occurrence of earthquakes. Further, following abrupt water-level changes, a considerable increase
in the number of earthquakes is observed. Finally, in terms of seismicity rate, vertical and horizontal migration, magnitude,
and distance, the earthquakes of Karun III behave differently. 相似文献
10.
熊坡背斜构造变形与蒲江-新津断裂活动特征 总被引：3，自引：0，他引：3
熊坡背斜位于龙门山构造带东南端的成都盆地内,是龙门山逆冲推覆构造向前推挤进入盆地内部的一个主要变形区域,与其配套发育的断裂为蒲江-新津断裂,断裂与背斜褶皱之间在构造变形模式上表现出明显的一致性。在褶皱和断裂的构造变形和活动特征上,熊坡背斜南段表现为一种不对称的褶皱,向NE方向发展表现为较为宽缓的对称褶皱形态,卷入的地层主要是中生代及其以前的地层,对蒲江-新津断裂的地貌调查结果表明,断裂没有对该区域内广泛发育的冲沟Ⅰ级阶地产生影响,而对山前发育的相当于南河(岷江Ⅰ级支流)Ⅳ级阶地的洪积台地有明显的控制作用,说明断裂活动时间应该为第四纪早期,到第四纪晚期活动减弱或是趋于静止 相似文献
11.
Estimating sediment deposition volume in a reservoir using the physiographic soil erosion-deposition model 总被引：1，自引：0，他引：1
A traditional way to estimate the amount of sediment deposition in a reservoir is simply based on the variation in the elevation measurements made at different times.In the present study,a physiographic soil erosion-deposition(PSED) model is applied to.Tseng-Wen reservoir watershed to compute the discharge and sediment concentration into the reservoir,and the deposition volume.The PSED model is verified using three storm events and then compared to the historic data recorded at the hydrological stations in the reservoir watershed,which shows close agreement between the simulated and recorded values of the flow discharge into the reservoir and the corresponding sediment yield. The numerical simulation was also conducted for all rainfall events from 1996 to 2001 to estimate the annual sediment deposition volume in the reservoir.Similarly,the simulated result was compared to that obtained from the cross-section elevation measurement.It is found that some measured data points had large uncertainty or were questionable.On the contrary,the model results and the measured data that renders the model credible.Therefore,regardless of the estimation for the flow discharge into the reservoir induced by a single rainfall event,the sediment yield,or the sediment deposition,the PSED model has successfully demonstrated its ability. 相似文献
12.
位于龙门山逆冲推覆构造带东侧的龙泉山背斜,构成了四川前陆盆地的前陆隆起。通过室内航空相片对凯江跨背斜段的地貌面的解译,结合野外考察可知凯江发育3级阶地,其中T1、T2为堆积阶地,T3为基座阶地。在野外用差分GPS测量了阶地的空间坐标信息,同时采集了各级阶地堆积物的测年样本,并经实验分析约束了阶地的形成年龄。另外,对石油地震剖面解译揭示出龙泉山背斜北段地壳缩短和隆升主要是通过褶皱膝折带迁移机制进行的,滑脱层的深度约6km。利用面积守恒准则计算出龙泉山背斜晚更新世以来的地壳缩短速率约为(1.36±0.41)mm/a、隆升速率为(0.64±0.19)mm/a。通过滑脱层的推覆抬升机制形成的龙泉山背斜,给青藏高原东缘变形模式中的逆断层推覆地壳缩短造山增加了证据。 相似文献
13.
位于帕米尔前缘逆冲推覆体(Pamir Front Thrust,PFT)东端的木什滑脱背斜,是帕米尔弧形推覆构造带最前缘和最新的变形带。对地形横剖面、纵剖面和水系发育特征的分析表明,木什背斜总体上具有由西向东扩展生长的特征。在背斜核部及北翼发育数级开阔平坦的沿轴向展布的河流阶地,阶地可划分为4期。利用阶地堆积细颗粒石英光释光测年获得阶地面T2a、T3和T4的形成年龄分别为(15.8±2.40)ka、(55.1±10.3)ka、(131.4±23.9)ka。伴随背斜的生长扩展,河流阶地面发生了横向和纵向掀斜,并形成断层陡坎和褶皱陡坎。木什背斜晚第四纪的缩短和隆升主要是通过褶皱翼旋转机制进行的,估算其最小缩短速率为(1.6±0.3)mm/a,最小隆升速率为(1.9±0.3)mm/a。与此同时,沿轴向背斜发生了向东的侧向迁移和旋转。根据背斜垂直隆升与侧向扩展之间的关系,估算背斜在131～16ka期间向东的侧向迁移扩展速率较快,为 (14.6±3.6)mm/a; 自16ka至今,侧向迁移扩展速率迅速减小至(1.7±0.3)mm/a,背斜向东的迁移扩展可能已基本停止,而以侧向旋转为主。 相似文献
14.
Tianshan is one of the longest and most active intracontinental orogenic belts in the world. Due to the collision between Indian and Eurasian plates since Cenozoic, the Tianshan has been suffering from intense compression, shortening and uplifting. With the continuous extension of deformation to the foreland direction, a series of active reverse fault fold belts have been formed. The Xihu anticline is the fourth row of active fold reverse fault zone on the leading edge of the north Tianshan foreland basin. For the north Tianshan Mountains, predecessors have carried out a lot of research on the activity of the second and third rows of the active fold-reverse faults, and achieved fruitful results. But there is no systematic study on the Quaternary activities of the Xihu anticline zone. How is the structural belt distributed in space?What are the geometric and kinematic characteristics?What are the fold types and growth mechanism?How does the deformation amount and characteristics of anticline change?In view of these problems, we chose Xihu anticline as the research object. Through the analysis of surface geology, topography and geomorphology and the interpretation of seismic reflection profile across the anticline, we studied the geometry, kinematic characteristics, fold type and growth mechanism of the structural belt, and calculated the shortening, uplift and interlayer strain of the anticline by area depth strain analysis.
In this paper, by interpreting the five seismic reflection profiles across the anticline belt, and combining the characteristics of surface geology and geomorphology, we studied the types, growth mechanism, geometry and kinematics characteristics, and deformation amount of the fold. The deformation length of Xihu anticline is more than 47km from west to east, in which the hidden length is more than 14km. The maximum deformation width of the exposed area is 8.5km. The Xihu anticline is characterized by small surface deformation, simple structural style and symmetrical occurrence. The interpretation of seismic reflection profile shows that the deep structural style of the anticline is relatively complex. In addition to the continuous development of a series of secondary faults in the interior of Xihu anticline, an anticline with small deformation amplitude(Xihubei anticline)is continuously developed in the north of Xihu anticline. The terrain high point of Xihu anticline is located about 12km west of Kuitun River. The deformation amplitude decreases rapidly to the east and decreases slowly to the west, which is consistent with the interpretation results of seismic reflection profile and the calculation results of shortening. The Xihu anticline is a detachment fold with the growth type of limb rotation. The deformation of Xihu anticline is calculated by area depth strain analysis method. The shortening of five seismic reflection sections A, B, C, D and E is(650±70) m, (1 070±70) m, (780±50) m, (200±40) m and(130±30) m, respectively. The shortening amount is the largest near the seismic reflection profile B of the anticline, and decreases gradually along the strike to the east and west ends of the anticline, with a more rapidly decrease to the east, which indicates that the topographic high point is also a structural high point. The excess area caused by the inflow of external material or outflow of internal matter is between -0.34km2 to 0.56km2. The average shortening of the Xihubei anticline is between(60±10) m and(130±40) m, and the excess area caused by the inflow of external material is between 0.50km2 and 0.74km2. The initial locations of the growth strata at the east part is about 1.9~2.0km underground, and the initial location of the growth strata at the west part is about 3.7km underground. We can see the strata overlying the Xihu anticline at 3.3km under ground, the strata above are basically not deformed, indicating that this section of the anticline is no longer active. 相似文献
In this paper, by interpreting the five seismic reflection profiles across the anticline belt, and combining the characteristics of surface geology and geomorphology, we studied the types, growth mechanism, geometry and kinematics characteristics, and deformation amount of the fold. The deformation length of Xihu anticline is more than 47km from west to east, in which the hidden length is more than 14km. The maximum deformation width of the exposed area is 8.5km. The Xihu anticline is characterized by small surface deformation, simple structural style and symmetrical occurrence. The interpretation of seismic reflection profile shows that the deep structural style of the anticline is relatively complex. In addition to the continuous development of a series of secondary faults in the interior of Xihu anticline, an anticline with small deformation amplitude(Xihubei anticline)is continuously developed in the north of Xihu anticline. The terrain high point of Xihu anticline is located about 12km west of Kuitun River. The deformation amplitude decreases rapidly to the east and decreases slowly to the west, which is consistent with the interpretation results of seismic reflection profile and the calculation results of shortening. The Xihu anticline is a detachment fold with the growth type of limb rotation. The deformation of Xihu anticline is calculated by area depth strain analysis method. The shortening of five seismic reflection sections A, B, C, D and E is(650±70) m, (1 070±70) m, (780±50) m, (200±40) m and(130±30) m, respectively. The shortening amount is the largest near the seismic reflection profile B of the anticline, and decreases gradually along the strike to the east and west ends of the anticline, with a more rapidly decrease to the east, which indicates that the topographic high point is also a structural high point. The excess area caused by the inflow of external material or outflow of internal matter is between -0.34km2 to 0.56km2. The average shortening of the Xihubei anticline is between(60±10) m and(130±40) m, and the excess area caused by the inflow of external material is between 0.50km2 and 0.74km2. The initial locations of the growth strata at the east part is about 1.9~2.0km underground, and the initial location of the growth strata at the west part is about 3.7km underground. We can see the strata overlying the Xihu anticline at 3.3km under ground, the strata above are basically not deformed, indicating that this section of the anticline is no longer active. 相似文献
15.
发育在帕米尔弧形推覆构造带最前缘的木什活动背斜是一南缓北陡的第四纪滑脱褶皱,背斜的最小地壳缩短量为0.7km,构造隆升幅度可达1.5km.木什背斜北翼逆断层由一系列坡向北的反向断层陡坎组成,不同断坎间垂直位移分布呈现此消彼长的特征,不论是整个北翼逆断层西段还是单条断坎,其垂直位移均呈东高西低的不对称分布,位移梯度东高西... 相似文献
16.
The Zafra de Záncara anticline (also known as the El Hito anticline), located in the Loranca Cenozoic Basin (part of the Tagus Basin, Central Spain), had been studied by several oil companies during the late 1960s and early 1970s. In 2009, within the ‘Plan for selection and characterization of suitable structures of CO2 geological storage’, this anticline was selected as a potential CO2 storage site. A preliminary three-dimensional geological model, based on five geological cross sections that were constrained with the interpretation of the available seismic profiles (that are rather old and do not have very good quality), was created. With the aim of improving the geological knowledge of the Zafra de Záncara anticline and helping to investigate the suitability of a nearby anticline, namely La Rambla, as another structural closure that might make it a possible CO2 storage site, a local gravity survey (1 station every km2) was carried out in the area, seven new geological cross sections, based on these existing seismic profiles and field geology, were build, and a new three-dimensional geological model that included both anticlines, improved through three-dimensional stochastic gravity inversion, was constructed. The densities needed for the geological formations of the model come from the analysis of rock samples, logging data from El Hito-1 drillhole and petrophysical information from Instituto Geológico y Minero de España database. The inversion has improved the knowledge about the geometry of the anticlines’ traps and seals as well as the geometry of the basement relief and the structural relationship between basement and cover. 相似文献
17.
18.
The ultimate goal of reservoir simulation in reservoir surveillance technology is to estimate long-term production forecasting and to plan development and management of petroleum fields. However, maintaining reliable reservoir models which honour available static and dynamic data, involve inherent risks due to the uncertainties in space and time of the distribution of hydrocarbons inside reservoirs. Recent applications have shown that these uncertainties can be reduced by quantitative integration of seismic data into the reservoir modelling workflows to identify which areas and reservoir attributes of the model should be updated. This work aims using seismic data to reduce ambiguity in calibrating reservoir flow simulation model with an uncertain petro-elastic model, proposing a circular workflow of inverted seismic impedance (3D and 4D) and engineering studies, with emphasis on the interface between static and dynamic models. The main contribution is to develop an updating procedure for adjusting reservoir simulation response before using it in the production forecasting and enhance the interpretive capability of reservoir properties. Accordingly, the workflow evaluates consistency of reservoir simulation model and inverted seismic impedance, assisted by production history data, to close the loop between reservoir engineering and seismic domains. The methodology is evaluated in a complex, faulted, sandstone reservoir, the Norne benchmark field, where a significant reservoir behaviour understanding (about the static and dynamic reservoir properties) is obtained towards the quantitative integration of seismic impedance data. This leads to diagnosis of the reservoir flow simulation reliability and generation of an updated simulation model consistent with observed seismic and well production history data, as well as a calibrated petro-elastic model. Furthermore, as Norne Field is a benchmark case, this study can be considered to enrich the discussions over deterministic or probabilistic history matching studies. 相似文献
19.
Probabilistic reservoir operation using Bayesian stochastic model and support vector machine 总被引：3，自引：0，他引：3
Available water resources are often not sufficient or too polluted to satisfy the needs of all water users. Therefore, allocating water to meet water demands with better quality is a major challenge in reservoir operation. In this paper, a methodology to develop operating strategies for water release from a reservoir with acceptable quality and quantity is presented. The proposed model includes a genetic algorithm (GA)-based optimization model linked with a reservoir water quality simulation model. The objective function of the optimization model is based on the Nash bargaining theory to maximize the reliability of supplying the downstream demands with acceptable quality, maintaining a high reservoir storage level, and preventing quality degradation of the reservoir. In order to reduce the run time of the GA-based optimization model, the main optimization model is divided into a stochastic and a deterministic optimization model for reservoir operation considering water quality issues.The operating policies resulted from the reservoir operation model with the water quantity objective are used to determine the released water ranges (permissible lower and upper bounds of release policies) during the planning horizon. Then, certain values of release and the optimal releases from each reservoir outlet are determined utilizing the optimization model with water quality objectives. The support vector machine (SVM) model is used to generate the operating rules for the selective withdrawal from the reservoir for real-time operation. The results show that the SVM model can be effectively used in determining water release from the reservoir. Finally, the copula function was used to estimate the joint probability of supplying the water demand with desirable quality as an evaluation index of the system reliability. The proposed method was applied to the Satarkhan reservoir in the north-western part of Iran. The results of the proposed models are compared with the alternative models. The results show that the proposed models could be used as effective tools in reservoir operation. 相似文献
20.
The behaviour of a magma plumbing system during a cycle of volcanic edifice growth is investigated with a simple physical model. Loading by an edifice at Earth's surface changes stresses in the upper crust and pressures in a magma reservoir. In turn, these changes affect magma ascent from a deep source to the reservoir and from reservoir to Earth's surface. The model plumbing system is such that a hydraulic connection is maintained at all times between the reservoir and a deep magma source at constant pressure. Consequently the input rate of magma into the reservoir is predicted by the model rather than imposed as an input parameter. The open hydraulic connection model is consistent with short-term measurements of deformation and seismicity at several active volcanoes. Threshold values for the reservoir pressure at the beginning and end of eruption evolve as the edifice grows and lead to long-term changes of eruption rate. Depending on the dimensions and depth of the reservoir, the eruption rate follows different trends as a function of time. For small reservoirs, the eruption rate initially increases as the edifice builds up and peaks at some value before going down. The edifice size at the peak eruption rate provides a constraint on the reservoir shape and depth. Edifice decay or destruction leads to resumption of eruptive activity and a new eruption cycle. A simple elastic model for country rock deformation is valid over a whole eruptive cycle extending to the cessation of eruptive activity. For large reservoirs, an elastic model is only valid over part of an eruptive cycle. Long-term stress changes eventually lead to reservoir instability in the form of either roof collapse and caldera formation or reservoir enlargement in the horizontal direction. 相似文献