首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 818 毫秒
1.
NMR relaxometry is a powerful tool for inferring porosity and permeability data. In practice, measured magnetization decay curves are inverted for relaxation time distributions. Subsequently, one presumes a linear relationship between the pore radius distribution and the T 1 and T 2 distribution, for longitudinal and transverse magnetization, respectively. The fundamental equations used are based on a pore model, in which pores are assumed to be isolated from each other with respect to the NMR process and have smooth walls. The present study is based on a geometrical pore space model with connected pores and structured pore walls. The physical processes of surface relaxation, irreversible dephasing of magnetic spins and diffusive proton exchange between pores, are described by a system of differential equations. The solution yields a set of exponential functions representing the relaxation time distribution. We describe the difference between the distributions obtained for diffusion coupling and for isolated pores. With diffusion coupling on, the spectral width of the T 1 distribution is strongly reduced, which indicates that the influence of large and small radii according to the T 1-pore radius relationship is mixed to some extent. For a fractal pore space structure, where large pores are surrounded by adjacent minor pores, the T 1 distribution does not resolve these substructures. Nevertheless, permeability values calculated from the logarithmic mean relaxation time T 1,LM are quite the same for diffusion coupling and for isolated pores. The T 2 distribution for diffusion coupling is little constricted and gives a better resolution of the pore wall structures than the corresponding T 1 distribution. The permeability values from T 2 distributions agree with the values from longitudinal magnetization, provided that we use a corrected relaxation time T 2,corr, accounting for the dependence of the surface relaxivity ρ 2 on pore radius. The study shows that radius distributions calculated from a T 1 and from a T 2 distribution differ from one another and both present an altered image of the true pore radius distribution. In practice, this has no serious influence on estimating permeability of medium- to high-permeability sandstones with the currently applied methods. The presented methodology of calculating the NMR response of pore space models with diffusion coupling may facilitate understanding porosity-permeability relationships of different rock types such as carbonate rocks with micro-porosity.  相似文献   

2.
核磁共振测井探测岩石内部磁场梯度的方法   总被引:3,自引:1,他引:2       下载免费PDF全文
岩石内部磁场梯度对核磁共振测井横向弛豫T2分布有较大的影响.本文提出了一种利用自旋回波脉冲序列探测岩石内部磁场梯度的二维核磁共振方法,发展了相应的反演方法,得到了含顺磁性物质的饱和水砂岩和泥质砂岩的内部磁场梯度分布规律.结果表明,随顺磁物质含量的增加,岩石内部磁场梯度增大.含顺磁物质的饱和水砂岩颗粒均匀,不含粘土矿物,内部磁场梯度呈单峰分布特征.含绿泥石粘土矿物的饱和水泥质砂岩,T2分布呈双峰特征,束缚水峰明显,小孔短T2对应的梯度值大于大孔长T2对应的梯度值.当绿泥石含量小于15%时,随绿泥石含量的增加,对应束缚水峰的短T2组分明显增多,T2分布展宽;当绿泥石含量大于15%以后,较大的岩石内部磁场梯度使短弛豫信息衰减迅速,表现出随绿泥石含量的增加,短T2组分减少,T2分布变窄.  相似文献   

3.
内部磁场梯度引起的扩散对NMR岩石测量响应的影响   总被引:2,自引:2,他引:0       下载免费PDF全文
低场NMR岩心分析能够刻度NMR测井响应,对于较为准确地预测储层的渗透率、束缚水体积等与产能密切相关的参数尤其重要.对来自南海东部油田的100%饱和盐水的砂岩岩心进行了变回波间隔的实验室NMR T2测量,随着回波间隔的增大对实验观测到的两种不同的T2 分布的移动进行了理论上的分析和解释.随着回波间隔的增加,T2谱向着长弛豫时间的移动可以用过优化NMR采集参数消除掉,而对另外一种由内部磁场梯度引起的移动即随着回波间隔增加T2谱向短弛豫时间的移动则复杂得多,至今为止也无法做到定量化内部磁场梯度的值,这种移动能够引起错误的NMR测井解释.因此理解内部磁场梯度对T2弛豫时间的影响是很有必要的.文中对由内部磁场梯度扩散引起的扩散弛豫对T2弛豫时间的影响进行了理论上的模拟计算,这有助于理解和解释岩石内部磁场梯度对NMR T2弛豫时间的影响.最后结合压汞毛管压力曲线,解释了具体的岩心实验结果,并且计算了具体的岩心的内部磁场梯度值,计算出的内部磁场梯度值应该被视为内部磁场梯度的几何平均值.  相似文献   

4.
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time (T 1) and transverse relaxation time (T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times (T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings (T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.  相似文献   

5.
The pore structure of the tight limestone in the Daanzhai Member of the Ziliujing Formation, Jurassic System, in central Sichuan Basin, China, is complex but essential to the exploration and development of tight oil. The pore structure of the tight limestone is studied by using scanning electron microscopy (SEM), nitrogen adsorption, high-pressure mercury intrusion, and nuclear magnetic resonance (NMR). The experimental results suggest that the pores are mainly slit pores and mesopores and macropores contribute to the pore volume and specific surface. The displacement pressure, average pore size, and homogeneity coefficient correlate with porosity and permeability and can be used to evaluate the pore structure. The full pore-size distribution was obtained by combining nitrogen adsorption and high-pressure mercury intrusion. We find that the limestone mainly contains mesopores with diameter of 2?50 nm. The T2 distribution was converted into pore-size distribution, well matching the full pore-size distribution. The relation between T2 and pore size obeys a power law and the geometric mean of T2 correlates with the pore structure and can be used in the pore structure evaluation.  相似文献   

6.
Modeling rock permeability from NMR relaxation data by PLS regression   总被引:1,自引:0,他引:1  
This study explores the application of the partial least squares regression (PLSR) technique to rock permeability prediction from nuclear magnetic resonance (NMR) relaxation data. A total of 68 Brazilian sandstone cores selected from reservoirs and outcrop analogs were fully saturated and analyzed by NMR. The permeability of the cores ranged from 0.007 to 9,800 mD. From their 1H transverse relaxation times (T2) measured at 2 MHz, two PLSR models were developed for the relaxation spectra and the raw relaxation curves. Both models led to more uniform and accurate predictions (RMSE = 0.47 and 0.50 log mD, respectively) compared with the classical Kenyon model (RMSE = 0.78 log mD).  相似文献   

7.
核磁共振共轭梯度解谱方法研究   总被引:1,自引:1,他引:0       下载免费PDF全文
核磁共振(NMR)解谱技术是核磁共振资料应用的基础和关键.本文将核磁共振解谱的混定线性反演问题转化为求目标函数极小值的最优化问题,然后利用共轭梯度算法具有二次终止性、收敛速度快的特点解决上述最优化问题.将该方法应用于无噪声和不同信噪比的理论数据解谱以及岩心NMR实验解谱并与真谱及实验室国外软件解谱结果对比表明:在信噪比SNR≥5时解谱结果和真谱符合得很好,谱线光滑连续;随着信噪比的降低对初始点的要求随之提高;两块岩心解谱结果与实验室结果符合得很好,利用解谱结果计算的核磁孔隙度与实验室氦孔隙度绝对误差分别为0.78%和0.57%.因此,本文方法有效、实用,具有较强的抗噪能力,对初始点的要求不高,能够应用于生产和科研实践中.  相似文献   

8.
Based on the analysis of mercury injection capillary pressure (MICP) and nuclear magnetic resonance (NMR) experimental data for core plugs, which were drilled from two Chinese tight sandstone reservoirs, permeability prediction models, such as the classical SDR, Timur-Coates, the Swanson parameter, the Capillary Parachor, the R10 and R35 models, are calibrated to estimating permeabilities from field NMR logs, and the applicabilities of these permeability prediction models are compared. The processing results of several field examples show that the SDR model is unavailable in tight sandstone reservoirs. The Timur-Coates model is effective once the optimal T 2cutoff can be acquired to accurately calculate FFI and BVI from field NMR logs. The Swanson parameter model and the Capillary Parachor model are not always available in tight sandstone reservoirs. The R35 based model cannot effectively work in tight sandstone reservoirs, while the R10 based model is optimal in permeability prediction.  相似文献   

9.
Several parameters of pore geometry are needed for estimating permeability which is a key parameter for the characterization of reservoir sandstones. Powerful techniques for probing the pore space are the self-diffusion and the relaxation time NMR methods. However, the quality of results depends on the petrophysical model which underlies the interpretation of measurements. We applied the pulsed field gradient nuclear magnetic resonance (PFG-NMR) technique and measured time-dependent self-diffusion coefficients, D(Δ), of water in anhydrite cemented sandstones with low porosity and high tortuosity. The conventional method of fitting data with a function (Padé approximation) of the surface-to-volume ratio and the tortuosity yielded uncertain results. As part of a novel approach, we developed a numerical simulation code based on physical principles and a fractal pore space model. We compared our method with the Padé approximation and tested with data from the literature. For porous media with low tortuosity values and simple geometry such as randomly packed glass beads, both methods are in good agreement and give similar results. In sedimentary rocks, however, the new method of calculation is able to determine more accurate details of the pore geometry. However, its main advantage occurs in porous media with fractal geometry or with constricted pores, where the D(Δ)/D0 curves deviate from the characteristic shape that is found for randomly packed glass beads.  相似文献   

10.
Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T2 data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T2 could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.  相似文献   

11.
储层渗透性与地层因素关系的实验研究与分析   总被引:3,自引:0,他引:3       下载免费PDF全文
本文对渤海湾盆地不同孔隙结构样品的孔、渗、核磁、岩电、压汞、X衍射及铸体薄片等配套岩石物理实验数据进行了综合分析,通过逐一考察同一套岩芯样品的地层因素与渗透率、压汞喉径均值、储层品质指数之间的实验关系,并分别与地层因素-孔隙度交绘图进行对比分析,发现储层渗透性与地层因素之间并非简单的单调函数关系,孔隙度相近但孔隙结构类型不同、渗透率差异明显的岩芯可以具有相近的地层因素,导电能力接近.在实验数据分析的基础上通过理论分析证明了这一实验关系的合理性,并指出孔隙度及导电能力相近的岩芯,其渗透率差异与喉径均值的平方比、孔隙曲折度及几何形态相关.  相似文献   

12.
The self-diffusion of water and hexadecane in medium and coarse sands from glacial sand deposits in central Germany were investigated by pulsed field gradient nuclear magnetic resonance (PFG NMR). Due to the restriction of the diffusion path at the pore/grain interface, the measured apparent self-diffusion coefficients (D(Δ)) in the pore space depend on the observation time (Δ) in the PFG NMR experiment. Although the bulk self-diffusion coefficients of water and hexadecane differ by about one order of magnitude, the apparent self-diffusion coefficients in the pore space obey the same characteristic time-behaviour, which depends only on geometrical properties of the pore system. Using the “short-time diffusion” model, surface-to-volume (S/V) ratios and inherent self-diffusion coefficients (D0) of the pore fluids were extracted from these diffusion measurements. The S/V ratios obtained are independent of the pore fluid used and agree with known geometrical properties of the sand grains. Moreover, the D0 values are consistent with the corresponding bulk self-diffusion coefficients measured separately. In contrast to these results of PFG NMR, simultaneous investigations of longitudinal (T1) nuclear magnetic relaxation reveal that the relaxation time of the pore fluid is a less suitable parameter for a quantitative estimation of geometrical properties of the pore/grain interface in these unconsolidated sediments since it depends on chemical properties of the fluid/grain interface.  相似文献   

13.
We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO2 = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO2 Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO2 used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8).  相似文献   

14.
Pore structure characteristics are important to oil and gas exploration in complex low-permeability reservoirs. Using multifractal theory and nuclear magnetic resonance (NMR), we studied the pore structure of low-permeability sandstone rocks from the 4th Member (ES4) of the Shahejie Formation in the south slope of the Dongying Sag. We used the existing pore structure data from petrophysics, core slices, and mercury injection tests to classify the pore structure into three categories and five subcategories. Then, the T2 spectra of samples with different pore structures were interpolated, and the one- and three-dimensional fractal dimensions and the multifractal spectrum were obtained. Parameters α (intensity of singularity) and f (α) (density of distribution) were extracted from the multifractal spectra. The differences in the three fractal dimensions suggest that the pore structure types correlate with α and f (α). The results calculated based on the multifractal spectrum is consistent with that of the core slices and mercury injection. Finally, the proposed method was applied to an actual logging profile to evaluate the pore structure of low-permeability sandstone reservoirs.  相似文献   

15.
A methodology for rock classification is presented that considers lithology and reservoir parameters on the basis of a combined fractal analysis of well logs and mercury porosimetry results with nuclear magnetic resonance outcomes. A sandy‐shaly thinly‐bedded Miocene gas bearing formation in the Carpathian Foredeep is investigated. Fractal correlation dimensions D2 calculated for standard logs are used to distinguish sandstone as the most homogeneous lithological group with the highest porosity. The fractal analysis also confirmed observations of gas accumulations in sandstone, shaly sandstone and sandy claystone as rocks of high porosity. The results of two laboratory methods are combined to improve reservoir properties assessment and evaluate movable media in pore space; this technique was based on the similarity of mercury porosimetry results plotted as cumulative intrusion volume versus pressure or pore diameter and also the curves of cumulative porosity and transverse relaxation time distributions using nuclear magnetic resonance. Close values of porosity from logs, recorded in situ and in laboratory measurements, provide the link between fractal analysis and porosimetry and nuclear magnetic resonance measurements.  相似文献   

16.
NMR Measurements in Carbonate Rocks: Problems and an Approach to a Solution   总被引:2,自引:0,他引:2  
Carbonate rocks are well known for their complex petrophysical behavior where, in contrast to siliciclastic rocks, different parameters, including porosity and permeability, usually are not directly related. This behavior is the result of thorough reorganization of porosity during diagenesis, and it turns prediction of reservoir quality of carbonate rocks into a challenge. The study presented here deals with the problem of utilizing NMR techniques in prediction of petrophysical properties in carbonates.We employ a visual porosity classification as a priori knowledge for better interpreting NMR data for prediction purposes. This allows for choice of suitable T2 cutoff values to differentiate movable from bound fluids adapted for the specific carbonate rock, thus resulting in better interpretation of NMR data. The approach of using a genetic pore type classification for adapting the conventional method for T2 cutoff determination, which originally was developed for siliciclastic rocks, is promising. Similarly, for permeability determination on the basis of NMR measurements, the classification of carbonate rocks based on porosity types also shows potential. The approach implemented here has the promise to provide a basis of standardized interpretation of NMR data from carbonate rocks.Acknowledgment We are grateful to Baker Hughes INTEQ, Celle, for permission to publish the results of this work. This study was part of I.S.s Masters Thesis at Hannover University that was carried out in collaboration with Baker Hughes INTEQ, Celle, Germany. We appreciate comments of an anonymous referee and PAGEOPH editor Brian J. Mitchell.  相似文献   

17.
Two time scales are distinguished in the geomagnetotail dynamics. The small scale (T 1) corresponds to disturbances propagating in the tail lobes, which have a relatively strong magnetic field and low plasma density. The larger scale (T 2) corresponds to plasma motions in the plasma sheet and has a relatively weak magnetic field and a relatively higher density. A disturbance, which is initiated by a localized burst of magnetic reconnection and appears in the geomagnetotail on the time scale T 1, generates the upset of equilibrium in the plasma sheet zones with intermediate spatial dimensions (about R E). The theoretical considerations and numerical simulation indicate that the relaxation process, which subsequently proceeds on the larger time scale (T 2), results in the appearance of extremely thin embedded current sheets and in the generation of fast plasma flows. This process gives an effective mechanism by which the magnetic energy stored in the geomagnetotail is transformed into the plasma flow kinetic energy. Such fast flows can also generate eddy plasma motions on smaller spatial scales. On the one hand, fast MHD components of this process carry a disturbance in other plasma sheet zones, where new magnetic reconnection bursts can originate at a large distance from the zone of an initial nonlinear disturbance. As a result, new recurrent processes of relaxation originate on the T 2 time scale. Alternation originating in such a way is apparently the characteristic feature of eddy disturbances actually observed in the plasma sheet.  相似文献   

18.
Using calculations of the magnetic field in the solar atmosphere in the potential approximation, it is shown that, (1) as distance R from the Sun’s center grows, the area of the positive magnetic field (S +field) in 10-deg latitude zones tends to 100% (0%) in the neighborhood of the solar minimum. At the distance R = 2.5R (R is the solar radius), these values of the positive field are observed during ≈(12–55) Carrington rotations (CRs) for solar minima between neighboring cycles; (2) polar magnetic field reversals can occur repeatedly. Note that a polar reversal at large heights ends by 6–16 Carrington rotations earlier than on the Sun’s surface. On the Sun’s surface, a field polar reversal begins earlier at lower latitudes than at high ones; (3) for each longitude at different Rs and separately for each solar hemisphere the radial component of the field was averaged on synoptic maps in the 0°–40° latitude range. It is established that the T R rotation periods of the boundaries between the sectors (areas of longitudes with the same sign of the averaged field) can be shorter than, longer than, and equal to Carrington solar rotation period T CR. It turned out that boundaries with T R < T CR are observed at all heights, while boundaries with T R > T CR are observed at relatively small heights.  相似文献   

19.
用地面核磁共振方法评估含水层涌水量的实例   总被引:7,自引:4,他引:3       下载免费PDF全文
地面核磁共振(NMR)方法是地球物理上采用的探测地下水的最新方法,能够探测地下含水层中的自由水.但是有时会出现对地下含水层的出水量判断失误的现象,这种失误发生的主要原因是由于IRIS仪器设计时假设利用NMR信号的弛豫时问就能够区分地下含水层中的自由水和束缚水,实际上弛豫时间取决于以下几个参数即:孔隙度、渗透率和导水率.地层孔隙水中氢质子弛豫时间不仅与其本身的弛豫特点有关,还与岩石孔隙结构、成份密切相关.为了更好的理解以上这些参数对地下含水层涌水量的影响,本文讨论了含水多孔介质的弛豫特性,研究有效孔隙度与含水量之间的关系,给出计算渗透率和导水率的方法,利用几个实测地点数据资料分析地下含水层岩性对涌水量的影响,结果表明地层中平均含水量大而且含水地层弛豫时间较长的地点才能获得较大的涌水量.  相似文献   

20.
Contaminant rebound and low contaminant removal are reported more frequently with in situ chemical oxidation than other in situ technologies. Although there are multiple causes for these results, a critical analysis indicates that low oxidant volume delivery is a key issue. The volume of oxidant injected is critical and porosity of the aquifer matrix can be used to estimate the pore volume. The total porosity (qT) is the volume of voids relative to the total volume of aquifer material. The mobile porosity (qM) is the fraction of voids that readily contributes to fluid displacement, and is less than qT leading to smaller estimates of oxidant volume. Injecting low‐oxidant volume may result in inadequate oxidant distribution and postinjection dispersal within the radius of influence, insufficient oxidant contact and oxidant loading, and incomplete treatment; whereas, greater oxidant volume achieves a greater oxidant footprint and may involve risk that the injected oxidant may migrate into nontarget areas and displacement of contaminated groundwater. Design guidelines and recommendations are provided that could help achieve more effective technology deployment, reduce the role of heterogeneities in the subsurface, and result in greater probability the oxidant is delivered to the targeted treatment zone.  相似文献   

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

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