Advances in seismological methods for characterizing fault zone structure

Large earthquakes frequently occur along complex fault systems.Understanding seismic rupture and long-term fault evolution requires constraining the geometric and material properties of fault zone structures.We provide a comprehensive overview of recent advancements in seismological methods used to study fault zone structures,including seismic tomography,fault zone seismic wave analysis,and seismicity analysis.Observational conditions limit our current ability to fully characterize fault zones,for example,insufficient imaging resolution to discern small-scale anomalies,incomplete capture of crucial fault zone seismic waves,and limited precision in event location accuracy.Dense seismic arrays can overcome these limitations and enable more detailed investigations of fault zone structures.Moreover,we present new insights into the structure of the Anninghe-Xiaojiang fault zone in the southeastern margin of the Qinghai-Xizang Plateau based on data collected from a dense seismic array.We found that utilizing a dense seismic array can identify small-scale features within fault zones,aiding in the interpretation of fault zone geometry and material properties.

Detailed sedimentary structure of the Mianning segment of the Anninghe fault zone revealed by H/V spectral ratio

The Anninghe fault is a major left-lateral strike-slip fault in southwest China and a seismic gap with a potential earthquake larger than MW 7.0 lies in the Mianning-Xichang segment according to recent observations.The shallow structure of this region can offer a glimpse into the geometry of the fault,which plays an important role in earthquake hazard mitigation.To further investigate the sedimentary structure of the Anninghe fault zone,two dense linear arrays with a station spacing of around 80 m were deployed across the fault.In this study,the H/V spectral ratio(HVSR),together with its peak frequency at each station site,was obtained by applying the Nakamura method.Our findings demonstrate that the peak frequency behaves in high correlation with lithology and is controlled by topography.HVSR in foothills or regions with magmatic intrusion shows a single peak at about 2–3 Hz.In locations with abundant Quaternary sedimentation,such as Anninghe valleys and fracture zones,another low-frequency peak around 0.4 Hz can be noticed in HVSR.By using the empirical relationship,the thickness of the sedimentary layer around the fault fracture zone is estimated to be 300–600 m.Furthermore,the sedimentary interface shows a downward dip to the east,possibly influenced by the east-west extrusion stress.Considering the resonance effect,buildings with 6–9 stories in the valley area of the Anninghe require additional attention in earthquake hazard prevention.

EVOLUTIONARY MODEL OF FREE ECONOMIC ZONES——Different Generations and Structural Features

Free economic zone (FEZ) has a long history and plays a more and more important role in the world economy. Most studies, however, focused on the theoretical analysis of benefit and cost as well as the economic role of FEZ in the less developed countries and little attention has been paid to the evolution of FEZ. This paper will improve the above-mentioned studies and put forward the structural and spatial evolutionary model of FEZ by analyzing the development of objectives, preferential policy, governance structure, industrial sectors and location of FEZs based on the international economic and political development. FEZs develop towards: 1) more comprehensive and macro objectives, 2) more industry-oriented and multi-preferential policies, 3) more cross-national and combination zones with administrative areas, 4) more technology-intensive and multi-industries, 5) more flexible location and larger spatial dimensions, 6) more rapid evolution and typologies, and 7) more economic integration to the host economy.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology.In this study,we generated a time-calibrated phylogeny of Scutiger boulengeri,an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang(Tibet)Plateau(QTP).Based on three mitochondrial DNA(mt DNA)genes,eight clades were assigned to three deeply divergent lineages.Analysis of nuclear DNA(nu DNA)genes revealed three distinct clusters without geographic structure,indicating significantly high rates of gene flow.Coalescent theory framework analysis(approximate Bayesian computation model DIYABC and Migrate-N)suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago(Ma).The ratio of mt DNA FST(fixation index)to nu DNA FST was 2.3,thus failing to show male-biased dispersal.Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China,without significant reproductive isolation but with strong introgression in S.boulengeri,suggesting high hybrid fitness.Furthermore,mt DNA genes exhibited isolation by distance(IBD)while nu DNA genes exhibited significant isolation by environment(IBE).Results suggested that mitonuclear discordance may have initially been caused by geographic isolation,followed by precipitation-mediated hybridization,producing a wide hybrid zone and geographic structure confusion of nu DNA genes in S.boulengeri.This study indicated that complicated historical processes may have led to specific genetic patterns,with a specific climate factor facilitating gene flow in the system.

Structural patterns of fault damage zones in carbonate rocks and their influences on petroleum accumulation in Tazhong Paleo-uplift, Tarim Basin, NW China

Based on the outcrop survey,3D seismic data interpretation,drilling data analysis,the structural patterns and distribution of fault damage zones in carbonate strata of Tazhong Paleo-uplift were established to reveal the oil and gas enrichment law in the fault damage zones.The following findings were reached:(1)Through the filed survey,the fault damage zone system consists of fault core,damage zone with branch fault and fracture network.Affected by the active nature of the major faults,the fault damage zones differ in planar pattern and scale along the major faults.(2)3D seismic profiles reveal that there are three types of fault damage zones in carbonate strata in Tazhong paleo-uplift,strike-slip fault damage zones,thrust fault damage zones and superimposed fault damage zones.Featuring3 flowers and 3 root belts in vertical,the strike-slip fault damage zone can be subdivided into linear type,oblique type,feather type and horsetail type in plane.Thrust fault damage zones can be further divided into fault anticline type,anticline type and slope type.As the superimposition result of the above two kinds of fault damage zones,superimposed fault damage zones appear in three patterns,intersect type,encompassment type and penetrating type.(3)Cores from wells and geochemical data show oil and gas may migrate along the major fault and laterally.The feather type in strike-slip fault system,fault anticline type in thrust fault damage zone and intersect type in superimposed fault damage zone are possible sites for high production and efficiency wells.

Relation between structural evolution of the Longmenshan orogenic zone and sedimentation of its foreland basin

In order to determine the area for oil and gas exploration in China’s north Sichuan basin,we have divided the time during which the Longmenshan foreland basin was formed into five periods,based on the sedimentary response relationship of the foreland basin to structural evolution: 1) a late Triassic Noric period;2) an early-Middle Jurassic period;3) a late Jurassic to early Cretaceous period;4) a late Cretaceous to Paleogene-Neogene period and 5) the Quaternary period. As well,we analyzed the sedimentary environment and lithologic features of every basin-forming period. The results show that there are several favorable source-reservoir-cap assemblages in our study area,making it a major region for future oil and gas exploration in China’s northern Sichuan basin.

Structural,Kinematic and Geochronological Study of the Main Tianshan Shear Zone in the Borohoro Ranges,NW China

The western segment of the Main Tianshan shear zone(MTSZ)branches into three major dextral strike-slip shear zones in the Borohoro Ranges,northwestern China.Mylonitic rocks within these zones exhibit vertical foliations,subhorizontal lineations,and dextral kinematics.The microstructures and CPO patterns of calcite and quartz in the tectonites suggest shearing at temperatures of 200℃ to 500℃ or above.The amount of dextral displacement along the zones is less than 100 km.LA-ICP-MS U-Pb dating of zircons from some pre-kinematic granites gave ages of 364-325 Ma and 316 Ma,and zircons from a syn-kinematic leucogranite gave an age of 291 Ma.We propose that dextral shearing along the MTSZ started in the earliest Permian.Our data indicate that collision between the Yili and Junggar plates might have finished before 325 Ma or 316 Ma,with a tectonic transformation from 325 Ma to 300 Ma,resulting in intracontinental transpression from 300 Ma or 291 Ma.We suggest that the formation of the Kazakhstan orocline(KO)took place in two stages,with the early stage of bending driven by convergence related to oblique subduction,and the late stage possibly associated with lateral disarticulation as a result of the dextral strike-slip shear zones.

Structural failure mechanism and strengthening method of fracture plugging zone for lost circulation control in deep naturally fractured reservoirs

Focused on the lost circulation control in deep naturally fractured reservoirs, the multiscale structure of fracture plugging zone is proposed based on the theory of granular matter mechanics, and the structural failure pattern of plugging zone is developed to reveal the plugging zone failure mechanisms in deep, high temperature, high pressure, and high in-situ stress environment. Based on the fracture plugging zone strength model, key performance parameters are determined for the optimal selection of loss control material(LCM). Laboratory fracture plugging experiments with new LCM are carried out to evaluate the effect of the key performance parameters of LCM on fracture plugging quality. LCM selection strategy for fractured reservoirs is developed. The results show that the force chain formed by LCMs determines the pressure stabilization of macro-scale fracture plugging zone. Friction failure and shear failure are the two major failure patterns of fracture plugging zone. The strength of force chain depends on the performance of micro-scale LCM, and the LCM key performance parameters include particle size distribution, fiber aspect ratio, friction coefficient, compressive strength, soluble ability and high temperature resistance. Results of lab experiments and field test show that lost circulation control quality can be effectively improved with the optimal material selection based on the extracted key performance parameters of LCMs.

Structural formation and evolution mechanisms of fracture plugging zone

A coupled CFD-DEM method is used to simulate the formation process of fracture plugging zone.A photo-elastic system characterizing mesoscale force chain network developed by our own is used to model the pressure evolution in fracture plugging zone to reveal the evolution mechanism of the structure of fracture plugging zone.A theoretical basis is provided for improving the lost circulation control effect in fractured reservoirs and novel methods are proposed for selecting loss control materials and designing loss control formula.CFD-DEM simulation results show that bridging probability is the key factor determining the formation of fracture plugging zone and fracture plugging efficiency.Critical and absolute bridging concentrations are proposed as the key indexes for loss control formula design.With the increase of absolute bridging concentration,the governing factor of bridging is changed from material grain size to the combination of material grain size and friction force.Results of photo-elastic experiments show that mesoscale force chain network is the intrinsic factor affecting the evolution of pressure exerting on the fracture plugging zone and determines the macroscopic strength of fracture plugging zone.Performance parameters of loss control material affect the force chain network structure and the ratio of stronger force chain,and further impact the stability and strength of fracture plugging zone.Based on the study results,the loss control formula is optimized and new-type loss control material is designed.Laboratory experiments results show that the fracture plugging efficiency and strength is effectively improved.

Comparison of structural diversity of tree-crop associations in Peripheral and Buffer zones of Gachabari Sal forest area,Bangladesh

The structural diversity of different tree-crop associations were studied at Gachabari Sal forest area of Madhupur Garh on Buffer and Peripheral Zone during 2006. The total density, basal area of trees in the Buffer and Peripheral Zone were 155.5 trees·hm^(-2), 795.4 trees·hm^(-2) and 3.9 m2·hm-2, 5.8 m^2·hm^(-2), respectively. No regeneration and natural trees were found in Peripheral Zone and the Zone is totally occupied by exotic species where the Buffer Zone comprised of both natural and exotic trees. The Peripheral Zone belonged to younger and smaller trees whereas the Buffer Zone belonged to mixture of smaller, taller, younger and mature trees simultaneously. For the practicing of different agroforestry systems both Zones have lost their original characters of Sal forest.

Comparison of structural diversity of tree-crop associations in Periph- eral and Buffer zones of Gachabari Sal forest area, Bangladesh

The original version of this article unfortunately contained three mistakes. （1） The spelling of the [A. Nishad＇] name was incorrect. The corrected spelling is given below.

Structural Relationships along a Neoarchean Arc-Continent Collision Zone, North China Craton

The Archean North China Craton is composed of the Western Block,Eastern Block and the intervening Central Orogenic Belt.A 4-10 km wide and 85 km long tectonic mélange belt informally called the Zanhuang tectonic

Estimation of main rheological parameters for Pangxidong-Jinkeng structural fracture zone and Qinzhou-Hangzhou metallogenetic belt in South China

The mylonites occurred in the fracture zones are studied by dynamically recrystallized quartz grains.The natural microstructures in mylonites are simulated and the deformation conditions of mylonitization are estimated by fractal analysis,recrystallized grain size paleopiezometer and flow laws of quartzite.Depending on fractal analysis,the deformation temperature of mylonitization is approximately 600℃,which presents high greenschist facies to low amphibolite facies.The mylonitization occurred at differential stresses of 9.1--10.7MPa(lower limits).Compared with extrapolation of quartzite flow laws and estimates of fractal analysis,the strain rate of mylonitization is under 10-13.8/s.

PCR-DGGE Analysis of Bacterial Communities Structure in Babylonia areolata Culture Systems of The Subtidal Zone and The Pond Mulched Plastic Film and Sand in Bottom

To know the bacterial communities structure in Babylonia areolata culture systems and to research and optimize the management pattem of Babylonia areola-ta culture systems of the pond mulched plastic film and sand in bottom, the bacte- rial communities in Babylonia areolata culture systems of the sub-tidal zone and the pond mulched plastic film and sand in bottom were analyzed at molecular level by adopting the denaturing gradient gel electrophoresis （DGGE）. The results indicated that the dominant bacterial communities in Babylonia areolata culture systems of the sub-tidal zone and the pond mulched plastic film and sand in bottom, which were built on the basis of the seawater in East-island of Zhanjiang, included Proteobac- teda Chloroflexi, Cyanobacteria and Actinobacteria. The dominant bacterial groups in the above pond culture system were Garnmaproteobacteria, Alphaproteobacteria, Deltaprotecbacteda, Epsilonproteobacteda, Anaerolineae, Cyanobacteria and Acti- nobacteda. The dominant bacterial communities in the subtidal zone culture system were Gammaprotecbacteda, Alphaproteobacteria, Deltaproteobacteria, Anaerolineae and Cyanobacteda, and there were less Epsilonproteobacteria and Actinobacteria in the culture system. The higher diversity was detected in the above two culture sys- tems. The results of unweighted pair group method with arithmetic average （UPG- MA） showed that the bacterial communities of the sediment samples S1 and S2 in the above two culture systems were a cluster, the similarity of bacterial communities was 54.5%. The bacterial communities of seawater samples S3 and S4 in the above culture systems were in clusters, and the similarity of the bacterial communi- ties was 84.0%. The results showed that the microorganism ecological level in the Babylonia areolata culture systems of the pond mulched plastic film and sand in bottom could be similar to the sub-tidal zone culture systems through changing the pond seawater and monitoring the microbial population.

Crustal structure in Xiaojiang fault zone and its vicinity

Based on the integrative interpretation of travel-time data and amplitude information obtained from the deep seismic sounding experiment on the Chuxiong-Luoping profile, eastern Yunnan province, carried out in January of 2005, we present a 2-D P wave velocity structure along the profile. The crustal structure shows remarkable contrasts between the two sides of the Xiaojiang fault zone, although the whole profile is situated within the Yangtze platform. The average P wave velocities of the crust on the west and east sides of the fault zone are 6.21 km/s and 6.32 km/s, respectively, and the crustal thicknesses are 41 km and 45 km, respectively. These results imply that the crust to the east of the Xiaojiang fault zone presents characteristics of crustal structure in a stable platform, while the crust to the west is complicated with a lower velocity zone in middle of the upper crust. The average velocity of 6.21 km/s is lower than the global continental crustal average （6.30 km/s）, indicating that the region is tectonically active. According to the lateral variation of velocity and depth of interfaces （including the Moho）, it is inferred that the Xiaojiang fault zone has cut through the whole crust. It is also deduced that existence of low velocity zone in middle of the upper crust is conducive to the south-southeastern sliding of the Sichuan- Yunnan （Chuan-Dian） rhombus block.

Rock Damage Structure of the South Longmen-Shan Fault in the 2008 M8 Wenchuan Earthquake Viewed with Fault-Zone Trapped Waves and Scientific Drilling

This article is to review results from scientific drilling and fault-zone trapped waves （FZTWs） at the south Longman-Shan fault （LSF） zone that ruptured in the 2008 May 12 M8 Wenchuan earthquake in Sichuan,China.Immediately after the mainshock,two Wenchuan Fault Scientific Drilling （WFSD） boreholes were drilled at WFSD-1 and WFSD-2 sites approximately 400 m and 1 km west of the surface rupture along the Yinxiu-Beichuan fault （YBF）,the middle fault strand of the south LSF zone.Two boreholes met the principal slip of Wenchuan earthquake along the YBF at depths of 589-m and 1230-m,respectively.The slip is accompanied with a 100-200-m-wide zone consisting of fault gouge,breccia,cataclasite and fractures.Close to WFSD-1 site,the nearly-vertical slip of ～4.3-m with a 190-m wide zone of highly fractured rocks restricted to the hanging wall of the YBF was found at the ground surface after the Wenchuan earthquake.A dense linear seismic array was deployed across the surface rupture at this venue to record FZTWs generated by aftershocks.Observations and 3-D finite-difference simulations of FZTWs recorded at this cross-fault array and network stations close to the YBF show a distinct low-velocity zone composed by severely damaged rocks along the south LSF at seismogenic depths.The zone is several hundred meters wide along the principal slip,within which seismic velocities are reduced by ～30-55％ from wall-rock velocities and with the maximum velocity reduction in the ～200-m-wide rupture core zone at shallow depth.The FZTW-inferred geometry and physical properties of the south LSF rupture zone at shallow depth are in general consistent with the results from petrological and structural analyses of cores and well log at WFSD boreholes.We interpret this remarkable low-velocity zone as being a break-down zone during dynamic rupture in the 2008 M8 earthquake.We examined the FZTWS generated by similar earthquakes before and after the 2008 mainshock and observed that seismic velocities within fault core zone was reduced by ～10％ due to severe damage of fault rocks during the M8 mainshock.Scientific drilling and locations of aftershocks generating prominent FZTWs also indicate rupture bifurcation along the YBF and the Anxian-Guangxian fault （AGF）,two strands of the south LSF at shallow depth.A combination of seismic,petrologic and geologic study at the south LSF leads to further understand the relationship between the fault-zone structure and rupture dynamics,and the amplification of ground shaking strength along the low-velocity fault zone due to its waveguide effect.

Crustal structure beneath the Qilian Orogen Zone from multiscale seismic tomography

The Qilian Orogen Zone(QOZ), located in the north margin of the Tibetan Plateau, is the key area for understanding the deformation and dynamics process of Tibet. Numerous geological and geophysical studies have been carried out on the mechanics of the Tibetan Plateau deformation and uplift; however, the detailed structure and deformation style of the Qilian Orogen Zone have remained uncertain due to poor geophysical data coverage and limited resolution power of inversion algorithms. In this study, we analyze the P-wave velocity structure beneath the Qilian Orogen Zone, obtained by applying multi-scale seismic tomography technique to P-wave arrival time data recorded by regional seismic networks. The seismic tomography algorithm used in this study employs sparsity constraints on the wavelet representation of the velocity model via L1-norm regularization. This algorithm can deal efficiently with uneven-sampled volumes, and can obtain multi-scale images of the velocity model. Our results can be summarized as follows:(1) The crustal velocity structure is strongly inhomogeneous and consistent with the surface geological setting. Significant low-velocity anomalies exist in the crust of northeastern Tibet, and slight high-velocity anomalies exist beneath the Qaidam Basin and Alxa terrane.(2)The Qilian Orogen Zone can be divided into two main parts by the Laji Shan Faults: the northwestern part with a low-velocity feature, and the southeastern part with a high-velocity feature at the upper and middle crust.(3) Our tomographic images suggest that northwestern and southeastern Qilian Orogen Zones have undergone different tectonic processes. In the northwest Qilian Orogen Zone, the deformation and growth of the Northern Tibetan Plateau has extended to the Heli Shan and Beida Shan region by northward overthrusting at the upper crust and thickening in the lower crust. We speculate that in the southeast Qilian Orogen Zone the deformation and growth of the Northern Tibet Plateau were of strike-slip style at the upper crust; in the lower crust, the evidence suggests ductile shear extrusion style and active frontage extension to the Alxa terrane.(4) The multi-scale seismic tomography technique provides multiscale analysis and sparse constraints, which has allowed to us obtain stable, high-resolution results.

Genetic models of structural traps related to normal faults in the Putaohua Oilfield,Songliao Basin

The Putaohua Oilfield is a fault-prolific area and the faults have close relation with structural traps. The genetic models of the structural traps in the Putaohua Oilfield can be divided into two types： individual fault model and multi-fault interaction model. This is based on the description of displacement distribution of typical individual normal faults, the geometry of the footwall and hanging wall, and the analysis of the interaction between faults and the corresponding change in geometry when the faults grow. The individual fault model is that the displacement reaches a maximum at or near the center of fault and decreases toward the fault tips, so a half-graben is formed on the hanging wall of the fault and a half- anticline is formed on the footwall because of the isostatic process. The multi-fault interaction model is that during the growth of faults, they overlap and interact with each other, and accommodation zones are formed in the overlapping segments. The accommodation zones are favorable targets for hydrocarbon exploration, and the trap characteristics are dependent The multi-fault interaction model can be subdivided on the extent of overlap and occurrence of faults. into three types： synthetic accommodation zone, convergent accommodation zone and divergent accommodation zone. Hydrocarbon migration and accumulation models of each type have been developed. The hydrocarbon migration and accumulation models of the traps with different genetic models have their own characteristics in the different stages of fault growth.

Internal structures and high-velocity frictional properties of Longmenshan fault zone at Shenxigou activated during the 2008 Wenchuan earthquake

This paper reports internal structures of a wide fault zone at Shenxigou,Dujiangyan,Sichuan province,China,and high-velocity frictional properties of the fault gouge collected near the coseismic slip zone during the 2008 Wenchuan earthquake.Vertical offset and horizontal displacement at the trench site were 2.8 m（NW side up）and 4.8 m（right-lateral）,respectively.The fault zone formed in Triassic sandstone,siltstone,and shale about 500 m away from the Yingxiu-Beichuan fault,a major fault in the Longmenshan fault system.A trench survey across the coseismic fault,and observations of outcrops and drill cores down to a depth of 57 m revealed that the fault zone consists of fault gouge and fault breccia of about0.5 and 250-300 m in widths,respectively,and that the fault strikes N62°E and dips 68° to NW.Quaternary conglomerates were recovered beneath the fault in the drilling,so that the fault moved at least 55 m along the coseismic slip zone,experiencing about 18 events of similar sizes.The fault core is composed of grayish gouge（GG） and blackish gouge（BG） with very complex slip-zone structures.BG contains low-crystalline graphite of about 30 %.High-velocity friction experiments were conducted at normal stresses of 0.6-2.1 MPa and slip rates of 0.1-2.1 m/s.Both GG and BG exhibit dramatic slip weakening at constant high slip rates that can be described as an exponential decay from peak friction coefficient lpto steadystate friction coefficient lssover a slip-weakening distance Dc.Deformation of GG and BG is characterized by overlapped slip-zone structures and development of sharp slickenside surfaces,respectively.Comparison of our data with those reported for other outcrops indicates that the high-velocity frictional properties of the Longmenshan fault zones are quite uniform and the high-velocity weakening must have promoted dynamic rupture propagation during the Wenchuan earthquake.

3-D shear wave velocity structure in the shallow crust of the Tan-Lu fault zone in Lujiang,Anhui,and adjacent areas,and its tectonic implications

The Tan-Lu fault zone is a large NNE-trending fault zone in eastern China.Investigations of the structures of the fault zone and its surrounding areas have attracted much attention.In this study,we used dense-array ambient noise tomography to construct a threedimensional shear wave velocity model of shallow crust in an area about 80km×70km in Lujiang,Anhui Province,eastern China.For approximately one month we collected continuous ambient noise signals recorded by 90 short-period seismographs in the region,and obtained the short-period Rayleigh wave empirical Green's functions between stations by the cross-correlation method;we also extracted 0.5–8 s fundamental mode Rayleigh wave group velocity and phase velocity dispersion curves.Based on the direct surface wave tomography method,we jointly inverted the group velocity and phase velocity dispersion data of all paths and obtained the 3-D shear wave velocity structure in the depth range of 0–5 km.The results revealed important geological structural features of the study area.In the north region,the sedimentary center of the Hefei Basin—the southwestern part of the Chaohu Lake—shows a significant low-velocity anomaly to a depth of at least 5 km.The southwestern and southeastern regions of the array are the eastern margin of the Dabie orogenic belt and the intrusion area of Luzong volcanic rocks,respectively,and both show obvious high-speed anomalies;the sedimentary area within the Tan-Lu fault zone(about 10 km wide)shows low-velocity anomalies.However,the volcanic rock intrusion area in the fault zone is shown as high velocity.Our shallow crustal imaging results reflect the characteristics of different structures in the study area,especially the high-speed intrusive rocks in the Tan-Lu fault zone,which were probably partially derived from the magmatic activity of Luzong volcanic basin.From the Late Cretaceous to Early Tertiary,the Tan-Lu fault zone was in a period of extensional activity;the special stress environment and the fractured fault zone morphology provided conditions for magma in the Luzong volcanic basin to intrude into the Tan-Lu fault zone in the west.Our 3-D model can also provide important information for deep resource exploration and earthquake strong ground motion simulation.