共查询到20条相似文献，搜索用时 593 毫秒
1.
2.
The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image. 相似文献
3.
利用逆重复M序列伪随机信号良好的自相关特性,提出相关检测地电传输特性的编码电磁勘探法.分析了编码电磁测深原理及参数提取方法.同时记录发送电流信号和多收发距电磁场响应,对源信号和场信号进行相关运算,解卷积分离接收系统响应后可得到大地的频率特性或时间特性,以此实现地电断面的精细探测.在时间域,大地冲激响应和阶跃响应含有丰富的地电信息.通过冲激响应的峰值时间或阶跃响应的晚期渐近值可估计地电阻率分布.基于层状模型的大地冲激响应和阶跃响应正演计算结果表明,编码电磁测深法对大埋深薄层目标体有精细的分辨能力.可应用于能源、矿产资源、水资源、环境地质及工程地质勘察,有良好的应用前景. 相似文献
4.
从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线. 相似文献
5.
Fast-decaying IP in frozen unconsolidated rocks and potentialities for its use in permafrost-related TEM studies 总被引：1，自引：0，他引：1
We investigate the early time induced polarization (IP) phenomenon in frozen unconsolidated rocks and its association with transient electromagnetic (TEM) signals measured in northern regions. The distinguishing feature of these signals is the distortion of the monotony or sign reversals in the time range from a few tens to a few hundreds of microseconds. In simulating TEM data, the IP effects in frozen ground were attributed to the dielectric relaxation phenomenon rather than to the frequency‐dependent conductivity. This enabled us to use laboratory experimental data available in the literature on dielectric spectroscopy of frozen rocks. In our studies we focused on simulating the transient response of a coincident‐loop configuration in three simple models: (i) a homogeneous frozen earth (half‐space); (ii) a two‐layered earth with the upper layer frozen; (iii) a two‐layered earth with the upper layer unfrozen. The conductivities of both frozen and unfrozen ground were assumed to exhibit no frequency dispersion, whereas the dielectric permittivity of frozen ground was assumed to be described by the Debye model. To simplify the presentation and the comparison analysis of the synthetic data, the TEM response of a frozen polarizable earth was normalized to that of a non‐polarizable earth having the same structure and resistivities as the polarizable earth. The effect of the dielectric relaxation on a TEM signal is marked by a clearly defined minimum. Its time coordinate tmin is approximately three times larger than the dielectric relaxation time constant τ. This suggests the use of tmin for direct estimation of τ, which, in turn, is closely associated with the temperature of frozen unconsolidated rock. The ordinate of the minimum is directly proportional to the static dielectric permittivity of frozen earth. Increasing the resistivity of a frozen earth and/or decreasing the loop size results in a progressively stronger effect of the dielectric relaxation on the TEM signal. In the case of unfrozen earth, seasonal freezing is not likely to have an appreciable effect on the TEM signal. However, for the frozen earth, seasonal thawing of a near‐surface layer may result in a noticeable attenuation of the TEM signal features associated with dielectric relaxation in a frozen half‐space. Forward calculations show that the dielectric relaxation of frozen unconsolidated rocks may significantly affect the transient response of a horizontal loop laid on the ground. This conclusion is in agreement with a practical example of inversion of the TEM data measured over the permafrost. 相似文献
6.
Abdelghani Boudhar Gilles Boulet Lahoucine Hanich Jean Emmanuel Sicart Abdelghani Chehbouni 《水文科学杂志》2013,58(5):931-943
ABSTRACTIn this study, we used an energy balance model and two simple methods based on readily available data to identify the processes driving the point-scale energy and mass balance of the snowpack. Data were provided from an experimental site located at 3200 m. All models were evaluated by comparing observed and modelled snow water equivalents. Performances are variable from one season to the next and the energy balance model gives better results (mean of root mean square error, RMSE = 25 mm and r2 = 0.90) than the two simplified approaches (mean of RMSE = 54 mm and r2 = 0.70). There are significant amounts of snow sublimation but they are highly variable from season to season, depending on wind conditions (between 7 and 20% of the total). While the main source of energy for melting is net radiation, the amount of heat brought by sensible heat flux is significant for two of the most windy snow seasons.
Editor Z.W. Kundzewicz Associate editor not assigned 相似文献
7.
将伪随机编码技术引入到人工源电磁法后,可以通过加大发射功率以及应用后续的相关处理技术来达到压制噪声、加大探测深度及提高分辨率的目的,因此引起了越来越多学者的关注及研究,但大多数研究集中在资料处理的相关技术上,对模拟电磁信号关注较少.然而,资料处理工作大多是从电磁信号出发的,模拟伪随机编码源激发下的电磁信号不但可以为资料处理环节提供理论数据,而且可以为检测资料处理的效果提供中间结果,因此,模拟电磁信号工作必不可少.本文根据获得接收信号的物理过程来实现伪随机编码源激发下的电磁信号合成.首先用解析公式获得特定地电结构的大地频率域响应,然后通过余弦变换得到时间域阶跃响应,接下来用阶跃响应的时间导数得到大地脉冲响应,通过将大地脉冲响应与伪随机编码源的褶积得到理想接收信号,最后,用低通滤波器来模拟发射设备和接收设备的频带限制,将之和噪声一起加到理想接收信号上,最终模拟出仿真的合成信号.通过和野外实际接收信号对比发现本文合成信号仿真度较高,可以服务于后续的数据处理环节. 相似文献
8.
Anwar H. Bhuyian Brd P. Thrane Martin Landr Stle E. Johansen 《Geophysical Prospecting》2010,58(3):505-533
A comprehensive controlled source electromagnetic (CSEM) modelling study, based on complex resistivity structures in a deep marine geological setting, is conducted. The study demonstrates the effects of acquisition parameters and multi‐layered resistors on CSEM responses. Three‐dimensional (3D) finite difference time domain (FDTD) grid‐modelling is used for CSEM sensitivity analysis. Interpolation of normalized CSEM responses provides attributes representing relative sensitivity of the modelled structures. Modelling results show that fine grid, 1 × 1 km receiver spacing, provides good correlations between CSEM responses and the modelled structures, irrespective of source orientation. The resolution of CSEM attributes decreases for receiver spacing >2 × 2 km, when using only in‐line data. Broadside data in the grid geometry increase data density by 100 – approximately 200% by filling in in‐line responses and improve the resolution of CSEM attributes. Optimized source orientation (i.e., oblique to the strike of an elongated resistor) improves the structural definition of the CSEM anomalies for coarse‐grid geometries (receiver spacing ≥3 × 3 km). The study also shows that a multi‐resistor anomaly is not simply the summation but a cumulative response with mutual interference between constituent resistors. The combined response of constituent resistors is approximately 50% higher than the cumulative response of the multi‐resistor for 0.5 Hz at 4000 m offset. A gradual inverse variation of offset and frequency allows differentiation of CSEM anomalies for multi‐layered resistors. Similar frequency‐offset variations for laterally persistent high‐resistivity facies show visual continuity with varying geometric expressions. 3D grid‐modelling is an effective and adequate tool for CSEM survey design and sensitivity analysis. 相似文献
9.
Huijuan Cui Vijay P. Singh 《Stochastic Environmental Research and Risk Assessment (SERRA)》2016,30(6):1545-1563
This paper develops a minimum relative entropy theory with frequency as a random variable, called MREF henceforth, for streamflow forecasting. The MREF theory consists of three main components: (1) determination of spectral density (2) determination of parameters by cepstrum analysis, and (3) extension of autocorrelation function. MREF is robust at determining the main periodicity, and provides higher resolution spectral density. The theory is evaluated using monthly streamflow observed at 20 stations in the Mississippi River basin, where forecasted monthly streamflows show the coefficient of determination (r 2) of 0.876, which is slightly higher in the Upper Mississippi (r 2 = 0.932) than in the Lower Mississippi (r 2 = 0.806). Comparison of different priors shows that the prior with the background spectral density with a peak at 1/12 frequency provides satisfactory accuracy, and can be used to forecast monthly streamflow with limited information. Four different entropy theories are compared, and it is found that the minimum relative entropy theory has an advantage over maximum entropy (ME) for both spectral estimation and streamflow forecasting, if additional information as a prior is given. Besides, MREF is found to be more convenient to estimate parameters with cepstrum analysis than minimum relative entropy with spectral power as random variable (MRES), and less information is needed to assume the prior. In general, the reliability of monthly streamflow forecasting from the highest to the lowest is for MREF, MRES, configuration entropy (CE), Burg entropy (BE), and then autoregressive method (AR), respectively. 相似文献
10.
Impact of tributary DOM and nutrient inputs on the nearshore ecology of a large, oligotrophic lake (Georgian Bay, Lake Huron, Canada) 总被引：1，自引：0，他引：1
Serghei A. Bocaniov David R. Barton Sherry L. Schiff Ralph E. H. Smith 《Aquatic Sciences - Research Across Boundaries》2013,75(2):321-332
Two adjacent bays in a large oligotrophic lake (Georgian Bay, Lake Huron) were compared to determine how the inputs from relatively pristine, but moderately humic, tributaries may influence phytoplankton, nutrients and system metabolism. Dissolved organic carbon (DOC) concentrations decreased from 4 to 5 gC m?3 at inner sites to 2 gC m?3 or less at outer sites. The concentration gradient from inner to outer was greater in the bay with a major tributary, and optical properties (intensity and slope of light absorption spectrum) showed there was a loss of material with allochthonous characteristics along the gradient. Chlorophyll a (Chl a) and total phosphorus (TP) were also higher at inner (2–4 mg Chl a m?3 and 8–12 mgP m?3, respectively) than outer sites (≤1 mg Chl a m?3 and 4–5 mgP m?3). Chl a and TP, as well as particulate nutrient ratios (C:P, C:N, N:P), indicated significant eutrophication at inner sites, especially in the bay with the tributary, and there was a strong positive Chl a-phosphorus relationship. The stable oxygen isotope ratio (18O:16O) of dissolved oxygen indicated greater influence of biological oxygen fluxes at inner sites (where ratios were 2–3 ppt below atmospheric equilibrium) than at outer sites (where ratios were within 0.5 ppt of equilibrium). Community photosynthesis:respiration ratios inferred from 18O:16O varied positively with Chl a and inorganic nutrients, but negatively with DOC. Altered loading of allochthonous organic matter can be expected under changing climate and development scenarios and will have significant influence on optical properties and system metabolism through changes in DOC in this coastal system. The effects will nonetheless be strongly modulated by any accompanying change in inorganic nutrients. 相似文献
11.
L. J. Tsay 《Pure and Applied Geophysics》1973,110(1):1982-1995
Summary In a four-element model, which consists of two Voigt elements in series, the impulse responses of spherical waves for particle displacement, velocity, acceleration and strain are computed numerically. These impulse responses are shown to be useful in synthesis for the desired wave response through convolution, between the input of the source function and the impulse responses. The rate of decay of peak amplitude versus distance, as waves radiate out, is examined in the region close to the source which may be nuclear or dynamite explosives and found to be linear in log-log scale. The amplitude spectra of the unit impulse response are obtained through Fourier analysis, and the model is found to behave as a low pass filter. With the properly chosen elastic constants for springs in the model and damping coefficients for dashpot, it is seen that, like the earth materials, the attenuation exponent (in neper/1000 ft) versus frequency in the model is approximately linear. One of the results observed byCollins andLee [1] is fitted by the model. The model thus simulates more nearly the wave in earth materials close to the source than the Voigt solid. 相似文献
12.
The attenuation of coda waves in the earth’s crust in southwest (SW) Anatolia is estimated by using the coda wave method, which is based on the decrease of coda wave amplitude in time and distance. A total of 159 earthquakes were recorded between 1997 and 2010 by 11 stations belonging to the KOERI array. The coda quality factor Q c is determined from the properties of scattered coda waves in a heterogeneous medium. Firstly, the quality factor Q 0 (the value of Q c at 1 Hz.) and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves for frequencies of 1.5, 3.0, 6.0, 8.0, 12 and 20 Hz. Secondly, the attenuation coefficients (δ) are estimated. The shape of the curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The average Q c values vary from 110 ± 15 to 1,436 ± 202 for the frequencies above. The Q 0 and η values vary from 63 ± 7 to 95 ± 10 and from 0.87 ± 0.03 to 1.04 ± 0.09, respectively, for SW Anatolia. In this region, the average coda Q–f relation is described by Q c = (78 ± 9)f 0.98±0.07 and δ = 0.012 km?1. The low Q 0 and high η are consistent with a region characterized by high tectonic activity. The Q c values were correlated with the tectonic pattern in SW Anatolia. 相似文献
13.
海洋可控源电磁法(mCSEM)的时间域冲激响应特征可以反映海底油气高阻薄层.本文计算了水平电偶极子源均匀大地半空间,海洋均匀双半空间和海洋四层模型的阶跃响应和冲激响应,提出了瞬变冲激时刻的概念.分析了水平电偶源瞬变冲激时刻与介质电导率的指示关系.对于海底油气高阻薄层宜采用多偏移距同时测量方式,由于在低电导率介质中电磁能量传播得要快,在适当的收发距瞬变冲激时刻会提前到达,提出的瞬变冲激时刻道间变化量可以明确指示高阻薄层的存在及埋深.文中还分析了海水深度对瞬变冲激时刻的影响.由于“天波”干扰,瞬变冲激响应受到一定收发距观测的限制.消除 “天波”影响是时间域和频率域mCSEM数据处理的研究热点. 相似文献
14.
Shivakumar Khaple Prashanth Reddy Hanmaiahgari Roberto Gaudio Subhasish Dey 《Acta Geophysica》2017,65(1):29-46
In this study, three kinds of pier arrangements were tested. They are (i) two piers in tandem, (ii) two piers in staggered arrangement, and (iii) three piers in symmetrically staggered arrangements. In the arrangement of two piers in tandem, the equilibrium scour depth at downstream pier decreases with an increase in downstream distance up to approximately eight times pier diameter and then increases with further increase in downstream distance. However, the scour depth at downstream pier is always smaller than that at upstream pier. In the arrangement of two staggered piers, the scour depth at the downstream pier for L/b = 4, where L is the offset distance and b is the pier diameter, is the same as that of the upstream pier at S = 8b, where S is the streamwise spacing or distance between piers. Further, for three piers in staggered arrangement, as the lateral spacing between downstream piers increases, the equilibrium scour depth at downstream pier decreases. 相似文献
15.
H. PASSALACQUA 《Geophysical Prospecting》1983,31(6):945-976
Using approximate boundary conditions, expressions for electromagnetic fields have been derived for a thin, highly resistive layer lying between two homogeneous layers excited by an electric dipole grounded on the surface of the earth. The variations of the fields with the parameter T/T1 (ratio of the transverse resistance of the thin layer to the transverse resistance of the first layer) were studied in relation to frequency, time, the normalized separation source—receiver, and the angle between the source and the radius to the observation point. For a value of h2/h1 (ratio of thickness of second layer to the thickness of the first layer) approximately equal to 0.2, the general three-layer medium case gives the same results as this approach. It was found that the electric fields have a very strong dependence on the parameter T (transverse resistance) which characterizes the thin, highly resistive layer. However, the magnetic fields depend only very weakly on this parameter. 相似文献
16.
Time series in the Earth Sciences are often characterized as self-affine long-range persistent, where the power spectral density, S, exhibits a power-law dependence on frequency, f, S(f) ~ f ?β , with β the persistence strength. For modelling purposes, it is important to determine the strength of self-affine long-range persistence β as precisely as possible and to quantify the uncertainty of this estimate. After an extensive review and discussion of asymptotic and the more specific case of self-affine long-range persistence, we compare four common analysis techniques for quantifying self-affine long-range persistence: (a) rescaled range (R/S) analysis, (b) semivariogram analysis, (c) detrended fluctuation analysis, and (d) power spectral analysis. To evaluate these methods, we construct ensembles of synthetic self-affine noises and motions with different (1) time series lengths N = 64, 128, 256, …, 131,072, (2) modelled persistence strengths β model = ?1.0, ?0.8, ?0.6, …, 4.0, and (3) one-point probability distributions (Gaussian, log-normal: coefficient of variation c v = 0.0 to 2.0, Levy: tail parameter a = 1.0 to 2.0) and evaluate the four techniques by statistically comparing their performance. Over 17,000 sets of parameters are produced, each characterizing a given process; for each process type, 100 realizations are created. The four techniques give the following results in terms of systematic error (bias = average performance test results for β over 100 realizations minus modelled β) and random error (standard deviation of measured β over 100 realizations): (1) Hurst rescaled range (R/S) analysis is not recommended to use due to large systematic errors. (2) Semivariogram analysis shows no systematic errors but large random errors for self-affine noises with 1.2 ≤ β ≤ 2.8. (3) Detrended fluctuation analysis is well suited for time series with thin-tailed probability distributions and for persistence strengths of β ≥ 0.0. (4) Spectral techniques perform the best of all four techniques: for self-affine noises with positive persistence (β ≥ 0.0) and symmetric one-point distributions, they have no systematic errors and, compared to the other three techniques, small random errors; for anti-persistent self-affine noises (β < 0.0) and asymmetric one-point probability distributions, spectral techniques have small systematic and random errors. For quantifying the strength of long-range persistence of a time series, benchmark-based improvements to the estimator predicated on the performance for self-affine noises with the same time series length and one-point probability distribution are proposed. This scheme adjusts for the systematic errors of the considered technique and results in realistic 95 % confidence intervals for the estimated strength of persistence. We finish this paper by quantifying long-range persistence (and corresponding uncertainties) of three geophysical time series—palaeotemperature, river discharge, and Auroral electrojet index—with the three representing three different types of probability distribution—Gaussian, log-normal, and Levy, respectively. 相似文献
17.
Claudia Beatriz Quinteros Cartaya Fidencio Alejandro Nava Pichardo Ewa Glowacka José Duglas Frez Cárdenas 《Pure and Applied Geophysics》2014,171(7):1367-1383
In order to analyze observed seismicity in central Japan and Venezuela, we applied a new method to identify semi-periodic sequences in the occurrence times of large earthquakes, which allows for the presence of multiple periodic sequences and/or events not belonging to any sequence in the time series. We also explored a scheme for diminishing the effects of a sharp cutoff magnitude threshold in selecting the events to analyze. A main four-event sequence with probability P c = 0.991 of not having occurred by chance was identified for earthquakes with M ≥ 8.0 in central Japan. Venezuela is divided, from West to East, into four regions; for each of these, the magnitude ranges and identified sequences are as follows. Region 1: M ≥ 6.0, a six-event sequence with P c = 0.923, and a four-event sequence with P c = 0.706. Region 2: M ≥ 5.6, a five-event sequence with P c = 0.942. Region 3: M ≥ 5.6, a four-event sequence with P c = 0.882. Region 4: M ≥ 6.0, a five-event sequence with P c = 0.891. Forecasts are made and evaluated for all identified sequences having four or more events and probabilities ≥0.5. The last event of all these sequences was satisfactorily aftcast by previous events. Whether the identified sequences do, in fact, correspond to physical processes resulting in semi-periodic seismicity is, of course, an open question; but the forecasts, properly used, may be useful as a factor in seismic hazard estimation. 相似文献
18.
Pierre Polsenaere Nicolas Savoye Henri Etcheber Mathieu Canton Dominique Poirier Steven Bouillon Gwenaël Abril 《Aquatic Sciences - Research Across Boundaries》2013,75(2):299-319
We measured spatial and temporal variations in carbon concentrations, isotopic compositions and exports during a complete hydrological cycle in nine watercourses draining a lowland forested podzolized catchment, flowing into the Arcachon lagoon (France). In addition, integrated fluxes of CO2 across the water-atmosphere interface were estimated to assess the relative importance of CO2 evasion versus lateral carbon transport at the catchment scale. Watercourse similarities and specificities linked to the local catchment characteristics are discussed and compared with other riverine systems. Low concentrations of suspended particulate matter and particulate organic carbon (POC) were generally measured in all the watercourses (8.4 ± 3.4 and 1.6 ± 0.6 mg L?1, respectively), reflecting limited mechanical soil erosion. The generally high POC content in the suspended matter (20 %), low Chl a concentrations (1.3 ± 1.4 μg L?1) and the relatively constant δ13C-POC value (near ?28 ‰) throughout the year reveal this POC originates from terrestrial C3 plant and soil detritus. The presence of podzols leads to high levels of dissolved organic carbon (DOC; 6.6 ± 2.2 mg L?1). Similarly, high dissolved inorganic carbon (DIC) concentrations were measured in the Arcachon lagoon catchment (5.9 ± 2.2 mg L?1). The δ13C-DIC value around ?20 ‰ throughout the year in many small watercourses reveals the predominance of terrestrial carbon mineralisation and silicate rock weathering in soils as the major DIC source. With pCO2 between 1,000 and 10,000 ppmv, all watercourses were a source of CO2 to the atmosphere, particularly during the low river stage. Organic carbon parameters remained relatively stable throughout the year, whereas DIC parameters showed strong seasonal contrasts closely linked to the hydrological regime and hyporheic flows. In total, the carbon export from the Arcachon watershed was estimated at 15,870 t C year?1 or 6 t C km?2 year?1, mostly exported to the lagoon as DOC (35 %), DIC (24 %) and lost as CO2 degassing to the atmosphere (34 %). 相似文献
19.
During geothermal power production using a borehole doublet consisting of a production and injection well, the reservoir conditions such as permeability k, porosity φ and Skempton coefficient B at the geothermal research site Gross Schoenebeck/Germany will change. Besides a temperature decrease at the injection well and a change of the chemical equilibrium, also the pore pressure p p will vary in a range of approximately 44 MPa ± 10 MPa in our reservoir at ?3850 to ?4258 m depth. This leads to a poroelastic response of the reservoir rocks depending on effective pressure p eff (difference between mean stress and pore pressure), resulting in a change in permeability k, porosity φ and the poroelastic parameter Skempton coefficient B. Hence, we investigated the effective pressure dependency of Flechtinger sandstone, an outcropping equivalent of the reservoir rock via laboratory experiments. The permeability decreased by 21% at an effective pressure range from 3 to 30 MPa, the porosity decreased by 11% (p eff = 6 to 65 MPa) and the Skempton coefficient decreased by 24% (p eff = 4 to 25 MPa). We will show which mechanisms lead to the change of the mentioned hydraulic and poroelastic parameters and the influence of these changes on the productivity of the reservoir. The most significant changes occur at low effective pressures until 15 to 20 MPa. For our in situ reservoir conditions p eff = 43 MPa a change of 10 MPa effective pressure will result in a change in matrix permeability of less than 4% and in matrix porosity of less than 2%. Besides natural fracture systems, fault zones and induced hydraulic fractures, the rock matrix its only one part of geothermal systems. All components can be influenced by pressure, temperature and chemical reactions. Therefore, the determined small poroelastic response of rock matrix does not significantly influence the sustainability of the geothermal reservoir. 相似文献