The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods fo...The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods for centroid measurement often necessitate custom equipment and specialized positioning devices,leading to high costs and limited accuracy.Here,we present a centroid measurement method that integrates 3D scanning technology,enabling accurate measurement of centroid across various types of objects without the need for specialized positioning fixtures.A theoretical framework for centroid measurement was established,which combined the principle of the multi-point weighing method with 3D scanning technology.The measurement accuracy was evaluated using a designed standard component.Experimental results demonstrate that the discrepancies between the theoretical and the measured centroid of a standard component with various materials and complex shapes in the X,Y,and Z directions are 0.003 mm,0.009 mm,and 0.105 mm,respectively,yielding a spatial deviation of 0.106 mm.Qualitative verification was conducted through experimental validation of three distinct types.They confirmed the reliability of the proposed method,which allowed for accurate centroid measurements of various products without requiring positioning fixtures.This advancement significantly broadened the applicability and scope of centroid measurement devices,offering new theoretical insights and methodologies for the measurement of complex parts and systems.展开更多
This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conduct...This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conductivity is proposed,given the high conductivity caused by small amounts of water distributed in a network within the low-resistivity reservoir,which mimics the eff ects of high water saturation.This approach signifi cantly improves the accuracy of hydrocarbon saturation calculations in these low-resistivity reservoirs.The corrected hydrocarbon saturation values highly matched the porosity and are consistent with experimental results.This study also establishes a discrimination process to determine whether corrections are required,verifying the eff ectiveness and accuracy of the method through an application example.展开更多
In order to improve the quality of 3D printed raspberry preserves after post-processing,microwave ovens combining infrared and microwave methods were utilized.The effects of infrared heating temperature,infrared heati...In order to improve the quality of 3D printed raspberry preserves after post-processing,microwave ovens combining infrared and microwave methods were utilized.The effects of infrared heating temperature,infrared heating time,microwave power,microwave heating time on the center temperature,moisture content,the chroma(C*),the total color difference(ΔE*),shape fidelity,hardness,and the total anthocyanin content of 3D printed raspberry preserves were analyzed by response surface method(RSM).The results showed that under combining with the two methods,infrared heating improved the fidelity and quality degradation of printed products,while microwave heating enhanced the efficiency of infrared heating.Infrared-microwave combination cooking could maintain relatively stable color appearance and shape of 3D printed raspberry preserves.The AHP–CRITIC hybrid weighting method combined with the response surface test to determine the comprehensive weights of the evaluation indicators optimized the process parameters,and the optimal process parameters were obtained:infrared heating temperature of 190℃,infrared heating time of 10 min and 30 s,microwave power of 300 W,and microwave heating time of 2 min and 6 s.The 3D printed raspberry cooking methods obtained under the optimal conditions seldom had color variation,porous structure,uniform texture,and high shape fidelity,which retained the characteristics of personalized manufacturing by 3D printing.This study could provide a reference for the postprocessing and quality control of 3D cooking methods.展开更多
The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock mater...The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(...In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(FDTD)method.The proposed structure has the hexagonal close packed of 16×20 circular rods that are suspended in the air substrate with a lattice constant of 606 nm.The plane wave expansion(PWE)method is used to obtain the band diagram for AODFF at an operating wavelength of 1550 nm.The proposed optical flip-flop achieves a low delay time of 0.2 ps and a high contrast ratio(CR)of 10.33 dB.The main advantage of this design is that the input power as low as 1 mW/μm^(2) is sufficient for its operation,since no nonlinear rods are included.In addition,the footprint of the proposed AODFF is 100μm^(2),which is smaller compared to the structures reported in the literature,and it has a fast switching frequency of 5 Tbit/s.展开更多
To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondar...To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.展开更多
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions a...A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.展开更多
Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is d...Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.展开更多
The carbonate reservoirs in the Tarim Basin are characterized by low matrix-porosity,heterogeneity and anisotropy,which make it difficult to predict and evaluate these reservoirs.The reservoir formations in Lundong ar...The carbonate reservoirs in the Tarim Basin are characterized by low matrix-porosity,heterogeneity and anisotropy,which make it difficult to predict and evaluate these reservoirs.The reservoir formations in Lundong area experienced a series of diagenesis and tectonic evolution stages.And secondary storage spaces such as fractures and dissolution caves were developed while nearly all the primary pores have disappeared.Based on a summary of different types of storage spaces and their responses in conventional logs,FMI and full waveform sonic logs which are sensitive to different reservoirs,the comprehensive probability index （CPI） method is applied to evaluating the reservoirs and a standard of reservoir classification is established.By comparing the evaluation results with actual welllogging results,the method has proven to be practical for formation evaluation of carbonate reservoirs,especially for the fractured carbonate reservoirs.In reservoir fluid identification,the multivariate stepwise discriminant analysis （MSDA） method is introduced.Combining the CPI method and MSDA method,comprehensive formation evaluation has been performed for fractured and caved carbonate reservoirs in the Tarim Basin.Additionally,on the basis of secondary pore inversion results,another new method of formation evaluation is also proposed in the discussion part of this paper.Through detailed application result analysis,the method shows a promising capability for formation evaluation of complex carbonate reservoirs dominated by various secondary pores such as holes,caves,and cracks.展开更多
The cephalopod beak is a vital hard structure with a stable configuration and has been widely used for the identification of cephalopod species. This study was conducted to determine the best standardization method fo...The cephalopod beak is a vital hard structure with a stable configuration and has been widely used for the identification of cephalopod species. This study was conducted to determine the best standardization method for identifying different species by measuring 12 morphological variables of the beaks of Illex argentinus, Ommastrephes bartramii, and Dosidicus gigas that were collected by Chinese jigging vessels. To remove the effects of size, these morphometric variables were standardized using three methods. The average ratios of the upper beak morphological variables and upper crest length of O. bartramii and D. gigas were found to be greater than those of I. argentinus. However, for lower beaks, only the average of LRL(lower rostrum length)/LCL(lower crest length), LRW(lower rostrum width)/LCL, and LLWL(lower lateral wall length)/LCL of O. bartramii and D. gigas were greater than those of I. argentinus. The ratios of beak morphological variables and crest length were found to be all significantly different among the three species(P < 0.001). Among the three standardization methods, the correct classification rate of stepwise discriminant analysis(SDA) was the highest using the ratios of beak morphological variables and crest length. Compared with hood length, the correct classification rate was slightly higher when using beak variables standardized by crest length using an allometric model. The correct classification rate of the lower beak was also found to be greater than that of the upper beak. This study indicates that the ratios of beak morphological variables to crest length could be used for interspecies and intraspecies identification. Meanwhile, the lower beak variables were found to be more effective than upper beak variables in classifying beaks found in the stomachs of predators.展开更多
An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential e...An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential equationwhich is regarded as an extended elliptic equation and whose degree Υ is expanded to the case of r>4.The efficiency ofthe method is demonstrated by the KdV equation and the variant Boussinesq equations.The results indicate that themethod not only offers all solutions obtained by using Fu's and Fan's methods,but also some new solutions.展开更多
For the gray attributes of the equipment program and its difficulty to carry out the quantitative assessment of the equipment program information, the gray relation projection method is simply reviewed. Combining the ...For the gray attributes of the equipment program and its difficulty to carry out the quantitative assessment of the equipment program information, the gray relation projection method is simply reviewed. Combining the super-data envelopment analysis（DEA） model and the gray system theory, a new super-DEA for measuring the weight is proposed, and a gray relation projection model is established to rank the equipment programs. Finally, this approach is used to evaluate the equipment program. The results are verified valid and can provide a new way for evaluating the equipment program.展开更多
Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured...Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.展开更多
The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical ...The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.展开更多
基金supported by National Natural Science Foundation of China(No.52176122).
文摘The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods for centroid measurement often necessitate custom equipment and specialized positioning devices,leading to high costs and limited accuracy.Here,we present a centroid measurement method that integrates 3D scanning technology,enabling accurate measurement of centroid across various types of objects without the need for specialized positioning fixtures.A theoretical framework for centroid measurement was established,which combined the principle of the multi-point weighing method with 3D scanning technology.The measurement accuracy was evaluated using a designed standard component.Experimental results demonstrate that the discrepancies between the theoretical and the measured centroid of a standard component with various materials and complex shapes in the X,Y,and Z directions are 0.003 mm,0.009 mm,and 0.105 mm,respectively,yielding a spatial deviation of 0.106 mm.Qualitative verification was conducted through experimental validation of three distinct types.They confirmed the reliability of the proposed method,which allowed for accurate centroid measurements of various products without requiring positioning fixtures.This advancement significantly broadened the applicability and scope of centroid measurement devices,offering new theoretical insights and methodologies for the measurement of complex parts and systems.
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
文摘This study discusses the challenges in logging the evaluation of low-resistivity oil reservoirs,especially the difficult problems involving their saturation calculation.A correction method for equivalent water conductivity is proposed,given the high conductivity caused by small amounts of water distributed in a network within the low-resistivity reservoir,which mimics the eff ects of high water saturation.This approach signifi cantly improves the accuracy of hydrocarbon saturation calculations in these low-resistivity reservoirs.The corrected hydrocarbon saturation values highly matched the porosity and are consistent with experimental results.This study also establishes a discrimination process to determine whether corrections are required,verifying the eff ectiveness and accuracy of the method through an application example.
基金Supported by the National Natural Science Foundation of China(32072352)。
文摘In order to improve the quality of 3D printed raspberry preserves after post-processing,microwave ovens combining infrared and microwave methods were utilized.The effects of infrared heating temperature,infrared heating time,microwave power,microwave heating time on the center temperature,moisture content,the chroma(C*),the total color difference(ΔE*),shape fidelity,hardness,and the total anthocyanin content of 3D printed raspberry preserves were analyzed by response surface method(RSM).The results showed that under combining with the two methods,infrared heating improved the fidelity and quality degradation of printed products,while microwave heating enhanced the efficiency of infrared heating.Infrared-microwave combination cooking could maintain relatively stable color appearance and shape of 3D printed raspberry preserves.The AHP–CRITIC hybrid weighting method combined with the response surface test to determine the comprehensive weights of the evaluation indicators optimized the process parameters,and the optimal process parameters were obtained:infrared heating temperature of 190℃,infrared heating time of 10 min and 30 s,microwave power of 300 W,and microwave heating time of 2 min and 6 s.The 3D printed raspberry cooking methods obtained under the optimal conditions seldom had color variation,porous structure,uniform texture,and high shape fidelity,which retained the characteristics of personalized manufacturing by 3D printing.This study could provide a reference for the postprocessing and quality control of 3D cooking methods.
基金supported by the National Natural Science Foundation of China(Grant Nos.42172312 and 52211540395)support from the Institut Universitaire de France(IUF).
文摘The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.
基金supported by the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
文摘In this study,we design and numerically investigate a novel all optical D flip-flop(AODFF)based on linear photonic crystal(LPhC)structure that is composed of optical waveguides using the finite difference time domain(FDTD)method.The proposed structure has the hexagonal close packed of 16×20 circular rods that are suspended in the air substrate with a lattice constant of 606 nm.The plane wave expansion(PWE)method is used to obtain the band diagram for AODFF at an operating wavelength of 1550 nm.The proposed optical flip-flop achieves a low delay time of 0.2 ps and a high contrast ratio(CR)of 10.33 dB.The main advantage of this design is that the input power as low as 1 mW/μm^(2) is sufficient for its operation,since no nonlinear rods are included.In addition,the footprint of the proposed AODFF is 100μm^(2),which is smaller compared to the structures reported in the literature,and it has a fast switching frequency of 5 Tbit/s.
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
文摘3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
基金Project (50975263) supported by the National Natural Science Foundation of ChinaProject (2010081015) supported by International Cooperation Project of Shanxi Province, China+1 种基金 Project (2010-78) supported by the Scholarship Council in Shanxi province, ChinaProject (2010420120005) supported by Doctoral Fund of Ministry of Education of China
文摘A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.
基金supported by the National Scientific and Technological Plan(Nos.2009BAB43B00 and 2009BAB43B01)
文摘Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘The carbonate reservoirs in the Tarim Basin are characterized by low matrix-porosity,heterogeneity and anisotropy,which make it difficult to predict and evaluate these reservoirs.The reservoir formations in Lundong area experienced a series of diagenesis and tectonic evolution stages.And secondary storage spaces such as fractures and dissolution caves were developed while nearly all the primary pores have disappeared.Based on a summary of different types of storage spaces and their responses in conventional logs,FMI and full waveform sonic logs which are sensitive to different reservoirs,the comprehensive probability index （CPI） method is applied to evaluating the reservoirs and a standard of reservoir classification is established.By comparing the evaluation results with actual welllogging results,the method has proven to be practical for formation evaluation of carbonate reservoirs,especially for the fractured carbonate reservoirs.In reservoir fluid identification,the multivariate stepwise discriminant analysis （MSDA） method is introduced.Combining the CPI method and MSDA method,comprehensive formation evaluation has been performed for fractured and caved carbonate reservoirs in the Tarim Basin.Additionally,on the basis of secondary pore inversion results,another new method of formation evaluation is also proposed in the discussion part of this paper.Through detailed application result analysis,the method shows a promising capability for formation evaluation of complex carbonate reservoirs dominated by various secondary pores such as holes,caves,and cracks.
基金supported by the National Natural Science Foundation of China(Nos.41306127 and 41276156)the National Science Foundation of Shanghai(No.13ZR1419700)+3 种基金the Innovation Program of Shanghai Municipal Education Commission(No.13YZ091)the Shanghai Leading Academic Discipline Project(Fisheries Discipline)supported by Shanghai Ocean University(SHOU)International Center for Marine StudiesShanghai Visiting 1000 Talent Program
文摘The cephalopod beak is a vital hard structure with a stable configuration and has been widely used for the identification of cephalopod species. This study was conducted to determine the best standardization method for identifying different species by measuring 12 morphological variables of the beaks of Illex argentinus, Ommastrephes bartramii, and Dosidicus gigas that were collected by Chinese jigging vessels. To remove the effects of size, these morphometric variables were standardized using three methods. The average ratios of the upper beak morphological variables and upper crest length of O. bartramii and D. gigas were found to be greater than those of I. argentinus. However, for lower beaks, only the average of LRL(lower rostrum length)/LCL(lower crest length), LRW(lower rostrum width)/LCL, and LLWL(lower lateral wall length)/LCL of O. bartramii and D. gigas were greater than those of I. argentinus. The ratios of beak morphological variables and crest length were found to be all significantly different among the three species(P < 0.001). Among the three standardization methods, the correct classification rate of stepwise discriminant analysis(SDA) was the highest using the ratios of beak morphological variables and crest length. Compared with hood length, the correct classification rate was slightly higher when using beak variables standardized by crest length using an allometric model. The correct classification rate of the lower beak was also found to be greater than that of the upper beak. This study indicates that the ratios of beak morphological variables to crest length could be used for interspecies and intraspecies identification. Meanwhile, the lower beak variables were found to be more effective than upper beak variables in classifying beaks found in the stomachs of predators.
基金National Natural Science Foundation of China under Grant No.10672053
文摘An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential equationwhich is regarded as an extended elliptic equation and whose degree Υ is expanded to the case of r>4.The efficiency ofthe method is demonstrated by the KdV equation and the variant Boussinesq equations.The results indicate that themethod not only offers all solutions obtained by using Fu's and Fan's methods,but also some new solutions.
基金supported by the National Natural Science Foundation of China(7107307971222106+2 种基金70901069)the Research Foundation of the National Excellent Doctoral Dissertation of Chinathe Research Fund for the Doctoral Program of Higher Education(20133402110028)
文摘For the gray attributes of the equipment program and its difficulty to carry out the quantitative assessment of the equipment program information, the gray relation projection method is simply reviewed. Combining the super-data envelopment analysis（DEA） model and the gray system theory, a new super-DEA for measuring the weight is proposed, and a gray relation projection model is established to rank the equipment programs. Finally, this approach is used to evaluate the equipment program. The results are verified valid and can provide a new way for evaluating the equipment program.
基金financially supported by the National Natural Science Foundation of China(No.41774125)Key Program of National Natural Science Foundation of China(No.41530320)+1 种基金the Key National Research Project of China(Nos.2016YFC0303100 and 2017YFC0601900)the Strategic Priority Research Program of Chinese Academy of Sciences Pilot Special(No.XDA 14020102)
文摘Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.
基金National Natural Science Foundation of China under Grant No.51078337,No.51108431 and No.91315301
文摘The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.
