Hengill,a high temperature thermal area in Iceland

The results of geological, geochemical and isotopic investigations in the Hengill thermal area in Iceland are presented. This area may be regarded as typical of the Icelandic high temperature areas. The Hengill area is mainly built up of pillow lavas and hyaloclastites, which were piled up in subglacial eruptions. Its western half is traversed by a very active NE-SW trending fault zone about 5 km broad, within which are several eruptive fissures of postglacial age. The volcanic rocks are basalts of various kinds, but minor occurrences of intermediate and rhyolitic rocks are also found. Solfataras and fumaroles occur mainly at higher altitude. Within the fault zone they have a linear distribution along faults. Hot water springs are found in the lowlands on the southern border of the thermal area. The concentrations of chlorine, tritium and deuterium in the water indicate that the hot springs on the lowland are fed by a deep lying aquifer which receives water from a recharge area situated in the western branch of the active volcanic zone upto 70 km to the northeast. The hot springs at higher altitude are the outlets of a secondary circulation system fed by local rain water. This water derives its heat partly from the heat content of the rocks penetrated and partly by steam and gas boiled off from the deeper lying aquifer. Preliminary results from other high temperature areas in Iceland indicate the presence of similar hydrological systems.     

北京南口—孙河断裂带水热系统特征与成因分析

中低温对流型地热资源在华北地区广泛分布,是一种清洁的替代能源.与活动断裂带相关的水热型地热资源是中低温地热系统的重要组成部分.本文基于高精度重力测量、微动测深及钻孔温度测量等数据,从热源、通道、储层和盖层四个方面探讨了南口—孙河断裂带水热系统特征.低重力异常揭示的燕山期花岗二长岩、闪长岩岩体范围为23.8 km~2和14.3 km~2,放射性测井数据计算得到其生热率均值为3.14μW·m~(-3),侏罗系火山岩生热率均值为1.65μW·m~(-3),隐伏岩体和火山岩均难以构成地热系统的附加热源.重力异常显示南口—孙河断裂带宽度约500~800 m,断裂带切割蓟县系雾迷山组白云岩热储层.钻井温度曲线显示断裂带内水热活动强烈,说明该断裂带是导水、导热的重要通道.断裂带南西侧马池口一带第四系松散层与侏罗系火山岩形成了热储盖层,微动测深显示火山岩最大厚度约1500 m.综上源、通、储、盖四个要素分析,该地热系统为热传导一对流复合型,来自京西北山区的大气降水经远距离径流深循环吸收地层热量后沿南口—孙河断裂上移到达裂隙发育的白云岩地层中形成热水.总之,沿南口—孙河断裂带具备了良好的地热地质条件,可达到规模开采的条件.     

Thermal structure of the yakedake volcano, Japan: Karukaya and Takara geothermal areas

The hydrothermal water balance and the thermal structure of Yakedake volcano and its vicinity are considered quantitatively. The hydrothermal activity is intense in the valleys at the western foot of the volcano and the Nakanoyu area. The total hot water flow from the discharge area amounts to 2.07 × 10⁴1/min, about 60% of which discharges from the Shinhodaka area alone. There are some large basins (Abodaira and others) in which the rocks are mainly tuff breccia and volcanic products showing very high permeability for water. The total area of the water recharge zone amounts to 18.2 × 10⁶m³. A model for the hydrothermal system within Yakedake volcano is proposed and from the results of boreholes, the thermal and geological structures of the Karukaya and Takara geothermal areas are also presented.Attempts were also made to estimate the subsurface temperature distribution from the observed near-surface ground temperatures. Results of three-dimensional conduction model calculations indicate that the subsurface temperatures are high in the central part of the crater and in the areas with self-flowing springs along the rivers. The obtained isotherms encircle the volcanic center of Yakedake.     

Geology of the platanares geothermal area, Departamento de Copán, Honduras

The Platanares geothermal area, Departamento de Copán, Honduras, is located within a graben that is complexly faulted. The graben is bounded on the north by a highland composed of Paleozoic (?) metamorphic rocks in contact with Cretaceous - Tertiary redbeds of unknown thickness. These are unconformably overlain by Tertiary andesitic lavas, rhyolitic ignimbrites, and associated sedimentary rocks. The volcanic rocks are mostly older than 14 Ma, and thus are too old to represent the surface expression of an active crustal magma body. Thermal fluids that discharge in the area are heated during deep circulation of meteoric water along faults in a region of somewhat elevated heat flow. Geothermometry based upon the chemical composition of thermal fluids from hot springs and from geothermal gradient coreholes suggests that the reservoir equilibrated at temperatures as high as 225 to 240°C, within the Cretaceous redbed sequence. Three continuously cored geothermal gradient holes have been drilled; fluids of about 165°C have been produced from two drilled along a NW-trending fault zone, from depths of 250 to 680 m. A conductive thermal gradient of 139°C/km, at a depth of 400 m, was determined from the third well, drilled 0.6 km west of that fault zone. These data indicate that the Platanares geothermal area holds considerable promise for electrical generation by moderate- to hightemperature geothermal fluids.     

青藏高原东南缘及邻区地壳介质散射强度变化特征

基于远震扰动场方法,利用47个数字地震观测台记录的10次远震记录对青藏高原东南缘及邻区地壳介质散射强度进行了研究.结果表明研究区内地壳介质横向非均匀性强烈,扬子块体地壳介质散射强度与青藏块体存在较为明显的差异;地壳介质散射强度梯度带与主要断裂带分布趋势一致,龙门山断裂带、安宁河断裂带与地壳介质散射强度梯度带具有较强的相关性,散射强度高值区偏向青藏块体一侧;强震震中沿散射强度梯度带分布,其震中偏向散射强度高值区一侧;散射强度呈现纵向非均匀性特征,下地壳介质散射强度高于上地壳.研究区强、弱散射强度的空间分布可能与区内强烈的地质构造运动、频繁的地震及火山活动所引起深部物质破碎、熔融及其沿断裂带的上涌与运移有关.     

Propagation of seasonal temperature signals into an aquifer upon bank infiltration

Infiltrating river water carries the temperature signal of the river into the adjacent aquifer. While the diurnal temperature fluctuations are strongly dampened, the seasonal fluctuations are much less attenuated and can be followed into the aquifer over longer distances. In one-dimensional model with uniform properties, this signal is propagated with a retarded velocity, and its amplitude decreases exponentially with distance. Therefore, time shifts in seasonal temperature signals between rivers and groundwater observation points may be used to estimate infiltration rates and near-river groundwater velocities. As demonstrated in this study, however, the interpretation is nonunique under realistic conditions. We analyze a synthetic test case of a two-dimensional cross section perpendicular to a losing stream, accounting for multi-dimensional flow due to a partially penetrating channel, convective-conductive heat transport within the aquifer, and heat exchange with the underlying aquitard and the land surface. We compare different conceptual simplifications of the domain in order to elaborate on the importance of different system elements. We find that temperature propagation within the shallow aquifer can be highly influenced by conduction through the unsaturated zone and into the underlying aquitard. In contrast, regional groundwater recharge has no major effect on the simulated results. In our setup, multi-dimensionality of the flow field is important only close to the river. We conclude that over-simplistic analytical models can introduce substantial errors if vertical heat exchange at the aquifer boundaries is not accounted for. This has to be considered when using seasonal temperature fluctuations as a natural tracer for bank infiltration.     

Three-dimensional temperature structure of a ridge-transform-ridge system

We have solved the problem of three-dimensional heat transfer by conduction and material transport within a ridge-transform-ridge system. The boundary conditions are those of a finite plate with temperatures fixed at the top and bottom. We do not solve for the pattern of flow, so must arbitrarily describe the material transport. Within triangular regions beneath the spreading centers, flow is assumed to be vertical and at constant velocity, with the volume of upwelling matching the volume of horizontal transport by the plates. The steady state solution is approached by repeated finite steps in time, each involving two stages. First, to represent material transport, temperatures in the grid of points are shifted horizontally and/or vertically to the point at which the associated particle would arrive given the time period and assumed velocity field. Second, a three-dimensional Fourier Transform technique is used to calculate the cooling which would occur in the static case during the same period.

We find that the geometry of flow in the upwelling region, here parameterized by velocity of upwelling or width of the upwelling zone, is more important than spreading rate in determining the temperature structure. The results indicate that conductive cooling cannot alone account for deepening of the median valley towards the fracture zone. For a 10 m.y. offset on the transform fault, conductive cooling across the fracture zone to the older lithosphere affects the upwelling material only to a distance of about 10 km from the fault, in contrast to the deepening of the median valley over several tens of kilometers observed on the Mid-Atlantic Ridge.     

Geological structure of hydrothermal fields in the taupo volcanic zone,New Zealand

Wairakei hydrothermal field is underlain by an acid volcanic sequence consisting of the following units: Recent pumice cover, Wairakei Breccia, Huka Falls Formation, Haparangi Rhyolite, Waiora Formation, Waiora Valley Andesite, Wairakei Ignimbrites, Ohakuri Group. The stratified volcanic sequence is draped over a basement horst and thickens eastward and westward into adjoining volcanotectonic depressions. These major depressions have grown progressively during the Quaternary by differential subsidence along active faults and were not the direct result of collapse following major ignimbrite and rhyolite eruptions. The bulk of the steam production is obtained from a thick aquifer of pumice breccias (Waiora aquifer), which is capped by lacustrine shales of the Huka Falls Formation. The best production is obtained from fault zones intersecting the lower part of this aquifer. The Ohakuri Group, underlying the Wairakei Ignimbrites, constitutes a lower, pumice breccia aquifer, not yet directly exploited by drill-holes. Hydrothermal water at 265°C is fed from the lower to the upper aquifer through linear fissures in the Wairakei Ignimbrites, principally at the cress of a small structural dome. These fissures are related to active, north-east striking, normal faults, having a small dextral transcurrent component. Major zones of heat liberation have been localised at the intersection of secondary, north-west cross faults. Fossil, hydrothermal, mud-flow, conglomerates intercalated in the mid-Pleistocene Huka Falls Formation, suggest that hydrothermal activity at Wairakei is at least 500,000 years old. Isopachs on the early Pleistocene Waiora Valley Andesite indicate that the andesite plug and associated flows were crupted at the same fault intersections, now controlling heat flow into the Waiora aquifer. The inferred great age of the hydrothermal system poses problems in maintaining heat flow, and a model of a large, semi-permanent granitic batholith leaking super-critical hydrothermal fluids up active fault zones to heat near-surface aquifers of meteoric water is proposed. Temperature and pressure equilibria of such a magma body may be maintained by gaseous diffusion. Structural and heat flow evidence suggest that the most probable location of the magma body is to the south of the hydrothermal field below a large, late Pleistocene, rhyolitic eruptive centre. Other hydrothermal fields in the Taupo Volcanic Zone — Orakeikorako, Waiotapu and Kawerau — have similar geological sequences and structures to Wairakei. Important features in common include:

1. 1)
   Proximity to a late Pleistocene, rhyolitic eruptive centre typified by ignimbrites, non-welded ash flows and rhyolite domes and sills, which is considered to directly overlie a large, semipermanent magma chamber.     
   
   

Ground water recharge and flow characterization using multiple isotopes

Stable isotopes of delta(18)O, delta(2)H, and (13)C, radiogenic isotopes of (14)C and (3)H, and ground water chemical compositions were used to distinguish ground water, recharge areas, and possible recharge processes in an arid zone, fault-bounded alluvial aquifer. Recharge mainly occurs through exposed stream channel beds as opposed to subsurface inflow along mountain fronts. This recharge distribution pattern may also occur in other fault-bounded aquifers, with important implications for conceptualization of ground water flow systems, development of ground water models, and ground water resource management. Ground water along the mountain front near the basin margins contains low delta(18)O, (14)C (percent modern carbon [pmC]), and (3)H (tritium units [TU]), suggesting older recharge. In addition, water levels lie at greater depths, and basin-bounding faults that locally act as a flow barrier may further reduce subsurface inflow into the aquifer along the mountain front. Chemical differences in ground water composition, attributed to varying aquifer mineralogy and recharge processes, further discriminate the basin-margin and the basin-center water. Direct recharge through the indurated sandstones and mudstones in the basin center is minimal. Modern recharge in the aquifer is mainly through the broad, exposed stream channel beds containing coarse sand and gravel where ground water contains higher delta(18)O, (14)C (pmC), and (3)H (TU). Spatial differences in delta(18)O, (14)C (pmC), and (3)H (TU) and occurrences of extensive mudstones in the basin center suggest sluggish ground water movement, including local compartmentalization of the flow system.     

The permeability of the Elkhorn fault zone,South Park,Colorado

The purposes of this study are to use both field and modeling approaches to characterize the permeability of a fault and to assess the role of the fault on regional ground water flow. The study subject is the Elkhorn fault, a low-angle reverse fault that brings Precambrian crystalline rocks over the sediments of Colorado's South Park Basin. The fault is hypothesized to act as a low-permeability barrier to flow, restricting interaction between the crystalline aquifer and the basin sediments. To test this hypothesis and to better predict the permeability structure of the fault, we synthesized geologic data to create a geologic model of the fault, conducted aquifer tests to estimate the hydrogeologic properties of the fault zone, and used ground water modeling to test the influence of a range of hydraulic properties for the fault zone on ground water flow in the region. Our study suggests that the fault is a low-permeability feature. Estimated heads are best matched to observations by modeling the fault as a 10-foot-thick interval of low-permeability fault gouge. Steady-state flow models show that much of the flow in the study area is topographically driven near land surface. Flow rates decrease with depth in the aquifers. In the footwall, ground water moves updip in the Michigan-San Isabel syncline to discharge in the South Park Basin. In the hanging wall, ground water moves east to a regional ground water divide. Sensitivity analyses indicate that hydraulic heads are most sensitive to changes in hydraulic conductivity and recharge.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei,Shaanxi Province,China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.

     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei,Shaanxi Province,China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei, Shaanxi Province, China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.     

Influence of ancient thrust faults on the hydrogeology of the Blue Ridge Province

The Blue Ridge Province contains ubiquitous northeast-southwest-trending thrust faults or smaller thrust "slivers" that greatly impact the nature and character of ground water flow in this region. Detailed investigations at a field site in Floyd County, Virginia, indicate that high-permeability zones occur in the brittle crystalline rocks above these thrust faults. Surface and borehole geophysics, aquifer tests, and chlorofluorocarbon and geochemical data reveal that the shallow saprolite aquifer is separated from the deeper fault-zone aquifer by a low-fracture permeability bedrock confining unit, the hydraulic conductivity of which has been estimated to be six orders of magnitude less than the conductivity of the fault zones at the test site. Within the Blue Ridge Province, these fault zones can occur at depths of 300 m or more, can contain a significant amount of storage, and yield significant quantities of water to wells. Furthermore, it is expected that these faults may compartmentalize the deep aquifer system. Recharge to and discharge from the deep aquifer occurs by slow leakage through the confining unit or through localized breach zones that occur where quartz accumulated in high concentrations during metamorphism and later became extensively fractured during episodes of deformation. The results of this investigation stress the importance of thrust faults in this region and suggest that hydrogeologic models for the Blue Ridge Province include these ancient structural features. Faults in crystalline-rock environments may not only influence the hydrology, they may dominate the flow characteristics of a region.     

Structural control on uranium mineralization in South China: Implications for fluid flow in continental strike-slip faults

South China is the most important uranium producer in the country. Much of the Mesozoic-Cenozoic geology of this area was dominated by NNE-trending intracontinental strike-slip faulting that resulted from oblique subduction of the paleo-Pacific plate underneath the eastern China continent. This strike-slip fault system was characterized by transpression in the early-mid Jurassic and by transtension from the latest Jurassic through Cretaceous to early Tertiary. Most uranium ore deposits in South China are strictly fault-hosted and associated with mid-late Mesozoic granitic intrusions and volcanic rocks, which formed under transpression and transtension regimes, respectively. Various data demonstrate that the NNE-trending strike-slip faults have played critical roles in the formation and distribution of hydrothermal uranium deposits. Extensive geochronological studies show that a majority of uranium deposits in South China formed during the time period of 140–40 Ma with peak ages between 87–48 Ma, coinciding well with the time interval of transtension. However, hydrothermal uranium deposits are not uniformly distributed along individual strike-slip fault. The most important ore-hosting segments are pull-apart stepovers, splay structures, extensional strike-slip duplexes, releasing bends and fault intersections. This non-uniform distribution of ore occurrences in individual fault zone reflects localization of hydrothermal fluids within those segments that were highly dilational and thus extremely permeable. The unique geometric patterns and structural styles of strike-slip faults may have facilitated mixing of deeply derived and near-surface fluids, as evidenced by stable isotopic data from many uranium deposits in South China. The identification of fault segments favorable for uranium mineralization in South China is important for understanding the genesis of hydrothermal ore deposits within continental strike-slip faults, and therefore has great implications for exploration strategies.     

Including Vertical Fault Structures in Layered Groundwater Flow Models

Accurate representation of groundwater flow and solute transport requires a sound representation of the underlying geometry of aquifers. Faults can have a significant influence on the structure and connectivity of aquifers, which may allow permeable units to connect, and aquifers to seal when juxtaposed against lower permeability units. Robust representation of groundwater flow around faults remains challenging despite the significance of faults for flow and transport. We present a methodology for the inclusion of faults utilizing the unstructured grid features of MODFLOW-USG and MODFLOW 6. The method focuses on the representation of fault geometries using non-neighbor connections between juxtaposed layers. We present an illustration of the method for a synthetic fluvial aquifer. The combined impact of the heterogeneous aquifer and fault offset is clearly visible where channel features at different depths in the aquifer were connected at the fault. These results highlight the importance of representing fault features in groundwater flow models.     

Modeling a Large-Scale Historic Aquifer Test: Insight into the Hydrogeology of a Regional Fault Zone

Faults can act as flow barriers or conduits to groundwater flow by introducing heterogeneity in permeability. We examine the hydrogeology of the Sandwich Fault Zone, a 137 km long zone of high-angle faults in northern Illinois, using a large-scale historic aquifer test. The fault zone is poorly understood at depth due to the majority of the faults being buried by glacial deposits and its near-vertical orientation which limits geologic sampling across faults. The aquifer test—perhaps one of the largest in terms of overall withdrawal in North American history—was conducted in 1942 at a facility adjacent to the fault zone. More than 34,000 m³/day was pumped for 37 days from nine multiaquifer wells open to the stratified Cambrian-Ordovician sandstone aquifer system. We modeled the aquifer test using a transient MODFLOW-USG model and simulated pumping wells with the CLN package. We tested numerous fault core/damage zone conceptualizations and calibrated to drawdown values recorded at production and observation wells. Our analysis indicates that the fault zone is a low-permeability feature that inhibits lateral movement of groundwater and that there is at least an order of magnitude decrease in horizontal hydraulic conductivity in the fault core compared to the undeformed sandstone. Large head declines have occurred north of the fault zone (over 300 m since predevelopment conditions) and modifying fault zone parameters significantly affects calibration to regional drawdown on a decadal scale. The flow-barrier behavior of the fault zone has important implications for future groundwater availability in this highly stressed region.     

Flow in an aquifer charged with hot water from a fault zone

Many geothermal anomalies are intersected by vertical fault zones (narrow zones of fractured material with large effective permeability). These conduits are probably responsible for much of the upwelling of hot water from depth. This paper considers a shallow aquifer intersected by a vertical fault. The fluid flow in the aquifer is numerically modeled as a two-dimensional problem. It is observed that the temperature distribution in the aquifer is governed primarily by lateral flow of hot water supplied from the intersecting vertical fault and only secondarily by conduction. The numerical results also provide a possible explanation for the local temperature maxima and inversions occasionally observed in borehole measurements. The present model is an alternative to that based on mushroom-shaped isotherm distributions found in high Rayleigh number large-scale circulation cell calculations.     

Combined electromagnetic, geochemical and thermal surveys of Taal volcano (Philippines) during the period 2005–2006

Since 1572, 33 phreatic to phreatomagmatic eruptions have occurred on Taal volcano (Philippines), some of them causing several hundred casualties. Considering the time delay between two consecutive eruptions, there is an 88% probability that Taal volcano should have already erupted. Since 1992, several phases of seismic activity have been recorded accompanied by ground deformation, opening of fissures, and surface activity. The volcanic activity of Taal appears to be controlled by dike injections and magma supply, buffered by a hydrothermal system that releases fluids and heat through boiling and subsequent steaming. In early 2005, a multidisciplinary project was launched for studying the hydrothermal activity. To map the hydrothermal system, combined surveys were carried out to investigate self-potential, total magnetic field, ground temperature and carbon dioxide soil degassing, along with satellite thermal imaging of the Main Crater Lake. The elevated temperatures and high concentrations of carbon dioxide, as well as electromagnetic anomalies, indicate large-scale hydrothermal degassing. This process is enhanced along the tectonic features (e.g., crater rim and faults) of the volcano, while active fissures opened along the E–W northern flank during the 1992–1994 seismic activity. Heat and fluids from the hydrothermal system are essentially released in the northern part of the crater, which is bounded to the South by a suspected NW–SE fault along which seismicity seems to take place, and dikes are thought to be intruded. During the January 2005 surveys, a new seismic crisis started, and the felt earthquakes prompted spontaneous evacuation of hundreds of inhabitants living on the volcano. Repeated surveys show changes of self-potential, total magnetic field, and ground temperature with time, without any noticeable spatial enlargement. These observations suggest that the northern flank located between the crater rim and the 1992–1994 fissures is connected with a deep thermal source in Main crater and is reactivated during seismic crises. This sector could be subjected to flank failure.     

SPECTRAL STRUCTURE OF VELOCITY INHOMOGENEITY OF CRUST MEDIUM BELOW THE SOUTHEASTERN MARGIN OF TIBETAN PLATEAU AND ITS ADJACENT REGIONS

In this paper we present results of spectral structure of crustal velocity inhomogeneity beneath the southeastern margin of Tibetan plateau and its adjacent region based on the S wave envelope broadening algorithm. The spectral structure of 8~16Hz band is selected to analyze the special character of crustal inhomogeneity and discuss the correlation between strong earthquakes and inhomogeneities. The result shows that strong and complex inhomogeneities of crustal medium are found in the southeastern margin of Tibetan plateau and its adjacent region. In the upper part of upper crust, the strong and small scale inhomogeneities are imaged in the Longmengshan fault zone and the north of the Anninghe fault zone, the weak and large scale inhomogeneites are imaged in the section from Huolu to Daofu of Xianshuihe fault zone and the south of the Anninghe fault zone. In the lower part of upper crust, strong inhomogeneites are found in the Longmengshan fault zone, Lianfeng fault zone, the north of the Anninghe fault zone and the sections from Huolu to Daofu of the Xianshuihe fault zone, weak inhomogeneites are found in the section from Daofu to Kangding of Xianshuihe fault zone. In the middle crust, strong inhomogeneities are observed in the section of the Baoxing to Dujiangyan, the Baoxing to Kangding, and Kangding to Shimian, and weak inhomogeneities are observed in the northwestern section from Huolu to Kangding, and the Lianfeng fault zone. Comparing the medium inhomogeneities with the location of the strong earthquakes, our results suggest existence of high correlation between them. Strong earthquakes are often located in the transitionary zone between the strong and the weak inhomogeneities. The spatial distribution of the strong and the weak medium inhomogeneities may be related to the broken medium from the strong movement of geological tectonic and the heat flow upwelling along active faults induced by frequent tectonic and volcanic activity.