Multi-scale permeability estimation for a tropical catchment     

Nick A. Chappell  Stewart W. Franks  Jonny Larenus 《水文研究》1998,12(9):1507-1523

Physically based and spatially distributed modelling of catchment hydrology involves the estimation of block or whole-hillslope permeabilities. Invariably these estimates are derived by calibration against rainfall–runoff response. Rarely are these estimates rigorously compared with parameter measurements made at the small scale. This study uses a parametrically simple model, TOPMODEL, and an uncertainty framework to derive permeability at the catchment scale. The utility of expert knowledge of the internal catchment dynamics (i.e. extent of saturated area) in constraining parameter uncertainty is demonstrated. Model-derived estimates are then compared with core-based measurements of permeability appropriately up-scaled. The observed differences between the permeability estimates derived by the two methods might be attributed to the role of intermediate scale features (natural soil pipes). An alternative method of determining block permeabilities at the intermediate or hillslope scale is described. This method uses pulse-wave tests and explicitly incorporates the resultant effects of phenomena such as soil piping and kinematic wave migration. The study aims to highlight issues associated with parameterizing or validating distributed models, rather than to provide a definitive solution. The fact that the permeability distribution within the Borneo study catchment is comparatively simple, assists the comparisons. The field data were collected in terrain covered by equatorial rainforest. Combined field measurement and modelling programmes are rare within such environments. © 1998 John Wiley & Sons, Ltd.  相似文献   

基于水体光学分类的鄱阳湖悬沙浓度反演方法研究     

况润元  赵艳福  罗卫  张刚华  陈彦兵 《水文》2017,37(6):23-28

鄱阳湖是我国重要的湿地生态系统,对调节流域的水沙变化有着重要作用。由于鄱阳湖湖区面积广,内部差异大,单一的悬浮泥沙反演模型不足以准确反演出湖区的悬浮泥沙浓度。以实测的反射光谱数据、泥沙浓度数据为基础,提出一种基于分类后的反演模型,即根据实测数据的光谱形态特征分类出5种典型的水体类型。在此基础上,将分类后的各类水体分别建立各自合适的反演模型进行反演。结果表明基于水体分类的经验模型反演达到满意效果,平均绝对误差为0.00217g/L,平均相对误差为3.022%。基于分类后的经验反演模型适用于鄱阳湖悬沙浓度分布的监测研究,有助于更加宏观、准确的掌握鄱阳湖泥沙浓度的空间分布和变化,为保持鄱阳湖资源的可持续开发与利用提供决策依据。  相似文献   

河南尉氏县西部地质地球物理综合解译及地热资源远景区预测     

刘志龙  朱怀亮  胥博文  叶高峰  胡志明  李文龙  王兴元 《地质论评》2020,66(5):1446-1455

尉氏县地热资源开发尚处于起步阶段,为适应大规模地热资源开发需求,当前急需探明该区地热资源赋存情况。研究区位于太康隆起次级构造单元中部鼻状构造带及尉氏—薛店次凹带结合部位,本文通过尉氏县城西部所采集的220个大地电磁测深点的二维及三维反演工作,结合二维反射地震剖面,分析了研究区地层、断裂展布构造,并对地热资源远景区进行了预测。维性分析结果表明研究区主构造走向为北西向,二维及三维大地电磁测深反演结果表明研究区4 km以浅地层由上到下可分为四层：新近系—第四系、三叠系、石炭系—二叠系、寒武系—奥陶系;综合二维反射地震及大地电磁测深资料,剖面位置处可识别出三条正断层。选择断裂带周边或穿过断裂带钻井,对水量有一定保证。  相似文献   

CASM: a unified state parameter model for clay and sand     

H. S. Yu 《国际地质力学数值与分析法杂志》1998,22(8):621-653

The purpose of this paper is to present a simple, unified critical state constitutive model for both clay and sand. The model, called CASM (Clay And Sand Model), is formulated in terms of the state parameter that is defined as the vertical distance between current state (v, p′) and the critical state line in v–ln p′ space. The paper first shows that the standard Cam-clay models (i.e. the original and modified Cam-clay models) can be reformulated in terms of the state parameter. Although the standard Cam-clay models prove to be successful in modelling normally consolidated clays, it is well known that they cannot predict many important features of the behavior of sands and overconsolidated clays. By adopting a general stress ratio-state parameter relation to describe the state boundary surface of soils, it is shown that a simple, unified constitutive model (CASM) can be developed for both clay and sand. It is also demonstrated that the standard Cam-clay yield surfaces can be either recovered or approximated as special cases of the yield locus assumed in CASM. The main feature of the proposed model is that a single set of yield and plastic potential functions has been used to model the behaviour of clay and sand under both drained and undrained loading conditions. In addition, it is shown that the behaviour of overconsolidated clays can also be satisfactorily modelled. Simplicity is a major advantage of the present state parameter model, as only two new material constants need to be introduced when compared with the standard Cam-clay models. © 1998 John Wiley & Sons, Ltd.  相似文献   

悬臂式支护结构设计中土层参数分算与合算的探讨   总被引：2，自引：0，他引：2  

杨佳平  郑毅  徐家斌 《世界地质》2005,24(3):293-296

针对基坑支护设计中土层参数的选取问题,以悬臂桩支护结构为例,在基坑支护设计中土层参数分别采用分算与合算两种方法计算,结果表明：当各土层差异较大时,采用两种方法计算得出的最大弯矩和埋深相差较大;采用合算的结果较小,分算的结果较大.采用m法对分算与合算的设计计算结果进行悬臂桩的位移和转角的计算,其结果都能满足工程技术相关规范要求.但合算的计算结果更经济,更符合工程实际需要.  相似文献   

沿107.6°E南北向剖面鄂尔多斯地块及周缘地区地壳结构          下载免费PDF全文

陈洁  陈永顺  郭震  杨挺 《地球物理学报》2020,63(7):2592-2604

鄂尔多斯地块紧邻青藏高原东北缘,位于华北克拉通的西部,在我国中生代、新生代以来东部地区的构造活动中起到了重要作用.对鄂尔多斯及其周缘地区的研究可以提供有关华北克拉通的形成、演化和破坏过程的重要信息.本文选取了纵贯鄂尔多斯的107.6°E附近南北剖面上的44个流动地震台站进行分析,采用接收函数方法,进行Kirchhoff偏移成像,并且结合在该区域内前人的地震面波频散进行联合反演,获得剖面下方的地壳内部精细结构.研究结果显示:(1)莫霍面在鄂尔多斯北部较平缓,约45km深;在鄂尔多斯南部有所加深,达到50km;其北边的河套盆地的地壳厚度约为50km;南边的渭河盆地到秦岭地区及四川盆地的地壳厚度从约为40km增厚到47～50km.(2)河套盆地下方存在大规模的低速异常,最深可达25km,反映了其显著的拉张构造和沉积历史.(3)秦岭造山带下方的低速异常对应于其主要为长英质的地壳组分,可能是由于中生代的拆沉作用导致的地壳下部基性岩石层的缺失.(4)以38°N为界的鄂尔多斯地块,南北部地壳速度结构存在差异,可能表明了这两部分经历的构造历史不同.  相似文献   

空-地-井重力异常正则化协同密度反演方法          下载免费PDF全文

王泰涵  马国庆  熊盛青  刘燕戌  曾昭发  李丽丽  孟庆发 《地球物理学报》2020,63(7):2737-2750

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.  相似文献   

Incorporating known petrophysical model in the seismic full-waveform inversion using the Gramian constraint     

M. Malovichko  N. Khokhlov  N. Yavich  M. S. Zhdanov 《Geophysical Prospecting》2020,68(4):1361-1378

In this paper, we develop a general approach to integrating petrophysical models in three-dimensional seismic full-waveform inversion based on the Gramian constraints. In the framework of this approach, we present an example of the frequency-domain P-wave velocity inversion guided by an electrical conductivity model. In order to introduce a coupling between the two models, we minimize the corresponding Gramian functional, which is included in the Tikhonov parametric functional. We demonstrate that in the case of a single-physics inversion guided by a model of different physical type, the general expressions of the Gramian functional and its gradients become simple and easy to program. We also prove that the Gramian functional has a non-negative quadratic form, so it can be easily incorporated in a standard gradient-based minimization scheme. The developed new approach of seismic inversion guided by the known petrophysical model has been validated by three-dimensional inversion of synthetic seismic data generated for a realistic three-dimensional model of the subsurface.  相似文献   

Multiscale phase inversion for 3D ocean-bottom cable data     

Lei Fu  Zongcai Feng  Gerard T. Schuster 《Geophysical Prospecting》2020,68(3):786-801

We invert three-dimensional seismic data by a multiscale phase inversion scheme, a modified version of full waveform inversion, which applies higher order integrations to the input signal to produce low-boost signals. These low-boost signals are used as the input data for the early iterations, and lower order integrations are computed at the later iterations. The advantages of multiscale phase inversion are that it (1) is less dependent on the initial model compared to full waveform inversion, (2) is less sensitive to incorrectly modelled magnitudes and (3) employs a simple and natural frequency shaping filtering. For a layered model with a three-dimensional velocity anomaly, results with synthetic data show that multiscale phase inversion can sometimes provide a noticeably more accurate velocity profile than full waveform inversion. Results with the Society of Exploration Geophysicists/European Association of Geoscientists and Engineers overthrust model shows that multiscale phase inversion more clearly resolves meandering channels in the depth slices. However, the data and model misfit functions achieve about the same values after 50 iterations. The results with three-dimensional ocean-bottom cable data show that, compared to the full waveform inversion tomogram, the three-dimensional multiscale phase inversion tomogram provides a better match to the well log, and better flattens angle-domain common image gathers. The problem is that the tomograms at the well log provide an incomplete low-wavenumber estimate of the log's velocity profile. Therefore, a good low-wavenumber estimate of the velocity model is still needed for an accurate multiscale phase inversion tomogram.  相似文献   

The effectiveness of a pseudo-inverse extended Born operator to handle lateral heterogeneity for imaging and velocity analysis applications     

Abdullah Alali  Bingbing Sun  Tariq Alkhalifah 《Geophysical Prospecting》2020,68(4):1154-1166

Wave equation–based migration velocity analysis techniques aim to construct a kinematically accurate velocity model for imaging or as an initial model for full waveform inversion applications. The most popular wave equation–based migration velocity analysis method is differential semblance optimization, where the velocity model is iteratively updated by minimizing the unfocused energy in an extended image volume. However, differential semblance optimization suffers from artefacts, courtesy of the adjoint operator used in imaging, leading to poor convergence. Recent findings show that true amplitude imaging plays a significant role in enhancing the differential semblance optimization's gradient and reducing the artefacts. Here, we focus on a pseudo-inverse operator to the horizontally extended Born as a true amplitude imaging operator. For laterally inhomogeneous models, the operator required a derivative with respect to a vertical shift. Extending the image vertically to evaluate such a derivative is costly and impractical. The inverse operator can be simplified in laterally homogeneous models. We derive an extension of the approach to apply the full inverse formula and evaluate the derivative efficiently. We simplified the implementation by applying the derivative to the imaging condition and utilize the relationship between the source and receiver wavefields and the vertical shift. Specifically, we verify the effectiveness of the approach using the Marmousi model and show that the term required for the lateral inhomogeneity treatment has a relatively small impact on the results for many cases. We then apply the operator in differential semblance optimization and invert for an accurate macro-velocity model, which can serve as an initial velocity model for full waveform inversion.  相似文献