首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Shear-wave velocity is a key parameter for calibrating monitoring time-lapse 4D seismic data during CO2-EOR (Enhanced Oil Recovery) and CO2 sequestration. However, actual S-wave velocity data are lacking, especially in 4D data for CO2 sequestration because wells are closed after the CO2 injection and seismic monitoring is continued but no well log data are acquired. When CO2 is injected into a reservoir, the pressure and saturation of the reservoirs change as well as the elastic parameters of the reservoir rocks. We propose a method to predict the S-wave velocity in reservoirs at different pressures and porosities based on the Hertz–Mindlin and Gassmann equations. Because the coordination number is unknown in the Hertz–Mindlin equation, we propose a new method to predict it. Thus, we use data at different CO2 injection stages in the Gao89 well block, Shengli Oilfield. First, the sand and mud beds are separated based on the structural characteristics of the thin sand beds and then the S-wave velocity as a function of reservoir pressure and porosity is calculated. Finally, synthetic seismic seismograms are generated based on the predicted P- and S-wave velocities at different stages of CO2 injection.  相似文献   

2.
We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO2 sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO2, methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.  相似文献   

3.
Causes and effects of global warming have been highly debated in recent years. Nonetheless, injection and storage of CO2 (CO2 sequestration) in the subsurface is becoming increasingly accepted as a viable tool to reduce the amount of CO2 from the atmosphere, which is a primary contributor to global warming. Monitoring of CO2 movement with time is essential to ascertain that sequestration is not hazardous. A method is proposed here to appraise CO2 saturation from seismic attributes using differential effective medium theory modified for pressure (PDEM). The PDEM theory accounts pressure-induced fluid flow between cavities, which is a very important investigation in the CO2-sequestered regime of heterogeneous microstructure. The study area is the lower Tuscaloosa formation at Cranfield in Mississippi, USA, which is one of the active enhanced oil recovery (EOR), and CO2 capture and storage (CCS) fields. Injection well (F1) and two observation wells (F2 and F3) are present close (within 112 m) to the detailed area of study for this region. Since the three wells are closely situated, two wells, namely injection well F1 and the furthest observation well F3, have been focused on to monitor CO2 movement. Time-lapse (pre- and post-injection) log, core and surface seismic data are used in the quantitative assessment of CO2 saturation from the PDEM theory. It has been found that after approximately 9 months of injection, average CO2 saturations in F1 and F3 are estimated as 50% in a zone of thickness ~ 25 m at a depth of ~ 3 km.  相似文献   

4.
The phase velocity of seismic waves varies with the propagation frequency, and thus frequency-dependent phenomena appear when CO2 gas is injected into a reservoir. By dynamically considering these phenomena with reservoir conditions it is thus feasible to extract the frequency-dependent velocity factor with the aim of monitoring changes in the reservoir both before and after CO2 injection. In the paper, we derive a quantitative expression for the frequency-dependent factor based on the Robinson seismic convolution model. In addition, an inversion equation with a frequency-dependent velocity factor is constructed, and a procedure is implemented using the following four processing steps: decomposition of the spectrum by generalized S transform, wavelet extraction of cross-well seismic traces, spectrum equalization processing, and an extraction method for frequency-dependent velocity factor based on the damped least-square algorithm. An attenuation layered model is then established based on changes in the Q value of the viscoelastic medium, and spectra of migration profiles from forward modeling are obtained and analyzed. Frequency-dependent factors are extracted and compared, and the effectiveness of the method is then verified using a synthetic data. The frequency-dependent velocity factor is finally applied to target processing and oil displacement monitoring based on real seismic data obtained before and after CO2 injection in the G89 well block within Shengli oilfield. Profiles and slices of the frequency-dependent factor determine its ability to indicate differences in CO2 flooding, and the predicting results are highly consistent with those of practical investigations within the well block.  相似文献   

5.
At present, using Eddy Covariance (EC) method to estimate the “true value” of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u*, we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u*c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.  相似文献   

6.
Tracking the migration of the CO2 plume is essential in order to better manage the operation of geologic sequestration of CO2. However, the large cost of most monitoring technologies, such as time-lapse seismic, limits their application. We investigated the application of a probabilistic history matching methodology using routine measurements of injection data, which are affected by the presence of large-scale heterogeneities, as an inexpensive alternative to track the migration of CO2 plume in an aquifer. The approach is demonstrated first through a synthetic example in which the ability to detect the presence of flow barriers was tested. In a second example, we applied our method to the In Salah field, one of the largest geological sequestration projects in the world, where the main direction of high permeability features was inferred. The accuracy and reproducibility of the results show that our approach for assisted history matching is an economic and viable option for plume monitoring during geologic CO2 sequestration.  相似文献   

7.
Use of eddy covariance (EC) techniques to map the spatial distribution of diffuse volcanic CO2 fluxes and quantify CO2 emission rate was tested at the Horseshoe Lake tree-kill area on Mammoth Mountain, California, USA. EC measurements of CO2 flux were made during September–October 2010 and ranged from 85 to 1,766 g m−2 day−1. Comparative maps of soil CO2 flux were simulated and CO2 emission rates estimated from three accumulation chamber (AC) CO2 flux surveys. Least-squares inversion of measured eddy covariance CO2 fluxes and corresponding modeled source weight functions recovered 58–77% of the CO2 emission rates estimated based on simulated AC soil CO2 fluxes. Spatial distributions of modeled surface CO2 fluxes based on EC and AC observations showed moderate to good correspondence (R 2 = 0.36 to 0.70). Results provide a framework for automated monitoring of volcanic CO2 emissions over relatively large areas.  相似文献   

8.
Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO2 sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO2 sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO2 sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO2 sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.  相似文献   

9.
The global rate of fossil fuel combustion continues to rise, but the amount of CO2 accumulating in the atmosphere has not increased accordingly. The causes for this discrepancy are widely debated. Particularly, the location and drivers for the interannual variability of atmospheric CO2 are highly uncertain. Here we examine links between global atmospheric CO2 growth rate (CGR) and the climate anomalies of biomes based on (1986–1995) global climate data of ten years and accompanying satellite data sets. Our results show that four biomes, the tropical rainforest, tropical savanna, C4 grassland and boreal forest, and their responses to climate anomalies, are the major climate-sensitive CO2 sinks/sources that control the CGR. The nature and magnitude by which these biomes respond to climate anomalies are generally not the same. However, one common influence did emerge from our analysis; the extremely high CGR observed for the one extreme El Niño year was caused by the response of the tropical biomes (rainforest, savanna and C4 grassland) to temperature.  相似文献   

10.
Soil temperature and gas (CO2 concentration and flux) have been investigated at Merapi volcano (Indonesia) during two inter-eruptive periods (2002 and 2007). Precise imaging of the summit crater and the spatial pattern of diffuse degassing along a gas traverse on the southern slope are interpreted in terms of summit structure and major caldera organization. The summit area is characterized by decreasing CO2 concentrations with distance from the 1932 crater rim, down to atmospheric levels at the base of the terminal cone. Similar patterns are measured on any transect down the slopes of the cone. The spatial distribution of soil gas anomalies suggests that soil degassing is controlled by structures identified as concentric historical caldera rims (1932, 1872, and 1768), which have undergone severe hydrothermal self-sealing processes that dramatically lower the permeability and porosity of soils. Temperature and CO2 flux measurements in soils near the dome display heterogeneous distributions which are consistent with a fracture network identified by previous geophysical studies. These data support the idea that the summit is made of isolated and mobile blocks, whose boundaries are either sealed by depositional processes or used as pathways for significant soil degassing. Within this context, self-sealing both prevents long-distance soil degassing and controls heat and volatile transfers near the dome. A rough estimate of the CO2 output through soils near the dome is 200–230 t day−1, i.e. 50% of the estimated total gas output from the volcano summit during these quiescent periods. On Merapi’s southern slope, a 2,500 m long CO2 traverse shows low-amplitude anomalies that fit well with a recently observed electromagnetic anomaly, consistent with a faulted structure related to an ancient avalanche caldera rim. Sub-surface soil permeability is the key parameter that controls the transfer of heat and volatiles within the volcano, allowing its major tectonic architecture to be revealed by soil gas and soil temperature surveys.  相似文献   

11.
On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO2 and SO2) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO2 flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO2 fluxes of the degassing plume and the CO2 fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO2 flux) and in soil degassing (CO2 flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.  相似文献   

12.
Soil H2 and CO2 surveys were carried out along seven active faults and around the aftershock region of the 2000 Tottori-ken Seibu earthquake in Japan. Diffuse CO2 effluxes were also measured along one fault and around the 2000 aftershock region. The results show highly variable H2 concentration in space and time and it seems that the maximum H2 concentration at each active fault correlates with fault activity as exemplified by the time of the latest big earthquakes. Even though observed H2 concentrations in four faults were markedly lower than those collected previously in the latter half of the 1970s, it is evident that the higher H2 concentrations in this study are due to the addition of the fault gases. Comparing the chemical composition of trapped gases (H2: 5–20% and CO2/H2: 0.5–12) in fractured rocks of drill cores bored at the Nojima fault, a soil gas sample with the highest H2 concentration showed large amounts of the trapped fault gas, diluted with atmospheric component. The profile experiment across a fracture zone at the Yamasaki fault showed higher H2 concentrations and lower CO2/H2 ratios as was observed in soil gas from the fracture zone. A few days after the 2000 Tottori-kei Seibu earthquake, no CO2 effluxes related to the occurrence of earthquakes were observed at the aftershock region. However, only above the epicenter zone, relatively high H2 concentrations in soil gases were observed.  相似文献   

13.
Carbon capture and storage is a viable greenhouse gas mitigation technology and the Sleipner CO2 sequestration site in the North Sea is an excellent example. Storage of CO2 at the Sleipner site requires monitoring over large areas, which can successfully be accomplished with time lapse seismic imaging. One of the main goals of CO2 storage monitoring is to be able to estimate the volume of the stored CO2 in the reservoir. This requires a parametrization of the subsurface as exact as possible. Here we use elastic 2D time‐domain full waveform inversion in a time lapse manner to obtain a P‐wave velocity constrain directly in the depth domain for a base line survey in 1994 and two post‐injection surveys in 1999 and 2006. By relating velocity change to free CO2 saturation, using a rock physics model, we find that at the considered location the aquifer may have been fully saturated in some places in 1999 and 2006.  相似文献   

14.
Seasonal metrics and environmental responses to forestry soil surface CO2 emission effluxes among three types of lower subtropical forests were consistently monitored over two years with static chamber-gas chromatograph techniques among three types of lower subtropical forests. Results showed that annual CO2 effluxes (S+L) reached 3942.20, 3422.36 and 2163.02 CO2 g·m?2·a?1, respectively in the monsoon evergreen broadleaf forest, mixed broadleaf-coniferous forest and coniferous forest. All the three types of forests revealed the same characteristics of seasonal changes with the CO2 effluxes peaking throughout June to August. During this peaking period, the effluxes were 35.9%, 38.1% and 40.2% of the total annual effluxes, respectively. The CO2 emission process responding to the environmental factors displayed significantly different patterns in forestry soils of the three types of forests. The coniferous forest (CF) was more sensitive to temperature than the other two types. The Q 10 values were higher, along with greater seasonal variations of the CO2 efflux, indicating that the structurally unique forestry ecosystem has disadvantage against interferences. All the three types of forestry CO2 effluxes showed significant correlation with the soil temperature (T s), soil water content (M s) and air pressure (P a). However, stepwise regression analysis indicated no significant correlation between air pressure and the soil CO2 efflux. With an empirical model to measure soil temperature and water content in 5 cm beneath the soil surface, the CO2 effluxes accounting for 75.7%, 77.8% and 86.5% of the efflux variability respectively in soils of BF, MF and PF were calculated. This model can be better used to evaluate the CO2 emission of soils under water stress and arid or semi-arid conditions.  相似文献   

15.
Long-term measurement of carbon metabolism of old-growth forests is critical to predict their behaviors and to reduce the uncertainties of carbon accounting under changing climate. Eddy covariance technology was applied to investigate the long-term carbon exchange over a 200 year-old Chinese broad-leaved Korean pine mixed forest in the Changbai Mountains (128°28′E and 42°24′N, Jilin Province, P. R. China) since August 2002. On the data obtained with open-path eddy covariance system and CO2 profile measurement system from Jan. 2003 to Dec. 2004, this paper reports (i) annual and seasonal variation of F NEE, F GPP and R E; (ii) regulation of environmental factors on phase and amplitude of ecosystem CO2 uptake and release Corrections due to storage and friction velocity were applied to the eddy carbon flux.LAI and soil temperature determined the seasonal and annual dynamics of FGPP and RE separately. VPD and air temperature regulated ecosystem photosynthesis at finer scales in growing seasons. Water condition at the root zone exerted a significant influence on ecosystem maintenance carbon metabolism of this forest in winter.The forest was a net sink of atmospheric CO2 and sequestered ?449 g C·m?2 during the study period; ?278 and ?171 gC·m?2 for 2003 and 2004 respectively. F GPP and F RE over 2003 and 2004 were ?1332, ?1294 g C·m?2. and 1054, 1124 g C·m?2 respectively. This study shows that old-growth forest can be a strong net carbon sink of atmospheric CO2.There was significant seasonal and annual variation in carbon metabolism. In winter, there was weak photosynthesis while the ecosystem emitted CO2. Carbon exchanges were active in spring and fall but contributed little to carbon sequestration on an annual scale. The summer is the most significant season as far as ecosystem carbon balance is concerned. The 90 days of summer contributed 66.9, 68.9% of F GPP, and 60.4, 62.1% of R E of the entire year.  相似文献   

16.
Lake Albano (Alban Hills volcanic complex, Central Italy) is located in a densely populated area near Rome. The deep lake waters have significant dissolved CO2 concentrations, probably related to sub-lacustrine fluid discharges fed by a pressurized CO2-rich reservoir. The analytical results of geochemical surveys carried out in 1989–2010 highlight the episodes of CO2 removal from the lake. The total mass of dissolved CO2 decreased from ∼5.8 × 107 kg in 1989 to ∼0.5 × 107 kg in 2010, following an exponential decreasing trend. Calculated values of both dissolved inorganic carbon and CO2 concentrations along the vertical profile of the lake indicate that this decrease is caused by CO2 release from the epilimnion, at depth <9 m, combined with (1) water circulation at depth <95 m and (2) CO2 diffusion from the deeper lake layers. According to this model, Lake Albano was affected by a large CO2 input that coincided with the last important seismic swarm at Alban Hills in 1989, suggesting an intimate relationship between the addition of deep-originated CO2 to the lake and seismic activity. In the case of a CO2 degassing event of an order of magnitude larger than the one that occurred in 1989, the deepest part of Lake Albano would become CO2-saturated, resulting in conditions compatible with the occurrence of a gas outburst. These results reinforce the idea that a sudden CO2 input into the lake may cause the release of a dense gas cloud, presently representing the major volcanic threat for this densely populated area.  相似文献   

17.
Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy   总被引：1，自引：0，他引：1  
At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO2–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO2 discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d–1. The steam associated with the CO2 output, computed assuming that the original H2O/CO2 ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×1012 J d–1 (17 MW). The location of the CO2 and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO2 degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO2 flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara  相似文献   

18.
Between the years 2008 and 2013, approximately 67 kilotons of CO2 have been injected at the Ketzin site, Germany. As part of the geophysical monitoring programme, time‐lapse electrical resistivity tomography has been applied using crosshole and surface‐downhole measurements of electrical resistivity tomography. The data collection of electrical resistivity tomography is partly based on electrodes that are permanently installed in three wells at the site (one injection well and two observation wells). Both types of ERT measurements consistently show the build‐up of a CO2‐related resistivity signature near the injection point. Based on the imaged resistivity changes and a petrophysical model, CO2 saturation levels are estimated. These CO2 saturations are interpreted in conjunction with CO2 saturations inferred from neutron‐gamma loggings. Apart from the CO2–brine substitution response in the observed resistivity changes, significant imprints from the dynamic behaviour of the CO2 in the reservoir are observed.  相似文献   

19.
We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO2 sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO2 injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO2 dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO2, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO2 plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

20.
We describe analytical details and uncertainty evaluation of a simple technique for the measurement of the carbon isotopic composition of CO2 in volcanic plumes. Data collected at Solfatara and Vulcano, where plumes are fed by fumaroles which are accessible for direct sampling, were first used to validate the technique. For both volcanoes, the plume-derived carbon isotopic compositions are in good agreement with the fumarolic compositions, thus providing confidence on the method, and allowing its application at volcanoes where the volcanic component is inaccessible to direct sampling. As a notable example, we applied the same method to Mount Etna where we derived a δ13C of volcanic CO2 between −0.9 ± 0.27‰ and −1.41 ± 0.27‰ (Bocca Nuova and Voragine craters). The comparison of our measurements to data reported in previous work highlights a temporal trend of systematic increase of δ13C values of Etna CO2 from ~ −4‰, in the 1970’s and the 1980’s, to ~ −1‰ at the present time (2009). This shift toward more positive δ13C values matches a concurrent change in magma composition and an increase in the eruption frequency and energy. We discuss such variations in terms of two possible processes: magma carbonate assimilation and carbon isotopic fractionation due to magma degassing along the Etna plumbing system. Finally, our results highlight potential of systematic measurements of the carbon isotopic composition of the CO2 emitted by volcanic plumes for a better understanding of volcanic processes and for improved surveillance of volcanic activity.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号