天山托木尔峰科其喀尔巴西冰川表面运动速度特征分析

天山托木尔峰科其喀尔巴西冰川是典型的树枝状山谷冰川,利用3组(6期)ASTER遥感影像通过COSI-corr软件反演了该冰川表面运动速度.与花杆测量数据进行对比,反演冰川表面运动速度平均绝对误差为3.1 m·a-1,相对误差为11.9%,二者在空间上的分布基本一致,表明其反演精度符合要求.在此基础上,分析冰川表面运动速...     

Multiscale modeling of ice deformation behavior

Understanding the flow of ice in glaciers and polar ice sheets is of increasing relevance in a time of potentially significant climate change. The flow of ice has hitherto received relatively little attention from the structural geological community. This paper aims to provide an overview of methods and results of ice deformation modeling from the single crystal to the polycrystal scale, and beyond to the scale of polar ice sheets. All through these scales, various models have been developed to understand, describe and predict the processes that operate during deformation of ice, with the aim to correctly represent ice rheology and self-induced anisotropy. Most of the modeling tools presented in this paper originate from the material science community, and are currently used and further developed for other materials and environments. We will show that this community has deeply integrated ice as a very useful “model” material to develop and validate approaches in conditions of a highly anisotropic behavior. This review, by no means exhaustive, aims at providing an overview of methods at different scales and levels of complexity.     

The subglacial thermal organisation (STO) of ice sheets

This paper examines ice-sheet wide variations in subglacial thermal regime and ice dynamics using the landform record exposed on the beds of former mid-latitude ice sheets (the Laurentide, Cordilleran, Fennoscandian and British-Irish Ice Sheets). We compare the landform patterns beneath these former ice sheets to the flow organisation beneath parts of the contemporary Antarctic Ice Sheet inferred from RADARSAT-1 Antarctic Mapping Project (RAMP) data. The evidence preserved in the landform record and observed on contemporary ice masses can be grouped into four major ice-dynamical components that collectively define the subglacial thermal organisation (STO) of ice sheets. These ice-dynamical components are frozen-bed patches, ice streams, ice-stream tributaries and lateral shear zones. Frozen-bed patches appear at a wide range of spatial scales, spanning four orders of magnitude. In some areas, frozen-bed zones comprise large proportions of the bed (e.g. near the ice divide in continental areas), whilst in other areas they constitute isolated “islands” in areas dominated by thawed-bed conditions. Ice streams, narrow zones of fast flow in ice sheets that are otherwise dominated by slow sheet flow, are also common features of Quaternary ice sheets. Tributaries to ice streams flow at velocities intermediate between full ice-stream and sheet flow, and may divert ice drainage from one primary ice-stream corridor to an adjacent one. Sharp lateral boundaries between landforms indicate sliding and non-sliding conditions, respectively. These lateral boundaries represent important discontinuities in the glacial landscape and mark the location of shear zones between thawed-bed ice streams and intervening frozen-bed areas. We use the landform evidence in the area around Great Bear Lake, Canada to trace the evolution of an ice-stream web through time, demonstrating that frozen-bed patches are integral components of this complex system. We conclude that frozen-bed patches are important for the stability of ice sheets because they laterally constrain and isolate peripheral drainage basins and their ice streams.     

东濮凹陷三叠系砂岩储层构造裂缝特征及形成期次

根据露头、 岩心、 薄片、 古地磁、 成像测井、 岩石力学实验等资料;对东濮凹陷北部三叠系砂岩储层构造裂缝特征及形成期次进行分析。东濮凹陷三叠系砂岩储层主要发育北北东、 北东和近东西向3组高角度构造裂缝;裂缝走向近于平行或垂直主断层;倾向为北西和南东。二马营组裂缝较为发育;和尚沟组和刘家沟组裂缝相对不发育。储层裂缝为3期;分别为燕山早期的北西向扩张裂缝和北北西、 北西西向剪切裂缝;燕山晚期形成北东向拉张裂缝;喜马拉雅期形成北北东、 北东向拉张裂缝、 北北东向剪切裂缝和东西向扩张裂缝。古、 今构造应力场和岩石力学性质的各向异性分析表明;现今构造应力场最大水平主应力和岩石抗压强度的各向异性均影响燕山早期北西向裂缝的张开度及有效性。储层裂缝表现为两期;分别为燕山晚期的北东向拉张裂缝、 喜马拉雅期的北北东、 北东向拉张裂缝;北北东向剪切裂缝和东西向扩张裂缝。研究结果阐明了东濮凹陷北部三叠系砂岩储层构造裂缝特征及发育规律。     

Characterization of pebble fabrics in modern terrestrial glacigenic sediments

ABSTRACT Pebble fabric data are available from several facies of glacigenic sediments deposited by modern glaciers, where sedimentary processes can be observed or inferred with relatively little ambiguity. Over 100 samples from contemporary environments illustrate fabrics characterizing melt-out till, deformed and undeformed lodgement till, sediment flow deposits and ice slope colluvium. Lodgement till fabric variability is related to the two-layer structure of these sediments; a structureless, friable upper layer with low shear strength and high consolidation coefficient, overlying a very compact material of horizontal platy structure. Fabric strength (assessed by eigenvalue analysis) is weaker and pebble dip is more dispersed in the upper structureless horizon. Stronger fabrics in the lower platy horizon may be primary depositional fabrics which are destroyed by subglacial shearing to give weaker fabrics in the upper horizon. Alternatively, upper horizon fabrics may be characteristic of all recently-deposited lodgement tills, with stronger fabrics developing at depth by dewatering and consolidation. There is a general reduction in fabric strength and an increase in particle dip associated with the transition from melt-out tills, through undeformed and deformed lodgement tills, to sediment flow deposits and ice slope colluvium. There is, however, considerable overlap in the fabric strengths characteristic of sediment flow deposits and deformed lodgement tills. Fabric data from modern glacial sedimentary facies are used to assist in interpreting the mode of deposition of some Quaternary glacial sediments. Relatively strong fabrics characteristic of melt-out tills and undeformed lodgement tills are more likely to be diagnostic of genesis than weaker fabrics associated with deformed sediments.     

Calving processes and the dynamics of calving glaciers

Calving of icebergs is an important component of mass loss from the polar ice sheets and glaciers in many parts of the world. Calving rates can increase dramatically in response to increases in velocity and/or retreat of the glacier margin, with important implications for sea level change. Despite their importance, calving and related dynamic processes are poorly represented in the current generation of ice sheet models. This is largely because understanding the ‘calving problem’ involves several other long-standing problems in glaciology, combined with the difficulties and dangers of field data collection. In this paper, we systematically review different aspects of the calving problem, and outline a new framework for representing calving processes in ice sheet models. We define a hierarchy of calving processes, to distinguish those that exert a fundamental control on the position of the ice margin from more localised processes responsible for individual calving events. The first-order control on calving is the strain rate arising from spatial variations in velocity (particularly sliding speed), which determines the location and depth of surface crevasses. Superimposed on this first-order process are second-order processes that can further erode the ice margin. These include: fracture propagation in response to local stress imbalances in the immediate vicinity of the glacier front; undercutting of the glacier terminus by melting at or below the waterline; and bending at the junction between grounded and buoyant parts of an ice tongue. Calving of projecting, submerged ‘ice feet’ can be regarded as a third-order process, because it is paced by first- or second-order calving above the waterline.First-order calving can be represented in glacier models using a calving criterion based on crevasse depth, which is a function of longitudinal strain rate. Modelling changes in terminus position and calving rates thus reduces to the problem of determining the ice geometry and velocity distribution. Realistic solutions to the problem of modelling ice flow therefore depend critically on an appropriate choice of sliding law. Models that assume that basal velocities are controlled by basal drag can replicate much of the observed behaviour of calving glaciers with grounded termini, but an important limitation is that they cannot be used to model floating glacier termini or ice shelves. Alternative sliding laws that parameterise drag from the glacier margins provide more flexible and robust ways of representing calving in ice sheet models. Such models can explain a remarkable range of observed phenomena within a simple, unifying framework, including: downglacier increases in velocity and strain rates where basal and/or lateral drag diminishes; flow acceleration in response to thinning through time; the tendency for glaciers to stabilise at ‘pinning points’ in relatively shallow water or fjord narrowings; the constraints on ice shelf stability; and the contrasts in calving rates between tidewater and freshwater calving glaciers. Many unresolved issues remain, however, including the role played by the removal of backstress in the acceleration of retreating calving glaciers, and the controls on melting at and below the waterline.     

The Role of Pre-existing Mechanical Anisotropy on Shear Zone Development within Oceanic Mantle Lithosphere: an Example from the Oman Ophiolite

Structural and fabric analysis of the well-exposed Hilti mantlesection, Oman ophiolite, suggests that shear zone development,which may have resulted from oceanic plate fragmentation, wasinfluenced by pre-existing mantle fabric present at the paleo-ridge.Detailed structural mapping in the mantle section revealed agently undulating structure with an east–west flow direction.A NW–SE strike-slip shear zone cuts across this horizontalstructure. The crystal preferred orientation (CPO) of olivinewithin the foliation is dominated by (010) axial patterns ratherthan more commonly observed (010)[100] patterns, suggestingthat the horizontal flow close to the Moho involved non-coaxialflow. Olivine CPO within the shear zone formed at low temperatureis characterized by (001)[100] patterns and a sinistral senseof shear. The olivine CPO becomes weaker with progressive mylonitizationand accompanying grain size reduction, and ultimately developsinto an ultra-mylonite with a random CPO pattern. The olivine[010]-axis is consistently sub-vertical, even where the horizontalfoliation has been rotated to a sub-vertical orientation withinthe shear zone. These observations suggest that the primarymechanical anisotropy (mantle fabric) has been readily transformedinto a secondary structure (shear zone) with minimum modification.This occurred as a result of a change of the olivine slip systemsduring oceanic detachment and related tectonics during cooling.We propose that primary olivine CPO fabrics may play a significantrole in the subsequent structural development of the mantle.Thus, the structural behavior of oceanic mantle lithosphereduring subduction and obduction may be strongly influenced byinitial mechanical anisotropy developed at an oceanic spreadingcenter. KEY WORDS: mantle lithosphere; anisotropy; shear zone; olivine CPO; Oman ophiolite     

Pleistocene and Holocene glacier thicknesses, transport histories and dynamics inferred from SEM microtextures on quartz particles

Recent analyses of microtextures on quartz particles (63–2000 μm) from Quaternary tills in Antarctica, Germany, southern Ontario, western Wyoming, Tibet, the Austrian Alps, and Mount Kenya show that glacial fracture and abrasion microfeatures may be used to infer the thickness, transport history and ice dynamics of Pleistocene and Holocene glaciers. Quartz sands emplaced by continental and mountain ice were studied by SEM after transport over variable distances in glaciers estimated to range from 150 m to 1500 m in thickness. Relative differences in ice thicknesses, distances of transport, and/or ice dynamics appear to have determined the frequency of occurrence and type of microtextures occurring on sand-size particles. Subparallel fracture microfeatures tend to increase in frequency over a greater proportion of particle surfaces with increasing ice thickness and distance of transport. Conchoidal fractures, the most typical in quartz, and to some degree crescentic gouges abound on fragments emplaced by continental ice. Other possible fracture and fragmentation mechanisms, producing features of generally glacial origin, involve low velocity impacts induced by stick-slip mechanisms, under variable cryostatic stresses, producing fracturing and abrasion across particle surfaces. Their generation implies high local contact stresses associated with high strain rates.     

Using volatiles in magma to decipher subglacial eruption dynamics

Eruptions beneath ice sheets and glaciers can generate hazardous ash plumes and powerful meltwater floods, as demonstrated by the recent Icelandic eruptions at Eyjafjallajökull in 2010 and at Grímsvötn in 2011. A key scientific objective for volcanologists is to better understand the factors controlling subglacial eruptions, but the eruptions are mostly hidden beneath ice that is often hundreds of metres thick, thereby preventing direct observation. New approaches are therefore needed to reconstruct the factors driving explosive activity and the response of the overlying ice. The dissolved volatile content, preserved in glassy volcanic rock, offers a useful means of reconstructing palaeo‐ice thicknesses. However, for subglacial rhyolite at least, there seems to be little or no correlation between loading pressure and eruptive style. Instead, there is a strong correlation between the pre‐eruptive volatile content, degassing path and eruptive behaviour. It seems that the style of many subglacial eruptions is controlled by the same mechanisms as subaerial eruptions, with explosivity strongly influenced by degassing and magmatic fragmentation.     

Anisotropy of magnetic susceptibility and rock magnetic applications in the Deccan volcanic province based on some case studies

Anisotropy of Magnetic Susceptibility (AMS) as a tool has been explored here to investigate the nature of petrofabrics in Deccan Volcanic Province (DVP) of west-central Indian region by representative sampling in typical pahoehoe and rubbly pahoehoe lava flows, dykes within flows, shear zone and the impact crater units. The rock magnetic analysis indicate varying degree of concentration of titanomagnetite compositions dominated by multi domain (MD) to pseudo single domain (PSD) grains favoring shape anisotropy of minerals that form primary fabrics. The pahoehoe type lava flows shows planar oblate fabrics without any preferred orientation of principle susceptibility axis (K₁) depicting crystal settling (of magnetic grains) as chief mechanism of fabric development. The rubbly pahoehoe type lava flow exhibit prolate fabrics with well clustered maximum susceptibility axis within horizontal to sub-horizontal planes depicting their response to viscosity shear. The dykes show well clustered K₁ parallel to it’s plane locked during rapid contractional cooling. The sampling at Lonar impact crater was unable to trace any clear fabric due to impact/shock induced deformation and rather preserve the primary fabrics. Further, the shear zone depict random fabrics demanding more detailed and systematic sampling in both the cases. The present investigation infer that the magnetic mineralogy and magnetic fabric variations in the DVP are controlled by the flow mechanism and style of cooling that is characteristic of the given flow unit or dyke and any secondary or superimposed fabric needs to be examined by critical sampling strategy. While more detailed attempts are required to establish the AMS as a tool to record various aspects including the flow dynamics and rate of effusion in the vast terrain of DVP; the present approach is useful to characterize and correlate the lava flow units and dyke occurrences.     

Rotating micro-structures in Antarctic cold basal ice: implications for glacier flow and its interpretation

Structural analyses were conducted in the basal zone of an Antarctic glacier. The studied basal ice sequence was retrieved from a 20-m-long subglacial tunnel dug at the margin of the glacier and is at the temperature of −17°C. For the first time, rotating clast systems embedded within debris-rich ice were thin-sectioned using specially designed cutting techniques. The observed structures reflect the occurrence of pervasive shearing at the base of the glacier, and can be used as shear sense indicators. In addition, some of these structures provide evidence for the presence of thin liquid films at the time of formation despite the marked freezing temperature of the ice. It is showed here that cautious analysis of deformation structures present in debris-bearing ice may bring insights not only into the flow dynamics of the embedding matrix, but also into the behaviour of the interstitial fluid network at the base of cold glaciers and ice sheets.     

Fracture step structure: geometrical characterization and effects on fluid flow and breakthrough curve

Groundwater flow and solute transport through fractured rock is highly responsive to the hydraulic anisotropy and heterogeneity that are specific to every major fracture. A major fracture is modeled as the combination of some primal master fractures and several splay fractures that branch out from primal master fractures: step structures (or jog parts). Step structures are commonly observed along a major fracture on various scales. Master fractures were formed and developed by shear movement while some splay fractures were formed by extension normal to their wall. This difference in fracturing process may lead to a permeability difference between master fractures and splay fractures which seems to be one of the major factors controlling flow and solute transport through the fracture networks due to its hydraulic anisotropic and heterogeneous features. This study is composed of two major components: (1) identification and characterization of a step structure from borehole data; (2) evaluation of effect of some idealized step structures on breakthrough curve by numerical simulations. The fracture data of four 1000-m boreholes were used to make clear fracture patterns in the Tono area of Japan. Some major fractures were identified using stereographic projection technique. On the basis of these results, several idealized models of a major fracture having a step was constructed for the numerical study. The obtained results from numerical simulations clearly imply that geometry of step structure plays an important role in flow and transport through the fracture networks.     

现代冰川过程与全球环境气候演变

文章从宏观和微观两个方面扼要阐述了现代冰川过程与全球变化之间的关系。南极冰盖和格陵兰冰盖冰川物质平衡值目前还没有确切结论，虽然它与全球海平面的升降密切相关。山地冰川末端进退变化和冰川物质平衡与全球升温对应较好。极地冰盖现代降水中的稳定同位素比率，主要阴、阳离子、生物有机酸、微粒、超痕量重金属元素、宇宙尘埃以及火山灰等杂质的含量，为认识地球现代环境气候状况提供了丰富的资料。极地冰盖冰芯的分析结果为重建过去气候环境提供了大信息量，高保真度和高分辨率的资料为预测未来气候环境奠定了坚实基础，具有其它任何载体无法取代的优越性。山地冰川的现代和过去气候环境记录，对研究全球和区域性气候环境状况与变迁意义重大     

Strength anisotropy of granular material consisting of perfectly round particles

The strength anisotropy of granular materials deposited under gravity has mostly been attributed to elongated particles' tendency to align long axes along the bedding plane direction. However, recent experiments on near‐spherical glass beads, for which preferred particle alignment is inapplicable, have exhibited surprisingly strong strength anisotropy. This study tests the hypothesis that certain amount of fabric anisotropy caused by the anisotropic stress during deposition under gravity can be locked in a circular‐particle deposit. Such locked‐in fabric anisotropy can withstand isotropic consolidation and leads to significant strength anisotropy. 2D discrete element method simulations of direct shear tests on circular‐particle deposits are conducted in this study, allowing for the monitoring of both stress and fabric. Simulations on both monodispersed and polydispersed circular‐particle samples generated under downward gravitational acceleration exhibit clear anisotropy in shear strength, thereby proving the hypothesis. When using contact normal‐based and void‐based fabric tensors to quantify fabric anisotropy in the material, we find that the intensity of anisotropy is discernible but low prior to shearing and is dependent on the consolidation process and the dispersity of the sample. The fact that samples with very low anisotropy intensity measurements still exhibit fairly strong strength anisotropy suggests that current typical contact normal‐based and void‐based second‐order fabric tensor formulations may not be very effective in reflecting the anisotropic peak shear strength of granular materials. Copyright © 2017 John Wiley & Sons, Ltd.     

Impacts of Climate Warming on Alpine Glacier Tourism and Adaptive Measures:A Case Study of Baishui Glacier No.1 in Yulong Snow Mountain,Southwestern China

Alpine glaciers usually feature with best hydrothermal condition in mountain climate,and present beautiful glacier scenery,various glacier landforms,rich biodiversity,and easier accessibility,compared with continental glaciers or ice sheets.Nevertheless,Alpine glaciers are more sensitive to climate warming,and climate warming has seriously affected Alpine glaciers and surroundingenvironment.The quality and attractiveness of Alpine glaciers to tourism has been and will continue to be diminished with tourists...     

Shear bands,crenulations and differentiated layering in ice-mica models

Thin sheets of composite ice-mica have been deformed in order to simulate the development of cleavages in quartz-mica rocks. A strong initial mica preferred orientation was variably oriented to the shortening direction. Deformation parallel to the foliation results in a crenulation type cleavage developing from shear bands initiated after a component of pure shear. Deformation oblique to the foliation produces a differentiated cleavage and involves a large component of shear strain subparallel to the original anisotropy. The strain is accommodated by intra- and intercrystalline processes that produce extensive grain elongation and rearrangement of the ductile matrix, thereby forming ice vs mica rich regions. On the other hand, there is no drastic morphological change when a sample is shortened perpendicular to an original foliation: that is, where the micas lie in the plane of no shear strain. Instead, the mica fabric is strengthened and the grains in the ductile matrix are flattened.Two models are presented for the initiation, propagation and evolution of the observed crenulation versus differentiated cleavage types. These depend on mica stacking and orientation relative to the transverse properties of the sample and also on the direction of anisotropy to the XY plane of the bulk strain ellipsoid. The models invoke shear on planes of high shear strain and rotation of the shear bands and rigid mica grains into a direction approximately parallel to the bulk extension direction.     

A shear‐dilation‐based model for evaluation of hydraulically stimulated naturally fractured reservoirs

The role of shear dilation as a mechanism of enhancing fluid flow permeability in naturally fractured reservoirs was mainly recognized in the context of hot dry rock (HDR) geothermal reservoir stimulation. Simplified models based on shear slippage only were developed and their applications to evaluate HDR geothermal reservoir stimulation were reported. Research attention is recently focused to adjust this stimulation mechanism for naturally fractured oil and gas reservoirs which reserve vast resources worldwide. This paper develops the overall framework and basic formulations of this stimulation model for oil and gas reservoirs. Major computational modules include: natural fracture simulation, response analysis of stimulated fractures, average permeability estimation for the stimulated reservoir and prediction of an average flow direction. Natural fractures are simulated stochastically by implementing ‘fractal dimension’ concept. Natural fracture propagation and shear displacements are formulated by following computationally efficient approximate approaches interrelating in situ stresses, natural fracture parameters and stimulation pressure developed by fluid injection inside fractures. The average permeability of the stimulated reservoir is formulated as a function of discretized gridblock permeabilities by applying cubic law of fluid flow. The average reservoir elongation, or the flow direction, is expressed as a function of reservoir aspect ratio induced by directional permeability contributions. The natural fracture simulation module is verified by comparing its results with observed microseismic clouds in actual naturally fractured reservoirs. Permeability enhancement and reservoir growth are characterized with respect to stimulation pressure, in situ stresses and natural fracture density applying the model to two example reservoirs. Copyright © 2002 John Wiley & Sons, Ltd.     

含裂隙煤层的地震记录模拟

采用线性滑移(Linear Slip Deformation)裂隙介质模型、各向异性弹性波模拟方法,对含裂隙煤层进行了2D2C(二维二分量)单炮记录模拟与分析。结果表明:煤层含裂隙后与不含裂隙时相比,反射时差的变化已不足以识别,但由于波速的改变会使煤层顶底界面的波阻抗发生变化,从而引起波的动力学特征的改变,AVO现象尤其明显。当裂隙为非垂直时,炮点两侧的反射波能量和振幅出现不对称现象。x分量和z分量的特征明显不同,转换波主要记录在x分量上,横波在x分量上有不可忽视的能量。因此,应用各向异性理论、利用多分量转换波研究煤层裂隙的特征是可能的。     

Folding by cataclastic flow: evolution of controlling factors during deformation

Folding at upper crustal levels occurs by bending of beds and flexural slip between beds. As a fold's interlimb angle decreases, changes in bed thickness and limb rotation are accommodated by various mechanisms, depending on deformation conditions. In the elastico-frictional (EF) regime, cataclastic flow may be the dominant mechanism for fold tightening. The Canyon Range (CR) syncline, located in the Sevier belt of central Utah, shows this type of deformation. The fold involves three thick quartzite units, with slight lithological variations between them. Fold tightening took place in the EF regime (<2 km overburden) by cataclastic flow, involving collective movement on a distributed network of fractures and deformation zones (DZs) from the micro- to the outcrop-scale. In detail, the degree of cataclastic deformation varies significantly across the fold due to minor variations in initial bedding thickness, grain size, matrix composition, etc. A cooperative relationship exists across different scales, and the fracture networks result in a fracture shape fabric that is relatively homogeneous at the outcrop-scale.The initial outcrop scale fracture/DZ network geometry is a product of the growth and linking of micro-scale cataclasite zones, which in turn is controlled by primary lithological variations. Once a fracture network forms, the material behavior of the fractured rock is unlike that of the original rock, with sliding of fracture-bound blocks accomplishing ‘block-controlled’ cataclastic flow. Thus, initial lithological variations at the micro-scale largely control the final deformation behavior at the largest scale. During progressive fold tightening, additional factors regulate cataclastic flow, such as fracture/DZ reactivation or healing, during folding. Although initial lithological variations in different units may produce unique network geometries, each unit's behavior may also depend upon the behavior of adjacent units. In the CR syncline, during the initial stages of cataclastic flow, the inherent nature of each quartzite unit results in unit-specific fracture network geometries. As deformation progresses, unit-specific networks begin to interact with those in surrounding units, resulting in feedback mechanisms regulating the later stages of network development. Thus, the nature of cataclastic flow changes dramatically from the initial to the final stages of folding.     

低温裂隙岩体水-热耦合模型研究及数值分析

裂隙渗流会引起裂隙周围岩体中的温度场变化,在低温岩体中其影响更为明显;此外,裂隙水与周围低温岩石介质发生热交换会引起裂隙中的水冰相变过程发生,而裂隙水冻结将阻碍裂隙渗流,引起裂隙渗流场的变化。因此,低温下的裂隙岩体水-热相互作用是一个强耦合过程。考虑裂隙中的水冰相变过程和渗流作用,建立了低温冻结条件下裂隙岩体水-热耦合模型;以冻结法施工为例,考察了低温冻结过程中裂隙水渗流对裂隙冻结交圈的影响。研究结果表明：由于裂隙渗流的存在,距裂隙较远处岩石先冻结,裂隙冻结所需时间远大于周围岩石;裂隙宽度和裂隙水压力差都会影响冻结交圈时间,裂隙越宽、水压力差越大,裂隙冻结需要时间越长;随着冻结时间的推进,裂隙水渗流速度逐渐降低,当裂隙冻结后裂隙渗流停止。最后通过构建随机裂隙网络模型,利用所建立的水-热耦合模型考察了裂隙网络渗流对冻结交圈的影响,说明了在冻结法施工中考虑裂隙的重要性。