Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
A new three-dimensional laterolog array sonde(3D-LS) is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferentia...A new three-dimensional laterolog array sonde(3D-LS) is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferential detection abilities. Six investigation modes are designed using the 3D finite-element method and different investigation depths are simulated based on the pseudo-geometrical factor of the six modes. The invasion profile is described using multi-array radial logs. From the analysis of the pseudo-geometrical factor, the investigation depth of the 3D-LS is about 1.5 m for conductive invasion, which is close to that of the dual laterolog tool but greater than that of the highly integrated azimuthal laterolog sonde. The vertical and azimuthal resolution is also analyzed with the same method. The 3DLS can detect low-resistivity anomalies of 0.5 m thickness and 15? around the borehole for infinitely thick formations. This study lays the foundation for more work on 3D laterolog array sonde for evaluating low-resistivity anomalies.展开更多
The aim of the study is to simulate actual operation of an elevator. First, it designed elevator scheduling algorithm for control ing operation of the elevator;second, it simulated elevator operation by the use of obj...The aim of the study is to simulate actual operation of an elevator. First, it designed elevator scheduling algorithm for control ing operation of the elevator;second, it simulated elevator operation by the use of object-oriented programming language, in which double buffering technology was used to solve the problem of elevator pictures flicker at refreshing; final y, test correctness and rationality of the system.展开更多
The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
We discuss the feasibility of using controlled-source electromagnetic (CSEM) in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program (ODP) Leg 164 log data, we have establi...We discuss the feasibility of using controlled-source electromagnetic (CSEM) in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program (ODP) Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.展开更多
The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient e...The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
To investigate the differences and the development trends of the 400 kA aluminum reduction cell, four representative cells were deeply analyzed. By using numerical simulation methods in ANSYS software, the structure p...To investigate the differences and the development trends of the 400 kA aluminum reduction cell, four representative cells were deeply analyzed. By using numerical simulation methods in ANSYS software, the structure parameters were firstly compared, and then three-dimensional models of electric-magnetic-flow field were built and solved with finite element method(FEM). The comparison of the structures reveals that the cell bodies are similar while the current flow path and distribution ratio of bus bars are different. It appears that most of the current(70%-80%) in side A are used as the magnetic field compensation current and flow through two ends. The numerical simulation results indicate that the distributions of magnetic fields are different but all satisfy with the magnetohydrodynamics(MHD) stabilization, and the flow patterns are all two or multi vortexes with appropriate velocities. The comparison shows that all studied cells can satisfy with the physical field requirement, and the commercial applications also verify that the 400 kA cells have become the product of the mature and world's leading technology.展开更多
In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position inf...In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.展开更多
The paper deals with analysis and synthesis of non-harmonic and non-linear sources and appliances, and their interaction with harmonic power supply network. Basic idea is based on knowledge of harmonic spectrum of the...The paper deals with analysis and synthesis of non-harmonic and non-linear sources and appliances, and their interaction with harmonic power supply network. Basic idea is based on knowledge of harmonic spectrum of the sources and/or appliances, respectively. Obviously, one need to know voltage harmonic components of voltage sources (renewable with inverters,...), and current harmonic components generated by non-linear appliances (rectifiers,...). Method of investigation lies on decomposition of real electric circuit into n-harmonic separated equivalent schemes for each harmonic component. Then transient analysis will be done for each scheme separately using "impedance harmonic matrices". The important fact is that each equivalent scheme is now linearized and therefore easily calculated. Finally, the effects of each investigated schemes arc summed into resulting quantities of real non-linear electric circuit.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
Based on mirror-blocks, a totally coded algorithm (TCA) for switched-current (SI) network analysis in frequency domain is presented. The algorithm is simple, available, and suitable for any swltched-current networ...Based on mirror-blocks, a totally coded algorithm (TCA) for switched-current (SI) network analysis in frequency domain is presented. The algorithm is simple, available, and suitable for any swltched-current networks. A basis of analysis and design for switched-current networks via this algorithm is provided.展开更多
The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred s...The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred simultaneously compared with pure water.A three-dimensional finite element numerical model of multi-feed microwave heating industrial liquids continuously flowing in a meter-scale circular tube is presented.The temperature field inside the applicator tube in the cavity is solved by COMSOL Multiphysics and professional programming to describe the momentum,energy and Maxwell's equations.The evaluations of the electromagnetic field,the temperature distribution and the velocity field are simulated for the fluids dynamically heated by singleand multi-feed microwave system,respectively.Both the pilot experimental investigations and numerical results of microwave with single-feed heating for fluids with different effective permittivity and flow rates show that the presented numerical modeling makes it possible to analyze dynamic process of multi-feed microwave heating the industrial liquid.The study aids in enhancing the understanding and optimizing of dynamic process in the use of multi-feed microwave heating industrial continuous flow for a variety of material properties and technical parameters.展开更多
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
基金sponsored by the National Oil and Gas Major Projects(No.2011ZX05020-009)
文摘A new three-dimensional laterolog array sonde(3D-LS) is presented. The 3DLS is based on existing high-resolution laterolog array and azimuthal resistivity imaging sondes with radial, longitudinal, and circumferential detection abilities. Six investigation modes are designed using the 3D finite-element method and different investigation depths are simulated based on the pseudo-geometrical factor of the six modes. The invasion profile is described using multi-array radial logs. From the analysis of the pseudo-geometrical factor, the investigation depth of the 3D-LS is about 1.5 m for conductive invasion, which is close to that of the dual laterolog tool but greater than that of the highly integrated azimuthal laterolog sonde. The vertical and azimuthal resolution is also analyzed with the same method. The 3DLS can detect low-resistivity anomalies of 0.5 m thickness and 15? around the borehole for infinitely thick formations. This study lays the foundation for more work on 3D laterolog array sonde for evaluating low-resistivity anomalies.
基金Supported by the Principal Fund Project of Tarim University(TDZKSSZD201207)~~
文摘The aim of the study is to simulate actual operation of an elevator. First, it designed elevator scheduling algorithm for control ing operation of the elevator;second, it simulated elevator operation by the use of object-oriented programming language, in which double buffering technology was used to solve the problem of elevator pictures flicker at refreshing; final y, test correctness and rationality of the system.
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
基金supported by the Program for New Century Excellent Talents in University(No.NCET-04-0370)
文摘We discuss the feasibility of using controlled-source electromagnetic (CSEM) in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program (ODP) Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.
基金supported by the National Key Scientific Instrument and Equipment Development Project(No.2011YQ03013307)the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education InstitutionsKey Laboratory of Coal Resources Exploration and Comprehensive Utilization,Ministry of Land and Resources
文摘The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(51104187,51274241,61321003) supported by the National Natural Science Foundation of ChinaProject(20100162120008) supported by Doctoral Fund of Ministry of Education of China
文摘To investigate the differences and the development trends of the 400 kA aluminum reduction cell, four representative cells were deeply analyzed. By using numerical simulation methods in ANSYS software, the structure parameters were firstly compared, and then three-dimensional models of electric-magnetic-flow field were built and solved with finite element method(FEM). The comparison of the structures reveals that the cell bodies are similar while the current flow path and distribution ratio of bus bars are different. It appears that most of the current(70%-80%) in side A are used as the magnetic field compensation current and flow through two ends. The numerical simulation results indicate that the distributions of magnetic fields are different but all satisfy with the magnetohydrodynamics(MHD) stabilization, and the flow patterns are all two or multi vortexes with appropriate velocities. The comparison shows that all studied cells can satisfy with the physical field requirement, and the commercial applications also verify that the 400 kA cells have become the product of the mature and world's leading technology.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
文摘In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.
文摘The paper deals with analysis and synthesis of non-harmonic and non-linear sources and appliances, and their interaction with harmonic power supply network. Basic idea is based on knowledge of harmonic spectrum of the sources and/or appliances, respectively. Obviously, one need to know voltage harmonic components of voltage sources (renewable with inverters,...), and current harmonic components generated by non-linear appliances (rectifiers,...). Method of investigation lies on decomposition of real electric circuit into n-harmonic separated equivalent schemes for each harmonic component. Then transient analysis will be done for each scheme separately using "impedance harmonic matrices". The important fact is that each equivalent scheme is now linearized and therefore easily calculated. Finally, the effects of each investigated schemes arc summed into resulting quantities of real non-linear electric circuit.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
文摘Based on mirror-blocks, a totally coded algorithm (TCA) for switched-current (SI) network analysis in frequency domain is presented. The algorithm is simple, available, and suitable for any swltched-current networks. A basis of analysis and design for switched-current networks via this algorithm is provided.
基金Project(KKSY201503006)supported by Scientific Research Foundation of Kunming University of Science and Technology,ChinaProject(2014FD009)supported by the Applied Basic Research Foundation(Youth Program)of ChinaProject(51090385)supported by the National Natural Science Foundation of China
文摘The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred simultaneously compared with pure water.A three-dimensional finite element numerical model of multi-feed microwave heating industrial liquids continuously flowing in a meter-scale circular tube is presented.The temperature field inside the applicator tube in the cavity is solved by COMSOL Multiphysics and professional programming to describe the momentum,energy and Maxwell's equations.The evaluations of the electromagnetic field,the temperature distribution and the velocity field are simulated for the fluids dynamically heated by singleand multi-feed microwave system,respectively.Both the pilot experimental investigations and numerical results of microwave with single-feed heating for fluids with different effective permittivity and flow rates show that the presented numerical modeling makes it possible to analyze dynamic process of multi-feed microwave heating the industrial liquid.The study aids in enhancing the understanding and optimizing of dynamic process in the use of multi-feed microwave heating industrial continuous flow for a variety of material properties and technical parameters.