首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
限制性克立格法在地学中的研究与应用   总被引:4,自引:2,他引:2       下载免费PDF全文
普通克立格法在矿产资源储量计算中通常会产生较多的平滑效应,本文运用限制性克立格法,通过其限制性条件来解决这一问题.阐述了限制性克里格的理论方法,设计出限制性克立格法的计算过程,结合实例,通过块段模型的结果对普通克里格法与限制性克里格法进行比较,结果表明,限制性克立格法在抑制平滑效应方面能够取得较为理想的结果.  相似文献   

2.
地应力对裂隙岩体渗透特性的影响   总被引:9,自引:0,他引:9  
提出了地应力场中岩体不连续面变形的本构关系,讨论了渗流与应力耦合条件下裂隙岩体渗透张量的计算方法,阐述了地应力对裂隙岩体渗透特性的影响.研究表明,地应力对岩体渗透特性的影响主要是通过改变不连续面的开度和不连续面网络的水力传导路径而产生的;裂隙岩体渗透性随岩体埋深的增加呈负指数减小,随侧压系数的增大呈双曲线减小.渗流场与应力场的耦合理论可望应用于水库诱发地震的研究中  相似文献   

3.
地应力对裂隙岩体渗透特性的影响   总被引:6,自引:0,他引:6       下载免费PDF全文
周创兵  熊文林 《地震学报》1997,19(2):154-163
提出了地应力场中岩体不连续面变形的本构关系,讨论了渗流与应力耦合条件下裂隙岩体渗透张量的计算方法,阐述了地应力对裂隙岩体渗透特性的影响.研究表明,地应力对岩体渗透特性的影响主要是通过改变不连续面的开度和不连续面网络的水力传导路径而产生的;裂隙岩体渗透性随岩体埋深的增加呈负指数减小,随侧压系数的增大呈双曲线减小.渗流场与应力场的耦合理论可望应用于水库诱发地震的研究中.   相似文献   

4.
韩家岩体与天井山金矿空间上密切相关.对韩家岩体的成因存在不同看法:依据“斑状”外貌,有人主张是燕山期花岗斑岩;依据隐约可见的片麻状构造,有人认为是晋宁期灵山岩体.鉴于其具有重要的找矿指示意义,本文通过岩石学、岩石地球化学、同位素年代学等方面的研究,认为韩家岩体是由灵山岩体蚀变而成的绢英岩.天井山金矿中的金多以自然金形式产于石英大脉的晚期裂隙中,并与石英-绢云母-含铁碳酸盐-硫化物细脉共生.研究表明,这些裂隙与绢英岩体顶部的石英网脉带贯通,两者具有相同蚀变矿物组合和流体来源.据此,本文认为与绢英岩化有关的岩浆-构造-蚀变作用是天井山金矿的核心成矿事件.由于韩家绢英岩体具有区域尺度的找矿标志,这对天井山矿区及其外围进一步找矿具有重要的指导意义.  相似文献   

5.
针对目前高坝坝址地震动参数确定方法的一些不足,本文尝试了一种基于地震学理论的坝址区三维地震动场生成方法.基本思路是将地震学和工程学结合,针对设定地震,建立震源-传播介质-坝址峡谷场地数值模型,通过超大规模的数值计算,模拟地震波从发震断层破裂开始到坝址场地的物理传播过程,生成坝址区的三维地震动参数.与传统基于衰减关系的地震动参数确定方法相比,这一方法可以考虑震源机制、传播介质和坝址峡谷场地效应等三大要素的影响.对于特定的坝址,可以生成符合实际地质构造、区域岩体动力特性以及坝址峡谷地质地形条件的地震动荷载分布,具有针对性,为重大高坝枢纽工程遭遇极端地震荷载作用时的抗震安全性提供了一种新的分析手段.  相似文献   

6.
裂隙储层具有强烈的非均质性及储层渗透性预测难度大的特点,利用逾渗理论来研究裂隙储层深层复杂介质渗透性,是当前研究的热点问题.基于连续逾渗,将裂隙网络模型合理简化.使用排除体积对裂隙密度进行无量纲化,从而使裂隙的渗透性与裂隙的形状无关.首先使用Monte Carlo方法得到裂隙在不同密度时的网络图,而后使用ComsolMultiphysics中的有限元求解器得到流体在裂隙中的速度和压力分布图,数值模拟表明:随着裂隙密度的增加,裂隙的连通性增加,进而流体在裂隙中的渗透性也增加;最后使用有限尺度转换定律和尺度放大思想,通过重复子区域的方法求解每个小区域的渗透性,得到裂隙储层在不同尺度下的渗透性,分析得到储层平均渗透性和无量纲化裂隙密度之间的关系.模拟得到的渗透率值和FRACA软件得到渗透率泊松相关系数为0.885.该规则的发现为逾渗理论从小尺度放大到宏观物理模型即储层尺度提供了一个理论依据.  相似文献   

7.
增强型地热系统(Enhanced Geothermal System,简称EGS)作为深层地热资源开采最有效方法已成为国际研究热点.其中多场耦合是EGS研究中的关键部分,但目前对最基础和最重要的热流耦合研究仍不充分,特别是对热流耦合下热传导过程和采热效率的研究还不透彻.本文构建了针对EGS单裂隙热流耦合的数学模型,模型主要控制热流耦合时水岩温度场的演化,以展现热流耦合下的热传导过程,并以此计算了裂隙宽度和裂隙水流速率对水岩温度场的影响,分析了裂隙特征对采热效率的影响.计算表明,在特定的开采条件下,都有相应的岩体温度场影响半径;热量从高温岩体传至裂隙流体的过程显示为进水口为低温区,以水流方向为轴,向两侧岩体对称展开的三角形热传导模式;裂隙宽度和裂隙流体速度对裂隙两侧的岩体温度场和采热半径产生明显影响,但裂隙长度总体上对采热率表现为正效应,而裂隙宽度和裂隙流体速度在一定范围内对采热率表现为负效应;裂隙宽度与裂隙流体速度对水岩温度场具有完全一致的影响效果.模型的对比验证确认了本模型的合理性和优越性.文中最大的创新点和亮点是首次推导给出了EGS热流耦合水岩温度场演化特征的简明数学表达式,另一亮点是基于岩体热传导系数来考虑耦合过程的热传导量.  相似文献   

8.
通过实地勘测研究,建立了黄河黑山峡大柳树坝址区岩体的滑移-弯曲拉裂松动变形模式。并运用2D-Sig-ma有限元程序对其进行分析,充分论证了坝址区岩体的顺层滑移-弯曲拉裂是该处岩体松动的一种重要模式,并给出了若干重要结论。  相似文献   

9.
采用声波测井、地震CT等综合物探方法,对沙坡头水利枢纽坝址基础岩体固结灌浆效果进行检测,基本查明了坝基岩体固结灌浆后的质量状况,取得了良好的应用效果.  相似文献   

10.
冰架崩解是南极洲物质损耗的主要途径,崩解是应力作用引起的裂隙(或裂缝)传播的结果,裂隙位置和深度的探测是理解崩解机理过程的重要基础.本文提出了一种利用卫星激光测高ICESat-1/GLAS高程数据产品提取冰架表面冰裂隙的方法,并以南极洲埃默里冰架为例验证了这种方法探测裂隙位置的准确性和深度探测精度.同时,基于提取的裂隙点深度分布特征提出了裂隙峰值应力点的探测方法,可用于追踪冰裂隙初始裂口位置和探测导致冰架崩解的高危区.利用2003~2008年间16个运行时期内132条ICESat-1/GLAS高程轨迹线分析了埃默里冰架冰裂隙深度的时空分布和裂隙峰值应力点的空间分布.结果显示,探测到的裂隙点深度在2.0~31.7 m均在海平面以上;裂隙深度变化未显示出随时间推移和冰流移动而增加的趋势,说明平流移动到冰架前缘的裂隙基本不会直接导致冰架的崩解;冰架局部应力集中区主要分布在冰流的缝合区内.  相似文献   

11.
针对裂隙岩体提出了一种基于双重逾渗模型的渗流系数分析方法。通过变参数研究,考察了岩体中渗流和孔压分布的不均匀;根据裂隙长度的不同将岩体分为孔隙控制渗流模式和裂隙控制渗流模式。这种模型不仅能够同时考虑孔隙和裂隙的作用,而且在模拟大规模裂隙网络时更加快捷简便。  相似文献   

12.
It is common in geostatistics to use the variogram to describe the spatial dependence structure and to use kriging as the spatial prediction methodology. Both methods are sensitive to outlying observations and are strongly influenced by the marginal distribution of the underlying random field. Hence, they lead to unreliable results when applied to extreme value or multimodal data. As an alternative to traditional spatial modeling and interpolation we consider the use of copula functions. This paper extends existing copula-based geostatistical models. We show how location dependent covariates e.g. a spatial trend can be accounted for in spatial copula models. Furthermore, we introduce geostatistical copula-based models that are able to deal with random fields having discrete marginal distributions. We propose three different copula-based spatial interpolation methods. By exploiting the relationship between bivariate copulas and indicator covariances, we present indicator kriging and disjunctive kriging. As a second method we present simple kriging of the rank-transformed data. The third method is a plug-in prediction and generalizes the frequently applied trans-Gaussian kriging. Finally, we report on the results obtained for the so-called Helicopter data set which contains extreme radioactivity measurements.  相似文献   

13.
In humid, well-vegetated areas, such as in the northeastern US, runoff is most commonly generated from relatively small portions of the landscape becoming completely saturated, however, little is known about the spatial and temporal behavior of these saturated regions. Indicator kriging provides a way to use traditional water table data to quantify probability of saturation to evaluate predicted spatial distributions of runoff generation risk, especially for the new generation of water quality models incorporating saturation excess runoff theory. When spatial measurements of a variable are transformed to binary indicators (i.e., 1 if above a given threshold value and 0 if below) and the resulting indicator semivariogram is modeled, indicator kriging produces the probability of the measured variable to exceed the threshold value. Indicator kriging gives quantified probability of saturation or, consistent with saturation excess runoff theory, runoff generation risk with depth to water table as the variable and the threshold set near the soil surface. The probability of saturation for a 120 m × 180 m hillslope based upon 43 measurements of depth to water table is investigated with indicator semivariograms for six storm events. The indicator semivariograms show high spatial structure in saturated regions with large antecedent rainfall conditions. The temporal structure of the data is used to generate interpolated (soft) data to supplement measured (hard) data. This improved the spatial structure of the indicator semivariograms for lower antecedent rainfall conditions. Probability of saturation was evaluated through indicator kriging incorporating soft data showing, based on this preliminary study, highly connected regions of saturation as expected for the wet season (April through May) in the Catskill Mountain region of New York State. Supplementation of hard data with soft data incorporates physical hydrology of the hillslope to capture significant patterns not available when using hard data alone for indicator kriging. With the need for water quality models incorporating appropriate runoff generation risk estimates on the rise, this manner of data will lay the groundwork for future model evaluation and development.  相似文献   

14.
岩石渗透特性的实验研究   总被引:2,自引:0,他引:2       下载免费PDF全文
本文介绍了在低压下(最大压力20MPa)对完整岩石和含自然裂面岩石渗透率的实验研究。在实验中对四种材料(花岗岩、变质岩、大理岩、水泥块)用六种方法进行了测试。实验结果表明,低压下完整岩石的渗透率为1μd量级,含裂面岩石的渗透率为1md量级,渗透率明显地与应力状态和加载方式有关。  相似文献   

15.
Geostatistical methods are well suited for analyzing the local and spatial uncertainties that accompany the modeling of highly heterogeneous three-dimensional (3D) geological architectures. The spatial modeling of 3D hydrogeological architectures is crucial for polluted site characterization, in regards to both groundwater modeling and planning remediation procedures. From this perspective, the polluted site of Porto Marghera, located on the periphery of the Venice lagoon, represents an interesting example. For this site, the available dense spatial sampling network, with 769 boreholes over an area of 6 km2, allows us to evaluate the high geological heterogeneity by means of indicator kriging and sequential indicator simulation. We show that geostatistical methodologies and ad hoc post processing of geostatistical analysis results allow us to effectively analyze the high hydrogeological heterogeneity of the studied site.  相似文献   

16.
Compositional Bayesian indicator estimation   总被引:1,自引:1,他引:0  
Indicator kriging is widely used for mapping spatial binary variables and for estimating the global and local spatial distributions of variables in geosciences. For continuous random variables, indicator kriging gives an estimate of the cumulative distribution function, for a given threshold, which is then the estimate of a probability. Like any other kriging procedure, indicator kriging provides an estimation variance that, although not often used in applications, should be taken into account as it assesses the uncertainty of the estimate. An alternative approach to indicator estimation is proposed in this paper. In this alternative approach the complete probability density function of the indicator estimate is evaluated. The procedure is described in a Bayesian framework, using a multivariate Gaussian likelihood and an a priori distribution which are both combined according to Bayes theorem in order to obtain a posterior distribution for the indicator estimate. From this posterior distribution, point estimates, interval estimates and uncertainty measures can be obtained. Among the point estimates, the median of the posterior distribution is the maximum entropy estimate because there is a fifty-fifty chance of the unknown value of the estimate being larger or smaller than the median; that is, there is maximum uncertainty in the choice between two alternatives. Thus in some sense, the latter is an indicator estimator, alternative to the kriging estimator, that includes its own uncertainty. On the other hand, the mode of the posterior distribution estimator, assuming a uniform prior, is coincidental with the simple kriging estimator. Additionally, because the indicator estimate can be considered as a two-part composition which domain of definition is the simplex, the method is extended to compositional Bayesian indicator estimation. Bayesian indicator estimation and compositional Bayesian indicator estimation are illustrated with an environmental case study in which the probability of the content of a geochemical element in soil being over a particular threshold is of interest. The computer codes and its user guides are public domain and freely available.  相似文献   

17.
The response of two arch dams to spatially varying ground motions recorded during earthquakes is computed by a recently developed linear analysis procedure, which includes dam–water–foundation rock interaction effects and recognizes the semi‐unbounded extent of the rock and impounded water domains. By comparing the computed and recorded responses, several issues that arise in analysis of arch dams are investigated. It is also demonstrated that spatial variations in ground motion, typically ignored in engineering practice, can have profound influence on the earthquake‐induced stresses in the dam. This influence obviously depends on the degree to which ground motion varies spatially along the dam–rock interface. Thus, for the same dam, this influence could differ from one earthquake to the next, depending on the epicenter location and the focal depth of the earthquake relative to the dam site. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

18.
Zhang F  Zhou Z  Huang Y  Chen Z 《Ground water》2004,42(4):509-515
A typical fractured rock mass is intersected by several sets of discontinuities, which provide the main flowpath for ground water. Due to the limitations of data obtained by conventional field measurements, it is often difficult to estimate the anisotropic permeability tensor associated with the joints existing in the rock mass. For that reason, determining permeability tensors for fractured rocks is an important topic in rock mass hydraulics. Based on field surveys, joint parameters can be analyzed by using probabilistic and statistical tools, and three-dimensional mapping of the jointed rock mass. Through analysis of a single joint's hydraulic characteristics, the principal value of the permeability tensor for the jointed rock mass can be determined by using Monte Carlo methods and the searching percolation trace method, which is developed in this paper. The study reports on practical examples demonstrating that results from the methods discussed in this paper are in agreement with those from field hydrogeological surveys and measurements.  相似文献   

19.
 Being a non-linear method based on a rigorous formalism and an efficient processing of various information sources, the Bayesian maximum entropy (BME) approach has proven to be a very powerful method in the context of continuous spatial random fields, providing much more satisfactory estimates than those obtained from traditional linear geostatistics (i.e., the various kriging techniques). This paper aims at presenting an extension of the BME formalism in the context of categorical spatial random fields. In the first part of the paper, the indicator kriging and cokriging methods are briefly presented and discussed. A special emphasis is put on their inherent limitations, both from the theoretical and practical point of view. The second part aims at presenting the theoretical developments of the BME approach for the case of categorical variables. The three-stage procedure is explained and the formulations for obtaining prior joint distributions and computing posterior conditional distributions are given for various typical cases. The last part of the paper consists in a simulation study for assessing the performance of BME over the traditional indicator (co)kriging techniques. The results of these simulations highlight the theoretical limitations of the indicator approach (negative probability estimates, probability distributions that do not sum up to one, etc.) as well as the much better performance of the BME approach. Estimates are very close to the theoretical conditional probabilities, that can be computed according to the stated simulation hypotheses.  相似文献   

20.
 This paper deals with the problem of spatial data mapping. A new method based on wavelet interpolation and geostatistical prediction (kriging) is proposed. The method – wavelet analysis residual kriging (WARK) – is developed in order to assess the problems rising for highly variable data in presence of spatial trends. In these cases stationary prediction models have very limited application. Wavelet analysis is used to model large-scale structures and kriging of the remaining residuals focuses on small-scale peculiarities. WARK is able to model spatial pattern which features multiscale structure. In the present work WARK is applied to the rainfall data and the results of validation are compared with the ones obtained from neural network residual kriging (NNRK). NNRK is also a residual-based method, which uses artificial neural network to model large-scale non-linear trends. The comparison of the results demonstrates the high quality performance of WARK in predicting hot spots, reproducing global statistical characteristics of the distribution and spatial correlation structure.  相似文献   

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

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