作为EMSLAB电磁感应国际合作计划的一部分,通过对海底观测的三分量变化地磁场的分析,构制了横跨Juan de Fuca洋脊、Juan de Fuca板块到板块俯冲带的三条地磁转换函数剖面,计算了各个测点的地磁场感应矢量。结果表明,Juan de Fuca板块的...作为EMSLAB电磁感应国际合作计划的一部分,通过对海底观测的三分量变化地磁场的分析,构制了横跨Juan de Fuca洋脊、Juan de Fuca板块到板块俯冲带的三条地磁转换函数剖面,计算了各个测点的地磁场感应矢量。结果表明,Juan de Fuca板块的东、西两部分的浅部电性结构是有差异的,近海岸线的东区更导电;而深部电性结构是相对均匀的。这一结果与由大地电磁资料计算得到的二维模型剖面是一致的。展开更多
To properly simulate hard rock with a high ratio of the uniaxial compressive strength to tensile strength(UCS/TS) and realistic strength-failure envelope,the rock deformation and mechanical characteristics were discus...To properly simulate hard rock with a high ratio of the uniaxial compressive strength to tensile strength(UCS/TS) and realistic strength-failure envelope,the rock deformation and mechanical characteristics were discussed in detail when the particle simulation method with the clump parallel-bond model(CPBM) was used to conduct a series of numerical experiments at the specimen scale.Meanwhile,the effects of the loading procedure and crack density on the mechanical behavior of a specimen,which was modeled by the particle simulation method with the CPBM,were investigated.The related numerical results have demonstrated that:1) The uniaxial compressive strength(UCS),tensile strength(TS) and elastic modulus are overestimated when the conventional loading procedure is used in the particle simulation method with the CPBM; 2) The elastic modulus,strength and UCS/TS decrease,while Poisson ratio increases with the increase of the crack density in the particle simulation method with the CPBM; 3) The particle simulation method with the CPBM can be used to reproduce a high value of UCS/TS(>10),as well as a high friction angle and reasonable cohesion strength; 4) As the confining pressure increases,both the peak strength of the simulated specimen and the number of microscopic cracks increase,but the ratio of tensile cracks number to shear cracks number decreases in the particle simulation method with the CPBM; 5) Compared with the conventional parallel-bond model,the CPBM can be used to reproduce more accurate results for simulating the rock deformation and mechanical characteristics.展开更多
The Ailao Shan-Red River fault zone is the boundary between the Yangtze block to the northeast and the Indochina block to the southwest.It is an important tectonic zone due to its role in the southeastward extrusion o...The Ailao Shan-Red River fault zone is the boundary between the Yangtze block to the northeast and the Indochina block to the southwest.It is an important tectonic zone due to its role in the southeastward extrusion of the Indochina block during and subsequent to the Indian-Eurasian collision.Diancang Shan（DCS） high-grade metamorphic complex,located at the northwest extension along the Ailao Shan-Red River（ASRR） shear zone,is a representative metamorphic complex of the ASRR tectonic belt.Structural and microstructural analysis of sheared rocks in the high-grade metamorphic rocks reveals that they are coherent with solid-state high-temperature ductile deformation,which is attributed to left-lateral shearing along the ASRR shear zone.New LA-ICP-MS zircon U-Pb geochronological and microstructural studies of the post-kinematic granitic plutons provide a straightforward time constraint on the termination ductile left-lateral shearing and exhumation of the metamorphic massif in the ASRR shear zone.It is suggested that the left-lateral shearing along the ASRR shear zone ended at ca.21 Ma at relative lower-temperature or decreasing temperature conditions.During or after the emplacement of the young dikes at ca.21 Ma,rapid brittle deformation event occurred,which makes the DCS massif start fast uplift/exhumation and cooling to a shallow crustal level.展开更多
Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation m...Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation method was presented to produce meaningful numerical solutions for static problems. The central part of the explicit dynamic relaxation method is to turn a time-independent static problem into an artificial time-dependent dynamic problem. The related numerical testing results demonstrate that: (1) the proposed effective loading algorithm is capable of enabling an applied load in a static problem to be propagated throughout the whole system within a given loading increment, so that the time-independent solution of the static problem can be obtained; (2) the proposed effective loading algorithm can be straightforwardly applied to the particle simulation method for solving a wide range of static problems.展开更多
Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process...Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.展开更多
Seven conodont biozones are recognized in the carbonate-dominated shelf-marine Middle Ordovician developed in the intracratonic sedimentary basins(Canning,Amadeus and Georgina)of central and north-western Australia,in...Seven conodont biozones are recognized in the carbonate-dominated shelf-marine Middle Ordovician developed in the intracratonic sedimentary basins(Canning,Amadeus and Georgina)of central and north-western Australia,in the Lachlan and New England orogens of New South Wales,and in the Takaka Terrane of New Zealand.A separate scheme identifying seven conodont biozones spanning the Middle Ordovician has also been developed for siliciclastic sequences deposited in slopebasinal environments in the Lachlan Orogen in New South Wales and Victoria.This biozonal classification consisting of two parallel biostratigraphic schemes for the shelf-marine and deep-marine successions respectively has significantly increased precision in regional and global biostratigraphic correlation and laid a solid foundation for the Middle Ordovician chronostratigraphy of Australia and New Zealand.Recognition of short-ranging pandemic species as the eponymous species of the biozones also supports direct correlation with the classical conodont successions established in Baltoscandia and the North American Midcontinent,and with those of the major Chinese terranes(South and North China and Tarim).The Lachlan Orogen appears to be globally unique in enabling correlation of contemporaneous conodont faunas over a considerable spectrum of water depths and biofacies ranging from carbonate shelves,slopes to deep-water basins.展开更多
Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-fron...Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-front instability problem in fluid-saturated porous rocks is no exception.Since this kind of instability problem has both the conventional(i.e.trivial)and the unconventional(i.e.nontrivial)solutions,it is necessary to examine the effects of different numerical algorithms,which are used to solve chemical dissolution-front instability problems in fluid-saturated porous rocks.Toward this goal,two different numerical algorithms associated with the commonly-used finite element method are considered in this paper.In the first numerical algorithm,the porosity,pore-fluid pressure and acid/solute concentration are selected as basic variables,while in the second numerical algorithm,the porosity,velocity of pore-fluid flow and acid/solute concentration are selected as basic variables.The particular attention is paid to the effects of these two numerical algorithms on the computational simulation results of unstable chemical dissolution-front propagation in fluid-saturated porous rocks.The related computational simulation results have demonstrated that:1)the first numerical algorithm associated with the porosity-pressure-concentration approach can realistically simulate the evolution processes of unstable chemical dissolution-front propagation in chemical dissolution systems.2)The second numerical algorithm associated with the porosity-velocity-concentration approach fails to simulate the evolution processes of unstable chemical dissolution-front propagation.3)The extra differential operation is the main source to result in the failure of the second numerical algorithm.展开更多
Ordovician conodonts representing 28 genera and 28 named and three unnamed species were identified from 740 chert and siliceous siltstone spot samples(>3000 thin sections)from deep-water turbiditic sequences of the...Ordovician conodonts representing 28 genera and 28 named and three unnamed species were identified from 740 chert and siliceous siltstone spot samples(>3000 thin sections)from deep-water turbiditic sequences of the Lachlan Orogen in central and southern New South Wales,Australia.Based on these faunas,a new conodont biozonal scheme has been established to divide the Ordovician turbiditic successions of the Lachlan Orogen into 12 superbiozones and biozones.They are(in ascending order)the Paracordylodus gracilis Superbiozone(including the Prioniodus oepiki Biozone),Periodon flabellum Superbiozone(including the Oepikodus evae Biozone in the lower part),Periodon hankensis Biozone,Periodon aculeatus Superbiozone(including the Histiodella labiosa,Histiodella holodentata,Histiodella kristinae,Pygodus serra and Pygodus anserinus biozones)and the Periodon grandis Biozone.The Pygodus anserinus Biozone is divided further into the lower and upper subbiozones.This new conodont biozonation scheme spanning the upper Tremadocian to middle Katian interval permits precise age-dating and correlation of deep-water siliciclastic rocks that characterize the Ordovician Deep-Sea Realm regionally and internationally.展开更多
Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys ...Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys the random uniform distribution,both the mean value and the variance of porosities in the heterogeneous porosity field can be used to reflect the overall heterogeneous characteristics of the porosity field.The main purpose of this work is to investigate the effects of porosity heterogeneity on chemical dissolution front instability in fluid-saturated rocks by the computational simulation method.The related computational simulation results have demonstrated that:1) since the propagation speed of a chemical dissolution front is inversely proportional to the difference between the final porosity and the mean value of porosities in the initial porosity field,an increase in the extent of the porosity heterogeneity can cause an increase in the mean value of porosities in the initial porosity field and an increase in the propagation speed of the chemical dissolution front.2) An increase in the variance of porosities in the initial porosity field can cause an increase in the instability probability of the chemical dissolution front in the fluid-saturated rock.3) The greater the mean value of porosities in the initial porosity field,the quicker the irregular morphology of the chemical dissolution front changes in the supercritical chemical dissolution systems.This means that the irregular morphology of a chemical dissolution front grows quicker in a porosity field of heterogeneity than it does in that of homogeneity when the chemical dissolution system is at a supercritical stage.展开更多
To look for gas hydrate, 22 multi-channel and 3 single-channel seismic lines on the East China Sea （ECS） shelf slope and at the bottom of the Okinawa Trough were examined. It was found that there was indeed bottom s...To look for gas hydrate, 22 multi-channel and 3 single-channel seismic lines on the East China Sea （ECS） shelf slope and at the bottom of the Okinawa Trough were examined. It was found that there was indeed bottom simulating reflector （BSR） occurrence, but it is very rare. Besides several BSRs, a gas seepage was also found. As shown by the data, both the BSR and gas seepage are all related with local geological structures, such as mud diapir, anticline, and fault-controlled graben-like structure. However, similar structural ＂anomalies＂ are quite common in the tectonically very active Okinawa Trough region, but very few of them have developed BSR or gas seepage. The article points out that the main reason is probably the low concentration of organic carbon of the sediment in this area. It was speculated that the rare occurrence of gas hydrates in this region is governed by structure-controlled fluid flow. Numerous faults and fractures form a network of high-permeability channels in the sediment and highly fractured igneous basement to allow fluid circulation and ventilation. Fluid flow in this tectonic environment is driven primarily by thermal buoyancy and takes place on a wide range of spatial scales. The fluid flow may play two roles to facilitate hydrate formation： to help gather enough methane into a small area and to modulate the thermal regime.展开更多
The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is ...The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.展开更多
Through integrating the state of the art scientific knowledge in different research fields, some potential mechanisms of large-scale movements of underground pore-fluids such as H2O and CO2 in the continental lithosph...Through integrating the state of the art scientific knowledge in different research fields, some potential mechanisms of large-scale movements of underground pore-fluids such as H2O and CO2 in the continental lithosphere were presented and discussed. The results show that the generation and propagation of porosity waves are important mechanisms to transport mass and heat fluxes from the continental lithospheric mantle into the lower continental crust; the generation and propagation of porosity waves, pore-fluid flow focusing through lower and middle crustal faults, advection of pore-fluids through the lower and middle crust, and whole-crust convection in some particular cases are important mechanisms to transport mass and heat fluxes from the lower into the upper continental crust; heat and mass transport through convective pore-fluid flow is the most effective mechanism of ore body formation and mineralization in hydrothermal systems; due to heat and mass exchange at the interface between the earth surface, hydrosphere and atmosphere, it is very important to consider the hydro-geological effect of the deep earth pore-fluids such as H2O and CO2 on the global warming and climate change in future investigations.展开更多
In order to simulate the instability phenomenon of a nonaqueous phase liquid(NAPL) dissolution front in a computational model, the intrinsic characteristic length is commonly used to determine the length scale at whic...In order to simulate the instability phenomenon of a nonaqueous phase liquid(NAPL) dissolution front in a computational model, the intrinsic characteristic length is commonly used to determine the length scale at which the instability of the NAPL dissolution front can be initiated. This will require a huge number of finite elements if a whole NAPL dissolution system is simulated in the computational model. Even though modern supercomputers might be used to tackle this kind of NAPL dissolution problem, it can become prohibitive for commonly-used personal computers to do so. The main purpose of this work is to investigate whether or not the whole NAPL dissolution system of an annular domain can be replaced by a trapezoidal domain, so as to greatly reduce the requirements for computer efforts. The related simulation results have demonstrated that when the NAPL dissolution system under consideration is in a subcritical state, if the dissolution pattern around the entrance of an annulus domain is of interest, then a trapezoidal domain cannot be used to replace an annular domain in the computational simulation of the NAPL dissolution system.However, if the dissolution pattern away from the vicinity of the entrance of an annulus domain is of interest, then a trapezoidal domain can be used to replace an annular domain in the computational simulation of the NAPL dissolution system. When the NAPL dissolution system under consideration is in a supercritical state, a trapezoidal domain cannot be used to replace an annular domain in the computational simulation of the NAPL dissolution system.展开更多
The asthenosphere upwelled on a large scale in the western Pacific and South China Sea during the Cenozoic, which formed strong upward throughflow and caused the thermal structure to be changed obviously. The mathemat...The asthenosphere upwelled on a large scale in the western Pacific and South China Sea during the Cenozoic, which formed strong upward throughflow and caused the thermal structure to be changed obviously. The mathematical analysis has demonstrated that the upward throughflow velocity may have varied from 3×10^11 to 6×10^12 m/s. From the relationship between the lithospheric thickness and the conductive heat flux, the lithospheric heat flux in the western Pacific should be above 30 mW/m^2, which is consistent with the observed data. The huge low-speed zone within the upper mantle of the marginal sea in the western Pacific reflects that the upper mantle melts partially, flows regionally in the regional stress field, forms the upward heat flux at its bottom, and causes the change of the lithospheric thermal structure in the region. The numerical simulation result of the expansion and evolution in the South China Sea has demonstrated that in the early expansion, the upward throughflow velocity was relatively fast, and the effect that it had on the thickness of the lithosphere was relatively great,resulting in the mid-ocean basin expanding rapidly. After the formation of the ocean basin in the South China Sea, the upward throughflow velocity decreased, but the conductive heat flux was relatively high, which is close to the actual situation. Therefore, from the heat transfer point of view, this article discusses how the upward heat flux affects the lithospheric thermal structure in the western Pacific and South China Sea. The conclusions show that the upward heat throughflow at the bottom of the lithospheric mantle resulted in the tectonic deformation at the shallow crust. The intensive uplifts and rifts at the crust led to the continent cracks and the expansion in the South China Sea.展开更多
The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper ...The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.展开更多
This paper presents numerical investigation on the ore-forming fluid migration driven by tectonic deformation and thermally-induced buoyancy force in the Chanziping ore district in South China.A series of numerical sc...This paper presents numerical investigation on the ore-forming fluid migration driven by tectonic deformation and thermally-induced buoyancy force in the Chanziping ore district in South China.A series of numerical scenarios are considered to examine the effect of meteoric water precipitation, the dip angle of the faults,unconformity surface,and thermal input on the ore genesis.Our computations reveal that the downward basinal fluid flow driven by extensional stress mixes with the upward basal fluid driven by the thermal input from depth at the junction of two faults at a temperature of about 200℃,triggering the precipitation of the Chanziping uranium deposit.展开更多
文摘作为EMSLAB电磁感应国际合作计划的一部分,通过对海底观测的三分量变化地磁场的分析,构制了横跨Juan de Fuca洋脊、Juan de Fuca板块到板块俯冲带的三条地磁转换函数剖面,计算了各个测点的地磁场感应矢量。结果表明,Juan de Fuca板块的东、西两部分的浅部电性结构是有差异的,近海岸线的东区更导电;而深部电性结构是相对均匀的。这一结果与由大地电磁资料计算得到的二维模型剖面是一致的。
基金Project(11272359) supported by the National Natural Science Foundation of China
文摘To properly simulate hard rock with a high ratio of the uniaxial compressive strength to tensile strength(UCS/TS) and realistic strength-failure envelope,the rock deformation and mechanical characteristics were discussed in detail when the particle simulation method with the clump parallel-bond model(CPBM) was used to conduct a series of numerical experiments at the specimen scale.Meanwhile,the effects of the loading procedure and crack density on the mechanical behavior of a specimen,which was modeled by the particle simulation method with the CPBM,were investigated.The related numerical results have demonstrated that:1) The uniaxial compressive strength(UCS),tensile strength(TS) and elastic modulus are overestimated when the conventional loading procedure is used in the particle simulation method with the CPBM; 2) The elastic modulus,strength and UCS/TS decrease,while Poisson ratio increases with the increase of the crack density in the particle simulation method with the CPBM; 3) The particle simulation method with the CPBM can be used to reproduce a high value of UCS/TS(>10),as well as a high friction angle and reasonable cohesion strength; 4) As the confining pressure increases,both the peak strength of the simulated specimen and the number of microscopic cracks increase,but the ratio of tensile cracks number to shear cracks number decreases in the particle simulation method with the CPBM; 5) Compared with the conventional parallel-bond model,the CPBM can be used to reproduce more accurate results for simulating the rock deformation and mechanical characteristics.
基金support from the State Key Research"973"Plan of China(No. 2009CB421001)National Natural Science Foundation of China(40872139)+1 种基金the 111 Project(B07011) of the Ministry of Education,State Key Laboratory of Geological Processes and Mineral Resources (GPMR200837)the Fundamental Research Funds for the Central Universities(GPMR2009PY01)
文摘The Ailao Shan-Red River fault zone is the boundary between the Yangtze block to the northeast and the Indochina block to the southwest.It is an important tectonic zone due to its role in the southeastward extrusion of the Indochina block during and subsequent to the Indian-Eurasian collision.Diancang Shan（DCS） high-grade metamorphic complex,located at the northwest extension along the Ailao Shan-Red River（ASRR） shear zone,is a representative metamorphic complex of the ASRR tectonic belt.Structural and microstructural analysis of sheared rocks in the high-grade metamorphic rocks reveals that they are coherent with solid-state high-temperature ductile deformation,which is attributed to left-lateral shearing along the ASRR shear zone.New LA-ICP-MS zircon U-Pb geochronological and microstructural studies of the post-kinematic granitic plutons provide a straightforward time constraint on the termination ductile left-lateral shearing and exhumation of the metamorphic massif in the ASRR shear zone.It is suggested that the left-lateral shearing along the ASRR shear zone ended at ca.21 Ma at relative lower-temperature or decreasing temperature conditions.During or after the emplacement of the young dikes at ca.21 Ma,rapid brittle deformation event occurred,which makes the DCS massif start fast uplift/exhumation and cooling to a shallow crustal level.
基金Projects(10872219 10672190) supported by the National Natural Science Foundation of China
文摘Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation method was presented to produce meaningful numerical solutions for static problems. The central part of the explicit dynamic relaxation method is to turn a time-independent static problem into an artificial time-dependent dynamic problem. The related numerical testing results demonstrate that: (1) the proposed effective loading algorithm is capable of enabling an applied load in a static problem to be propagated throughout the whole system within a given loading increment, so that the time-independent solution of the static problem can be obtained; (2) the proposed effective loading algorithm can be straightforwardly applied to the particle simulation method for solving a wide range of static problems.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.
基金funded by the Geological Survey of New South Wales。
文摘Seven conodont biozones are recognized in the carbonate-dominated shelf-marine Middle Ordovician developed in the intracratonic sedimentary basins(Canning,Amadeus and Georgina)of central and north-western Australia,in the Lachlan and New England orogens of New South Wales,and in the Takaka Terrane of New Zealand.A separate scheme identifying seven conodont biozones spanning the Middle Ordovician has also been developed for siliciclastic sequences deposited in slopebasinal environments in the Lachlan Orogen in New South Wales and Victoria.This biozonal classification consisting of two parallel biostratigraphic schemes for the shelf-marine and deep-marine successions respectively has significantly increased precision in regional and global biostratigraphic correlation and laid a solid foundation for the Middle Ordovician chronostratigraphy of Australia and New Zealand.Recognition of short-ranging pandemic species as the eponymous species of the biozones also supports direct correlation with the classical conodont successions established in Baltoscandia and the North American Midcontinent,and with those of the major Chinese terranes(South and North China and Tarim).The Lachlan Orogen appears to be globally unique in enabling correlation of contemporaneous conodont faunas over a considerable spectrum of water depths and biofacies ranging from carbonate shelves,slopes to deep-water basins.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-front instability problem in fluid-saturated porous rocks is no exception.Since this kind of instability problem has both the conventional(i.e.trivial)and the unconventional(i.e.nontrivial)solutions,it is necessary to examine the effects of different numerical algorithms,which are used to solve chemical dissolution-front instability problems in fluid-saturated porous rocks.Toward this goal,two different numerical algorithms associated with the commonly-used finite element method are considered in this paper.In the first numerical algorithm,the porosity,pore-fluid pressure and acid/solute concentration are selected as basic variables,while in the second numerical algorithm,the porosity,velocity of pore-fluid flow and acid/solute concentration are selected as basic variables.The particular attention is paid to the effects of these two numerical algorithms on the computational simulation results of unstable chemical dissolution-front propagation in fluid-saturated porous rocks.The related computational simulation results have demonstrated that:1)the first numerical algorithm associated with the porosity-pressure-concentration approach can realistically simulate the evolution processes of unstable chemical dissolution-front propagation in chemical dissolution systems.2)The second numerical algorithm associated with the porosity-velocity-concentration approach fails to simulate the evolution processes of unstable chemical dissolution-front propagation.3)The extra differential operation is the main source to result in the failure of the second numerical algorithm.
文摘Ordovician conodonts representing 28 genera and 28 named and three unnamed species were identified from 740 chert and siliceous siltstone spot samples(>3000 thin sections)from deep-water turbiditic sequences of the Lachlan Orogen in central and southern New South Wales,Australia.Based on these faunas,a new conodont biozonal scheme has been established to divide the Ordovician turbiditic successions of the Lachlan Orogen into 12 superbiozones and biozones.They are(in ascending order)the Paracordylodus gracilis Superbiozone(including the Prioniodus oepiki Biozone),Periodon flabellum Superbiozone(including the Oepikodus evae Biozone in the lower part),Periodon hankensis Biozone,Periodon aculeatus Superbiozone(including the Histiodella labiosa,Histiodella holodentata,Histiodella kristinae,Pygodus serra and Pygodus anserinus biozones)and the Periodon grandis Biozone.The Pygodus anserinus Biozone is divided further into the lower and upper subbiozones.This new conodont biozonation scheme spanning the upper Tremadocian to middle Katian interval permits precise age-dating and correlation of deep-water siliciclastic rocks that characterize the Ordovician Deep-Sea Realm regionally and internationally.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘Homogeneity and heterogeneity are two totally different concepts in nature.At the particle length scale,rocks exhibit strong heterogeneity in their constituents and porosities.When the heterogeneity of porosity obeys the random uniform distribution,both the mean value and the variance of porosities in the heterogeneous porosity field can be used to reflect the overall heterogeneous characteristics of the porosity field.The main purpose of this work is to investigate the effects of porosity heterogeneity on chemical dissolution front instability in fluid-saturated rocks by the computational simulation method.The related computational simulation results have demonstrated that:1) since the propagation speed of a chemical dissolution front is inversely proportional to the difference between the final porosity and the mean value of porosities in the initial porosity field,an increase in the extent of the porosity heterogeneity can cause an increase in the mean value of porosities in the initial porosity field and an increase in the propagation speed of the chemical dissolution front.2) An increase in the variance of porosities in the initial porosity field can cause an increase in the instability probability of the chemical dissolution front in the fluid-saturated rock.3) The greater the mean value of porosities in the initial porosity field,the quicker the irregular morphology of the chemical dissolution front changes in the supercritical chemical dissolution systems.This means that the irregular morphology of a chemical dissolution front grows quicker in a porosity field of heterogeneity than it does in that of homogeneity when the chemical dissolution system is at a supercritical stage.
基金the National Natural Science Foundation of China (No.40776032)Key Discipline Program of Chinese Academy of Sciences (No.KZCX2-211-01)National Basic Research Program of China (973 Program,No.207CB411702)
文摘To look for gas hydrate, 22 multi-channel and 3 single-channel seismic lines on the East China Sea （ECS） shelf slope and at the bottom of the Okinawa Trough were examined. It was found that there was indeed bottom simulating reflector （BSR） occurrence, but it is very rare. Besides several BSRs, a gas seepage was also found. As shown by the data, both the BSR and gas seepage are all related with local geological structures, such as mud diapir, anticline, and fault-controlled graben-like structure. However, similar structural ＂anomalies＂ are quite common in the tectonically very active Okinawa Trough region, but very few of them have developed BSR or gas seepage. The article points out that the main reason is probably the low concentration of organic carbon of the sediment in this area. It was speculated that the rare occurrence of gas hydrates in this region is governed by structure-controlled fluid flow. Numerous faults and fractures form a network of high-permeability channels in the sediment and highly fractured igneous basement to allow fluid circulation and ventilation. Fluid flow in this tectonic environment is driven primarily by thermal buoyancy and takes place on a wide range of spatial scales. The fluid flow may play two roles to facilitate hydrate formation： to help gather enough methane into a small area and to modulate the thermal regime.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.
基金Project(10672190) supported by the National Natural Science Foundation of China
文摘Through integrating the state of the art scientific knowledge in different research fields, some potential mechanisms of large-scale movements of underground pore-fluids such as H2O and CO2 in the continental lithosphere were presented and discussed. The results show that the generation and propagation of porosity waves are important mechanisms to transport mass and heat fluxes from the continental lithospheric mantle into the lower continental crust; the generation and propagation of porosity waves, pore-fluid flow focusing through lower and middle crustal faults, advection of pore-fluids through the lower and middle crust, and whole-crust convection in some particular cases are important mechanisms to transport mass and heat fluxes from the lower into the upper continental crust; heat and mass transport through convective pore-fluid flow is the most effective mechanism of ore body formation and mineralization in hydrothermal systems; due to heat and mass exchange at the interface between the earth surface, hydrosphere and atmosphere, it is very important to consider the hydro-geological effect of the deep earth pore-fluids such as H2O and CO2 on the global warming and climate change in future investigations.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘In order to simulate the instability phenomenon of a nonaqueous phase liquid(NAPL) dissolution front in a computational model, the intrinsic characteristic length is commonly used to determine the length scale at which the instability of the NAPL dissolution front can be initiated. This will require a huge number of finite elements if a whole NAPL dissolution system is simulated in the computational model. Even though modern supercomputers might be used to tackle this kind of NAPL dissolution problem, it can become prohibitive for commonly-used personal computers to do so. The main purpose of this work is to investigate whether or not the whole NAPL dissolution system of an annular domain can be replaced by a trapezoidal domain, so as to greatly reduce the requirements for computer efforts. The related simulation results have demonstrated that when the NAPL dissolution system under consideration is in a subcritical state, if the dissolution pattern around the entrance of an annulus domain is of interest, then a trapezoidal domain cannot be used to replace an annular domain in the computational simulation of the NAPL dissolution system.However, if the dissolution pattern away from the vicinity of the entrance of an annulus domain is of interest, then a trapezoidal domain can be used to replace an annular domain in the computational simulation of the NAPL dissolution system. When the NAPL dissolution system under consideration is in a supercritical state, a trapezoidal domain cannot be used to replace an annular domain in the computational simulation of the NAPL dissolution system.
基金supported jointly by the Important Direction Project of the CAS Knowledge Innovation Program (Nos. KZCX2-YW-203-01, KZCX2-YW-128-4)the Program of Key Laboratory of Marginal Sea Geology (No. MSGL04-8)
文摘The asthenosphere upwelled on a large scale in the western Pacific and South China Sea during the Cenozoic, which formed strong upward throughflow and caused the thermal structure to be changed obviously. The mathematical analysis has demonstrated that the upward throughflow velocity may have varied from 3×10^11 to 6×10^12 m/s. From the relationship between the lithospheric thickness and the conductive heat flux, the lithospheric heat flux in the western Pacific should be above 30 mW/m^2, which is consistent with the observed data. The huge low-speed zone within the upper mantle of the marginal sea in the western Pacific reflects that the upper mantle melts partially, flows regionally in the regional stress field, forms the upward heat flux at its bottom, and causes the change of the lithospheric thermal structure in the region. The numerical simulation result of the expansion and evolution in the South China Sea has demonstrated that in the early expansion, the upward throughflow velocity was relatively fast, and the effect that it had on the thickness of the lithosphere was relatively great,resulting in the mid-ocean basin expanding rapidly. After the formation of the ocean basin in the South China Sea, the upward throughflow velocity decreased, but the conductive heat flux was relatively high, which is close to the actual situation. Therefore, from the heat transfer point of view, this article discusses how the upward heat flux affects the lithospheric thermal structure in the western Pacific and South China Sea. The conclusions show that the upward heat throughflow at the bottom of the lithospheric mantle resulted in the tectonic deformation at the shallow crust. The intensive uplifts and rifts at the crust led to the continent cracks and the expansion in the South China Sea.
基金Project(10872219) supported by the National Natural Science Foundation of China
文摘The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.
基金supported by"The Program to Sponsor Teams for Innovation in the Construction of Talent Highlands in Guangxi Institutions of Higher Learning"partially supported by the Natural Sciences and Engineering Research Council of Canada(NSERC) through a Discovery Grant to J.Yang (Grant No.RGPIN 261283)
文摘This paper presents numerical investigation on the ore-forming fluid migration driven by tectonic deformation and thermally-induced buoyancy force in the Chanziping ore district in South China.A series of numerical scenarios are considered to examine the effect of meteoric water precipitation, the dip angle of the faults,unconformity surface,and thermal input on the ore genesis.Our computations reveal that the downward basinal fluid flow driven by extensional stress mixes with the upward basal fluid driven by the thermal input from depth at the junction of two faults at a temperature of about 200℃,triggering the precipitation of the Chanziping uranium deposit.
