In the middle and second half of oilfield development,further exploration and development of petroleum reservoirs in complex fault blocks are major tasks.Based on the characteristics of fully developed farewell faults...In the middle and second half of oilfield development,further exploration and development of petroleum reservoirs in complex fault blocks are major tasks.Based on the characteristics of fully developed farewell faults in the Wenweigu oil field and the data from 3D earthquakes,RFT,HDT and other techniques,farewell faults and their distribution were identified.Conformation of the fault blocks has been provided and the precision of describing the farewell faults improved.Research technology of farewell faults has applied these methods in this region.展开更多
In industrial processes,there exist faults that have complex effect on process variables.Complex and simple faults are defined according to their effect dimensions.The conventional approaches based on structured resid...In industrial processes,there exist faults that have complex effect on process variables.Complex and simple faults are defined according to their effect dimensions.The conventional approaches based on structured residuals cannot isolate complex faults.This paper presents a multi-level strategy for complex fault isolation.An extraction procedure is employed to reduce the complex faults to simple ones and assign them to several levels.On each level,faults are isolated by their different responses in the structured residuals.Each residual is obtained insensitive to one fault but more sensitive to others.The faults on different levels are verified to have different residual responses and will not be confused.An entire incidence matrix containing residual response characteristics of all faults is obtained,based on which faults can be isolated.The proposed method is applied in the Tennessee Eastman process example,and the effectiveness and advantage are demonstrated.展开更多
The spatial and temproal evolution of strain. fault displacement and acoustic emissions during deformation of fault systems with different geometrical textures are studied experimentally under biaxial compresison, and...The spatial and temproal evolution of strain. fault displacement and acoustic emissions during deformation of fault systems with different geometrical textures are studied experimentally under biaxial compresison, and the characteristics of typical instability events are analysed. The results show that fault systems with different geometrical textures have different evolutional images of physical field during deformation. Based on the characteristics of physical field and the deformation mechanism, various types of instability - two types of stick-slip, fracturing type and mixed type instability can be recognized. Different types of instability differ clearly in their precursors, and the instability type is closely related with the geometrical texture and the deformation stage of the fault system. Therefore, it is very significant for earthquake prediction and precursor analysis to investigatethe geometrical textures of natural active faults.展开更多
The geological background of the Subei basin is that of small relief subsidence, low geothermal gradient, multi-sedimentary hiatuses, intense reconstruction of the basin, frequent magmatic activity, and a unique combi...The geological background of the Subei basin is that of small relief subsidence, low geothermal gradient, multi-sedimentary hiatuses, intense reconstruction of the basin, frequent magmatic activity, and a unique combination of source rock and reservoir. This geological background resulted in characteristics such as many small fault blocks, multiple oil-bearing formations, scattered oil distribution, mini- to small-sized reservoirs, and difficulties in exploration. Aimed at such characteristics, an effective exploration strategy was adopted, and the oil reserves, production and economic benefits of the Jiangsu oilfield were significantly increased. This exploration strategy included understanding the hydrocarbon generation mechanism of source rocks, progressive evaluation of oil resources, comprehensive research on the faulted systems, the distribution of oil reservoirs and their controlling factors. The techniques used included integration of acquisition, processing and interpretation with 3-D seismics as the core technology, trap description and evaluation, directional drilling and cluster drilling, integration of cuttings logging, gas chromatographic logging and geochemical logging, and integration of early reservoir description and progressive exploration and development. This strategy could be guidance for other complex fault blocks.展开更多
Because the extract of the weak failure information is always the difficulty and focus of fault detection. Aiming for specific statistical properties of complex wavelet coefficients of gearbox vibration signals, a new...Because the extract of the weak failure information is always the difficulty and focus of fault detection. Aiming for specific statistical properties of complex wavelet coefficients of gearbox vibration signals, a new signal-denoising method which uses local adaptive algorithm based on dual-tree complex wavelet transform (DT-CWT) is introduced to extract weak failure information in gear, especially to extract impulse components. By taking into account the non-Gaussian probability distribution and the statistical dependencies among wavelet coefficients of some signals, and by taking the advantage of near shift-invariance of DT-CWT, the higher signal-to-noise ratio (SNR) than common wavelet denoising methods can be obtained. Experiments of extracting periodic impulses in gearbox vibration signals indicate that the method can extract incipient fault feature and hidden information from heavy noise, and it has an excellent effect on identifying weak feature signals in gearbox vibration signals.展开更多
Take the lacustrine delta in the second member of Paleogene Shahejie Formation in block Wang43, Dongying depression, Bohai Bay Basin as an example, the deposition architectural characteristics of lacustrine deltaic sa...Take the lacustrine delta in the second member of Paleogene Shahejie Formation in block Wang43, Dongying depression, Bohai Bay Basin as an example, the deposition architectural characteristics of lacustrine deltaic sand bodies controlled by syn-depositional faults in complex fault blocks of rift basin are examined from the aspect of the tectonic-deposition response, using cores, well logs and three-dimensional seismic data. The small-scale syn-depositional faults in complex fault blocks are dense and different in dip, the activity along the strike of syn-depositional fault varies in different positions, and all these control the sedimentary process of deltaic sand bodies. Influenced by syn-depositional faults, the deltaic distributary channel is more likely to pass through the position with weak fault activity, and be deflected or limited at the position with strong fault activity. In downthrown side of a single syn-depositional fault or micro-graben areas, sand bodies increase in thickness and planar scale, and sand bodies of multiple stages are likely to stack over each other vertically. In micro-horst areas controlled by syn-depositional faults, the sand bodies decrease in abundance, and appear in intermittent superimposed pattern vertically. This study can provide new research ideas and theoretical basis for exploration and development research in complex fault blocks.展开更多
The boundary integral equation method (BIEM) is now widely used in numerical studies on earthquake rupture dynamics, and is proved to be a powerful tool to deal with problems on complex fault system. However, since ...The boundary integral equation method (BIEM) is now widely used in numerical studies on earthquake rupture dynamics, and is proved to be a powerful tool to deal with problems on complex fault system. However, since this method heavily lies on the specific forms of Green's function and only the Green's function in full-space has a closed analytic expression, it is usually limited to a full-space medium. In this study, as a first step to extend this method to an arbitrary complex fault system in half-space, the boundary integral equations (BIEs) for dynamic strike-slip on vertical complex fault system in half-space are derived based on exact Green's function for isotropic and homogeneous half-space. Effect of the geometry of the complex fault system are dealt with carefully. Final BIEs is composed of two parts: contribution from full-space, which has been thoroughly investigated by Aochi and his co-workers by using the Green's function for full-space, and that from free surface, which is studied in detail in this study.展开更多
The first-principles method based on the projector augmented wave method within the generalized gradient approximation was employed to calculate the superlattice intrinsic stacking fault(SISF) and complex stacking f...The first-principles method based on the projector augmented wave method within the generalized gradient approximation was employed to calculate the superlattice intrinsic stacking fault(SISF) and complex stacking fault(CSF) energies of the binary Ni3Al alloys with different Al contents and the ternary Ni3Al intermetallic alloys with addition of alloying elements,such as Pd,Pt,Ti,Mo,Ta,W and Re.The results show that the energies of SISF and CSF increase significantly with increase of Al contents in Ni3Al.Addition of Pd and Pt occupying the Ni sublattices does not change the SISF and CSF energies of Ni3Al markedly in comparison with the Ni-23.75Al alloy.While addition of alloying elements,such as Ti,Mo,Ta,W and Re,occupying the Al sublattices dramatically increases the SISF and CSF energies of Ni3Al.The results suggest that the energies of SISF and CSF are dependent both on the Al contents and on the site occupancy of the ternary alloying element in Ni3Al intermetallic alloys.展开更多
Analysis of the in situ stress orientation and magnitude in the No.4 Structure of Nanpu Sag was performed on the basis of data obtained from borehole breakout and acoustic emission measurements.On the basis of mechani...Analysis of the in situ stress orientation and magnitude in the No.4 Structure of Nanpu Sag was performed on the basis of data obtained from borehole breakout and acoustic emission measurements.On the basis of mechanical experiments,logging interpretation,and seismic data,a 3 D geological model and heterogeneous rock mechanics field of the reservoir were constructed.Finite element simulation techniques were then used for the detailed prediction of the 3 D stress field.The results indicated that the maximum horizontal stress orientation in the study area was generally NEE-SWW trending,with significant changes in the in situ stress orientation within and between fault blocks.Along surfaces and profiles,stress magnitudes were discrete and the in situ stress belonged to theⅠa-type.Observed inter-strata differences were characterized as five different types of in situ stress profile.Faults were the most important factor causing large distributional differences in the stress field of reservoirs within the complex fault blocks.The next important influence on the stress field was the reservoir’s rock mechanics parameters,which impacted on the magnitudes of in situ stress magnitudes.This technique provided a theoretical basis for more efficient exploration and development of low-permeability reservoirs within complex fault blocks.展开更多
This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our alg...This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our algorithm has very high computation precision and efficiency; therefore, it is a very practical tool to investigate strong ground motion problems. Additionally, the equations given in this study can also be used to invert the fault rupture process.展开更多
The Queshan MCC is an important example of a crustal extensional structure in the eastern Jiaodong Peninsula along the southeastern margin of the NCC in the Early Cretaceous. The MCC is a typical Cordilleran-type core...The Queshan MCC is an important example of a crustal extensional structure in the eastern Jiaodong Peninsula along the southeastern margin of the NCC in the Early Cretaceous. The MCC is a typical Cordilleran-type core complex with a three-layered structure:(1) the upper plate is constituted by the Cretaceous supradetachment basin and Paleoproterozoic basement;(2) the lower plate comprises the Neoarchean high-grade metamorphic complexes and late Mesozoic granitic intrusions; and(3) the two plates are separated by a master detachment fault. A series of late NEN-oriented brittle faults superimposed on and destructed the early MCC. Petrology, geometry, kinematics, macro- and micro-structures and quartz c-axis fabrics imply that the MCC has a progressive exhumation history from middle-lower to subsurface level(via middle-upper crustal level) under the nearly WNW-ESE regional extensional regime. We present structural and geochronological evidence to constrain the exhumation of the Queshan MCC from ca. 135 to 113 Ma. Based on the comprehensive analysis of the different patterns of extensional structures in the Jiaodong and Liaodong Peninsula, we have defined the Jiao-Liao Early Cretaceou extensional province and further divided the crustal extension of it into two stages: the first stage was the intense flow of the middle-lower crust and the second stage was the extension of the middle-upper crust. Combining the tectonic setting, the lithosphere thinning in the Jiao-Liao Early Cretaceous extensional province can be considered a typical model for the response of crust-mantle detachment faulting under regional extension in East Asia.展开更多
Fluid flow in hydrocarbon reservoirs and consequently optimum scenario for hydrocarbon production,is heavily influenced by reservoir heterogeneities.Faults are one of the most common types of heterogeneity found in re...Fluid flow in hydrocarbon reservoirs and consequently optimum scenario for hydrocarbon production,is heavily influenced by reservoir heterogeneities.Faults are one of the most common types of heterogeneity found in reservoirs.Leaky faults,baffles(limited extent faults)and complex multiple fault geometries are among the most complicated and important types of faults that are difficult to characterize.Leaky faults,unlike the sealing faults,are in partial communication with other portions of the reservoir.Because of faults'effect on reservoir connectivity and possible infill drilling plan for accessing all parts of the reservoirs,possible communication across the fault must be precisely modeled.In order to detect the effect of a fault on communication within the reservoir,we need to analyze dynamic data.There are a few analytical methods for modelling partially communicating faults,however,these methods may not be accurate enough and may be limited in application,especially in complex situations.Numerical methods(i.e.finite difference or finite element)are also not computationally economical when a large number of grid blocks are simulated.In the current work,the Fast Marching Method(FMM)is applied to effectively mimic fluid flow in the heterogeneous areas,such as complex faults.It is shown that FMM can capture the effect of different fault configurations on the bottom hole pressure and is also able to capture different linear,radial and spherical flows.展开更多
Molecular dynamics simulation is carried out to study the nanometric machining of single crystal Nickel(Ni). Through an investigation of atomic displacement and the variation of cutting force, it is found that the lat...Molecular dynamics simulation is carried out to study the nanometric machining of single crystal Nickel(Ni). Through an investigation of atomic displacement and the variation of cutting force, it is found that the latter is in accordance with the number variation of elastic displaced atoms in the workpiece. It is further found that the generation of complex stacking faults is the predominant cause of cutting force fluctuation, and the stacking faults with complex structures lead to work-hardening. The temperature of the cutting tool and workpiece is studied during the machining process. It is concluded that the selection of averaging steps has a significant influence on the system temperature distribution. Thus, the time-spatial averaging method, which has a high accuracy and consistency in temperature distribution, is proposed.展开更多
文摘In the middle and second half of oilfield development,further exploration and development of petroleum reservoirs in complex fault blocks are major tasks.Based on the characteristics of fully developed farewell faults in the Wenweigu oil field and the data from 3D earthquakes,RFT,HDT and other techniques,farewell faults and their distribution were identified.Conformation of the fault blocks has been provided and the precision of describing the farewell faults improved.Research technology of farewell faults has applied these methods in this region.
基金Supported by the National Natural Science Foundation of China(60574047)the National High Technology Research and Development Program of China(2007AA04Z168,2009AA04Z154)the Research Fund for the Doctoral Program of Higher Education in China(20050335018)
文摘In industrial processes,there exist faults that have complex effect on process variables.Complex and simple faults are defined according to their effect dimensions.The conventional approaches based on structured residuals cannot isolate complex faults.This paper presents a multi-level strategy for complex fault isolation.An extraction procedure is employed to reduce the complex faults to simple ones and assign them to several levels.On each level,faults are isolated by their different responses in the structured residuals.Each residual is obtained insensitive to one fault but more sensitive to others.The faults on different levels are verified to have different residual responses and will not be confused.An entire incidence matrix containing residual response characteristics of all faults is obtained,based on which faults can be isolated.The proposed method is applied in the Tennessee Eastman process example,and the effectiveness and advantage are demonstrated.
文摘The spatial and temproal evolution of strain. fault displacement and acoustic emissions during deformation of fault systems with different geometrical textures are studied experimentally under biaxial compresison, and the characteristics of typical instability events are analysed. The results show that fault systems with different geometrical textures have different evolutional images of physical field during deformation. Based on the characteristics of physical field and the deformation mechanism, various types of instability - two types of stick-slip, fracturing type and mixed type instability can be recognized. Different types of instability differ clearly in their precursors, and the instability type is closely related with the geometrical texture and the deformation stage of the fault system. Therefore, it is very significant for earthquake prediction and precursor analysis to investigatethe geometrical textures of natural active faults.
文摘The geological background of the Subei basin is that of small relief subsidence, low geothermal gradient, multi-sedimentary hiatuses, intense reconstruction of the basin, frequent magmatic activity, and a unique combination of source rock and reservoir. This geological background resulted in characteristics such as many small fault blocks, multiple oil-bearing formations, scattered oil distribution, mini- to small-sized reservoirs, and difficulties in exploration. Aimed at such characteristics, an effective exploration strategy was adopted, and the oil reserves, production and economic benefits of the Jiangsu oilfield were significantly increased. This exploration strategy included understanding the hydrocarbon generation mechanism of source rocks, progressive evaluation of oil resources, comprehensive research on the faulted systems, the distribution of oil reservoirs and their controlling factors. The techniques used included integration of acquisition, processing and interpretation with 3-D seismics as the core technology, trap description and evaluation, directional drilling and cluster drilling, integration of cuttings logging, gas chromatographic logging and geochemical logging, and integration of early reservoir description and progressive exploration and development. This strategy could be guidance for other complex fault blocks.
基金Beijing Municipal Natural Science Foundation of China (No. 3062012).
文摘Because the extract of the weak failure information is always the difficulty and focus of fault detection. Aiming for specific statistical properties of complex wavelet coefficients of gearbox vibration signals, a new signal-denoising method which uses local adaptive algorithm based on dual-tree complex wavelet transform (DT-CWT) is introduced to extract weak failure information in gear, especially to extract impulse components. By taking into account the non-Gaussian probability distribution and the statistical dependencies among wavelet coefficients of some signals, and by taking the advantage of near shift-invariance of DT-CWT, the higher signal-to-noise ratio (SNR) than common wavelet denoising methods can be obtained. Experiments of extracting periodic impulses in gearbox vibration signals indicate that the method can extract incipient fault feature and hidden information from heavy noise, and it has an excellent effect on identifying weak feature signals in gearbox vibration signals.
基金Supported by the China National Science and Technology Major Project(2016ZX05011-002,2016ZX05010-001,2016ZX05011-001)the National Basic Research Program(973)Program(2015CB250901)the National Natural Science Foundation of China(41902122)
文摘Take the lacustrine delta in the second member of Paleogene Shahejie Formation in block Wang43, Dongying depression, Bohai Bay Basin as an example, the deposition architectural characteristics of lacustrine deltaic sand bodies controlled by syn-depositional faults in complex fault blocks of rift basin are examined from the aspect of the tectonic-deposition response, using cores, well logs and three-dimensional seismic data. The small-scale syn-depositional faults in complex fault blocks are dense and different in dip, the activity along the strike of syn-depositional fault varies in different positions, and all these control the sedimentary process of deltaic sand bodies. Influenced by syn-depositional faults, the deltaic distributary channel is more likely to pass through the position with weak fault activity, and be deflected or limited at the position with strong fault activity. In downthrown side of a single syn-depositional fault or micro-graben areas, sand bodies increase in thickness and planar scale, and sand bodies of multiple stages are likely to stack over each other vertically. In micro-horst areas controlled by syn-depositional faults, the sand bodies decrease in abundance, and appear in intermittent superimposed pattern vertically. This study can provide new research ideas and theoretical basis for exploration and development research in complex fault blocks.
基金supported by the President Fund of GUCAS(No. O85101CM03)National Natural Science Foundation of China(Nos.90715019 and 40821062)partially by National Basic Research Program of China (No.2004CB418404)
文摘The boundary integral equation method (BIEM) is now widely used in numerical studies on earthquake rupture dynamics, and is proved to be a powerful tool to deal with problems on complex fault system. However, since this method heavily lies on the specific forms of Green's function and only the Green's function in full-space has a closed analytic expression, it is usually limited to a full-space medium. In this study, as a first step to extend this method to an arbitrary complex fault system in half-space, the boundary integral equations (BIEs) for dynamic strike-slip on vertical complex fault system in half-space are derived based on exact Green's function for isotropic and homogeneous half-space. Effect of the geometry of the complex fault system are dealt with carefully. Final BIEs is composed of two parts: contribution from full-space, which has been thoroughly investigated by Aochi and his co-workers by using the Green's function for full-space, and that from free surface, which is studied in detail in this study.
基金Project(50871065) supported by the National Natural Science Foundation of ChinaProjects(08DJ1400402,09JC1407200,10DZ2290904) supported by the Science and Technology Committee of Shanghai Municipality,China
文摘The first-principles method based on the projector augmented wave method within the generalized gradient approximation was employed to calculate the superlattice intrinsic stacking fault(SISF) and complex stacking fault(CSF) energies of the binary Ni3Al alloys with different Al contents and the ternary Ni3Al intermetallic alloys with addition of alloying elements,such as Pd,Pt,Ti,Mo,Ta,W and Re.The results show that the energies of SISF and CSF increase significantly with increase of Al contents in Ni3Al.Addition of Pd and Pt occupying the Ni sublattices does not change the SISF and CSF energies of Ni3Al markedly in comparison with the Ni-23.75Al alloy.While addition of alloying elements,such as Ti,Mo,Ta,W and Re,occupying the Al sublattices dramatically increases the SISF and CSF energies of Ni3Al.The results suggest that the energies of SISF and CSF are dependent both on the Al contents and on the site occupancy of the ternary alloying element in Ni3Al intermetallic alloys.
基金financially supported by the National Oil and Gas Major Project(2016ZX05047-003,2016ZX05014002-006)the National Natural Science Foundation of China(41572124)the Fundamental Research Funds for the Central Universities(17CX05010)
文摘Analysis of the in situ stress orientation and magnitude in the No.4 Structure of Nanpu Sag was performed on the basis of data obtained from borehole breakout and acoustic emission measurements.On the basis of mechanical experiments,logging interpretation,and seismic data,a 3 D geological model and heterogeneous rock mechanics field of the reservoir were constructed.Finite element simulation techniques were then used for the detailed prediction of the 3 D stress field.The results indicated that the maximum horizontal stress orientation in the study area was generally NEE-SWW trending,with significant changes in the in situ stress orientation within and between fault blocks.Along surfaces and profiles,stress magnitudes were discrete and the in situ stress belonged to theⅠa-type.Observed inter-strata differences were characterized as five different types of in situ stress profile.Faults were the most important factor causing large distributional differences in the stress field of reservoirs within the complex fault blocks.The next important influence on the stress field was the reservoir’s rock mechanics parameters,which impacted on the magnitudes of in situ stress magnitudes.This technique provided a theoretical basis for more efficient exploration and development of low-permeability reservoirs within complex fault blocks.
基金National Natural Science Foundation of China (40474011 and 40521002).
文摘This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our algorithm has very high computation precision and efficiency; therefore, it is a very practical tool to investigate strong ground motion problems. Additionally, the equations given in this study can also be used to invert the fault rupture process.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41430211, 90814006 & 91214301)Doctoral Foundation of Ministry of Education of China (Grant No. 20110022130001)
文摘The Queshan MCC is an important example of a crustal extensional structure in the eastern Jiaodong Peninsula along the southeastern margin of the NCC in the Early Cretaceous. The MCC is a typical Cordilleran-type core complex with a three-layered structure:(1) the upper plate is constituted by the Cretaceous supradetachment basin and Paleoproterozoic basement;(2) the lower plate comprises the Neoarchean high-grade metamorphic complexes and late Mesozoic granitic intrusions; and(3) the two plates are separated by a master detachment fault. A series of late NEN-oriented brittle faults superimposed on and destructed the early MCC. Petrology, geometry, kinematics, macro- and micro-structures and quartz c-axis fabrics imply that the MCC has a progressive exhumation history from middle-lower to subsurface level(via middle-upper crustal level) under the nearly WNW-ESE regional extensional regime. We present structural and geochronological evidence to constrain the exhumation of the Queshan MCC from ca. 135 to 113 Ma. Based on the comprehensive analysis of the different patterns of extensional structures in the Jiaodong and Liaodong Peninsula, we have defined the Jiao-Liao Early Cretaceou extensional province and further divided the crustal extension of it into two stages: the first stage was the intense flow of the middle-lower crust and the second stage was the extension of the middle-upper crust. Combining the tectonic setting, the lithosphere thinning in the Jiao-Liao Early Cretaceous extensional province can be considered a typical model for the response of crust-mantle detachment faulting under regional extension in East Asia.
文摘Fluid flow in hydrocarbon reservoirs and consequently optimum scenario for hydrocarbon production,is heavily influenced by reservoir heterogeneities.Faults are one of the most common types of heterogeneity found in reservoirs.Leaky faults,baffles(limited extent faults)and complex multiple fault geometries are among the most complicated and important types of faults that are difficult to characterize.Leaky faults,unlike the sealing faults,are in partial communication with other portions of the reservoir.Because of faults'effect on reservoir connectivity and possible infill drilling plan for accessing all parts of the reservoirs,possible communication across the fault must be precisely modeled.In order to detect the effect of a fault on communication within the reservoir,we need to analyze dynamic data.There are a few analytical methods for modelling partially communicating faults,however,these methods may not be accurate enough and may be limited in application,especially in complex situations.Numerical methods(i.e.finite difference or finite element)are also not computationally economical when a large number of grid blocks are simulated.In the current work,the Fast Marching Method(FMM)is applied to effectively mimic fluid flow in the heterogeneous areas,such as complex faults.It is shown that FMM can capture the effect of different fault configurations on the bottom hole pressure and is also able to capture different linear,radial and spherical flows.
基金supported by the National Natural Science Foundation of China(Grant No.51375082)
文摘Molecular dynamics simulation is carried out to study the nanometric machining of single crystal Nickel(Ni). Through an investigation of atomic displacement and the variation of cutting force, it is found that the latter is in accordance with the number variation of elastic displaced atoms in the workpiece. It is further found that the generation of complex stacking faults is the predominant cause of cutting force fluctuation, and the stacking faults with complex structures lead to work-hardening. The temperature of the cutting tool and workpiece is studied during the machining process. It is concluded that the selection of averaging steps has a significant influence on the system temperature distribution. Thus, the time-spatial averaging method, which has a high accuracy and consistency in temperature distribution, is proposed.