The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.Howeve...The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.展开更多
Benzene derivatives are volatile organic compounds commonly present in the atmospheric environment,which are toxic and complex in composition.They have become a key regulatory object in China s atmospheric environment...Benzene derivatives are volatile organic compounds commonly present in the atmospheric environment,which are toxic and complex in composition.They have become a key regulatory object in China s atmospheric environment management.In this paper,Shimadzu Nexis GC-2030 gas chromatography was used to simultaneously detect eight benzene derivatives.According to the Environmental Monitoring—Technical Guideline on Drawing and Revising Analytical Method Standards(HJ 168-2010),the monitoring methods for benzene,toluene,ethylbenzene,p-xylene,m-xylene,isopropylbenzene,o-xylene,and styrene in the Stationary Source Emission—Determination of Benzene and Its Analogies—Bags Sampling/Direct Injection—Gas Chromatography(HJ 1261-2022)are verified,and their linear relationships,detection limits,precision and accuracy are analyzed.展开更多
Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and f...Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.展开更多
The transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical...The transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical leads.This fact further influences the measurement accuracy of thermal conductivity.To solve this problem,the influence of heat loss through the electrical leads on the heating power is studied theoretically.The mathematical formula of heat loss is deduced,and the corresponding correction model is presented.A series of measurement experiments on different materials have been conducted by using the hot disk thermal constant analyzer.The results show that the influence of the heat loss on the measurement is sensitive to different test materials and probes with different sizes.When the thermal conductivity of the material is greater than 0.2 W/(m·K),the influence of the heat loss is less than 0.16%,which can be ignored.As to the lower thermal conductivity materials,it is necessary to compensate the heat loss through the electrical leads,and the accuracy of thermal conductivity measurement can be effectively improved.展开更多
The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field ...The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.展开更多
The estimation of non-point source pollution loads into the Danjiangkou Reservoir is highly significant to environmental protection in the watershed. In order to overcome the drawbacks of traditional watershed numeric...The estimation of non-point source pollution loads into the Danjiangkou Reservoir is highly significant to environmental protection in the watershed. In order to overcome the drawbacks of traditional watershed numerical models, a base flow separation method was established coupled with a digital filtering method and a flux method. The digital filtering method has been used to separate the base flows of the Hanjiang,Tianhe, Duhe, Danjiang, Laoguan, and Qihe rivers. Based on daily discharge, base flow, and pollutant concentration data, the flux method was used to calculate the point source pollution load and non-point source pollution load. The results show that:(1) In the year 2013, the total inflow of the six rivers mentioned above accounted for 95.9% of the total inflow to the Danjiangkou Reservoir. The total pollution loads of chemical oxygen demand(CODMn) and total phosphorus(TP) from the six rivers were 58.20 103 t and 1.863 10~3 t, respectively, and the non-point source pollution loads were 39.82 10~3 t and 1.544 10~3 t, respectively, indicating that the non-point source pollution is a major factor(with a contribution rate of 68.4% for CODMnand 82.9% for TP).(2) The Hanjiang River is the most significant contributor of pollution loads to the Danjiangkou Reservoir, and its CODMnand TP contribution rates reached 79.3% and 83.2%, respectively. The Duhe River took the second place.(3) Non-point source pollution mainly occurred in the wet season in 2013, accounting for 80.8% and 90.9% of the total pollution loads of CODMnand TP, respectively. It is concluded that the emphasis of pollution control should be placed on non-point source pollution.展开更多
The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatib...The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatibility can be examined in terms of the distribution of station residuals.For an ideal distribution,the input error is held at the station where it takes place as the station residual and the error is not permitted to spread to other stations.A comparison study of two optimization methods,namely the least squares method and the absolute value method,shows that the distribution with this character constrains the input errors and minimizes their impact,which explains the much more robust performance by the absolute value method in dealing with large and isolated input errors.When the errors in the input data are systematic and/or extreme in that the basic data structure is altered by these errors,none of the optimization methods are able to function.The only means to resolve this problem is the early detection and correction of these errors through a data screening process.An efficient data screening process is of primary importance for AE/MS source location.In addition to its critical role in dealing with those systematic and extreme errors,data screening creates a favorable environment for applying optimization methods.展开更多
Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregula...Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.展开更多
The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computat...The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computational fluid dynamics(CFD)and other mathematical techniques.In this paper,the hull form optimization method based on the Rankine source method and nonlinear programming(NLP)is discussed;in the optimization process,a hull form modification function is introduced to represent an improved hull surface and to generate a new smooth hull surface by changing its frame lines and bow stem profiles under the prescribed design constraints. Numerical example is given for a practical container hull form.Finally,shape optimization of bow bulls is shown for non-protruding and protruding bow bulls.This study presents a simplified and practical design method to the select frame lines of bow bulls.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
Finite water depth effect for wave-body problems are studied by continuous Rankine source method and non-desingularized technique.Free surface and seabed surface profiles are represented by continuous panels rather th...Finite water depth effect for wave-body problems are studied by continuous Rankine source method and non-desingularized technique.Free surface and seabed surface profiles are represented by continuous panels rather than a discretization by isolated points.These panels are positioned exactly on the fluid boundary surfaces and therefore no desingularization technique is required.Space increment method is applied for both free surface source and seabed source arrangements to reduce computational cost and improve numerical efficiency.Fourth order Runge-Kutta iteration scheme is adopted on the free surface updating at every time step.The finite water depth effect is studied quantitatively for a series of cylinders with different B/T ratios.The accuracy and efficiency of the proposed model are validated by comparison with published numerical results and experimental data.Numerical results show that hydrodynamic coefficients vary for cylinder bodies with different ratios of B/T.For certain set of B/T ratios the effect of finite water depth increases quickly with the increase of motion frequency and becomes stable when frequency is relatively large.It also shows that water depths have larger hydrodynamic effects on cylinder with larger breadth to draft ratios.Both the heave added mass and damping coefficients increase across the frequency range with the water depths decrease for forced heave motion.The water depths have smaller effects on sway motion response than on heave motion response.展开更多
The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitut...The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method （VSM）, and steepest descent integration method （SDIM） were used to evaluate this type of Green＇s function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves.展开更多
Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly ...Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.展开更多
Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the...Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.展开更多
In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various accelerati...In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.展开更多
A method of source depth estimation based on the multi-path time delay difference is proposed. When the minimum time arrivals in all receiver depths are snapped to a certain time on time delay-depth plane, time delay ...A method of source depth estimation based on the multi-path time delay difference is proposed. When the minimum time arrivals in all receiver depths are snapped to a certain time on time delay-depth plane, time delay arrivals of surface-bottom reflection and bottom-surface reflection intersect at the source depth. Two hydrophones deployed vertically with a certain interval are required at least. If the receiver depths are known, the pair of time delays can be used to estimate the source depth. With the proposed method the source depth can be estimated successfully in a moderate range in the deep ocean without complicated matched-field calculations in the simulations and experiments.展开更多
A panel method is described for calculating potential flow around near-surface submarines. The method uses Havelock sources which automatically satisfy the linearized free-surface boundary condition. Outputs from the ...A panel method is described for calculating potential flow around near-surface submarines. The method uses Havelock sources which automatically satisfy the linearized free-surface boundary condition. Outputs from the method include pressure field, pressure drag, wave resistance, vertical force, trim moment and wave pattern. Comparisons are made with model tests for wave resistance of Series 58 and DARPA SUBOFF hulls, as well as with wave resistance, lift force and trim moment of three length-to-diameter variants of the DSTO Joubert submarine hull. It is found that the Havelock source panel method is capable of determining with reasonable accuracy wave resistance, vertical force and trim moment for submarine hulls. Further experimental data are required in order to assess the accuracy of the method for pressure field and wave pattern prediction. The method is implemented in the computer code“HullWave”and offers potential advantages over RANS-CFD codes in terms of speed, simplicity and robustness.展开更多
Field experiments were carried out at the National Research Center farm, Nubaria area, Behura Governorate, Egypt, to study the effect of nitrogen source, flushing schedule and irrigation method on emitter clogging. Pe...Field experiments were carried out at the National Research Center farm, Nubaria area, Behura Governorate, Egypt, to study the effect of nitrogen source, flushing schedule and irrigation method on emitter clogging. Peanut Giza 5 variety (Arachishy pogaea L.) was planted in sandy soil during two successive growing seasons (2010-2011) in the 1st week of May and harvested after 130 days. Treatments used are: 1) two irrigation methods: surface drip irrigation and sub-surface drip irrigation (SDI;SSDI), 2) nitrogen source (NS):NH4NO3, (NH4)2SO4 and Ca(NO3)2 (NS1, NS2 and NS3) and 3) flushing number (FL) 0, 1 and 4 (FL1, FL2;FL3). The experiments design was split-split plot and three replicates were used. Data obtained were subjected to statistical analysis. The main effects of treatments used on clogging per cent could be written the following ascending orders: SDI 3 2 1, NS1 2 3. Concerning the 1st interaction the following ascending orders denote their effects on clogging percent: SDI × FL3 2 1, SDI × NS1 2 3, SSDI × FL3 2 1, SSDI × NS1 2 3, FL1 × NS1 1 × NS2 1 × NS3, FL2 × NS1 2 × NS2 2 × NS3 and FL3 × NS1 3 × NS2 3 × NS3. The differences between any two treatments and/or any two interactions in clogging percent were significant at the 5% level. The effect of the 2nd interaction on clogging percent was significant at the 5% level. The maximum value of clogging (20.18%) and the lowest one (3.9%) were archived in the interactions: SSDI × FL1 × NS3 and SDI × FL3 × NS1, respectively.展开更多
Environmental uncertainty represents the limiting factor in matched-field localization. Within a Bayesian framework, both the environmental parameters, and the source parameters are considered to be unknown variables....Environmental uncertainty represents the limiting factor in matched-field localization. Within a Bayesian framework, both the environmental parameters, and the source parameters are considered to be unknown variables. However, including environmental parameters in multiple-source localization greatly increases the complexity and computational demands of the inverse problem. In the paper, the closed-form maximumlikelihood expressions for source strengths and noise variance at each frequency allow these parameters to be sampled implicitly, substantially reducing the dimensionality and difficulty of the inversion. This paper compares two Bayesian-point-estimation methods： the maximum a posteriori（MAP） approach and the marginal posterior probability density（PPD） approach to source localization. The MAP approach determines the sources locations by maximizing the PPD over all source and environmental parameters. The marginal PPD approach integrates the PPD over the unknowns to obtain a sequence of marginal probability distribution over source range or depth.Monte Carlo analysis of the two approaches for a test case involving both geoacoustic and water-column uncertainties indicates that：（1） For sensitive parameters such as source range, water depth and water sound speed, the MAP solution is better than the marginal PPD solution.（2） For the less sensitive parameters, such as,bottom sound speed, bottom density, bottom attenuation and water sound speed, when the SNR is low, the marginal PPD solution can better smooth the noise, which leads to better performance than the MAP solution.Since the source range and depth are sensitive parameters, the research shows that the MAP approach provides a slightly more reliable method to locate multiple sources in an unknown environment.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42250103 and 42174090)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(Grant No.GLAB2023ZR02)the MOST Special Fund from the State Key Laboratory of Geological Processes and Mineral Resources(Grant No.MSFGPMR2022-4).
文摘The equivalent source(ES)method in the spherical coordinate system has been widely applied to processing,reduction,field modeling,and geophysical and geological interpretation of satellite magnetic anomaly data.However,the inversion for the ES model suffers from nonuniqueness and instability,which remain unresolved.To mitigate these issues,we introduce both the minimum and flattest models into the model objective function as an alternative regularization approach in the spherical ES method.We first present the methods,then analyze the accuracy of forward calculation and test the proposed ES method in this study by using synthetic data.The experimental results from simulation data indicate that our proposed regularization effectively suppresses the Backus effect and mitigates inversion instability in the low-latitude region.Finally,we apply the proposed method to magnetic anomaly data from China Seismo-Electromagnetic Satellite-1(CSES-1)and Macao Science Satellite-1(MSS-1)magnetic measurements over Africa by constructing an ES model of the large-scale lithospheric magnetic field.Compared with existing global lithospheric magnetic field models,our ES model demonstrates good consistency at high altitudes and predicts more stable fields at low altitudes.Furthermore,we derive the reduction to the pole(RTP)magnetic anomaly fields and the apparent susceptibility contrast distribution based on the ES model.The latter correlates well with the regional tectonic framework in Africa and surroundings.
文摘Benzene derivatives are volatile organic compounds commonly present in the atmospheric environment,which are toxic and complex in composition.They have become a key regulatory object in China s atmospheric environment management.In this paper,Shimadzu Nexis GC-2030 gas chromatography was used to simultaneously detect eight benzene derivatives.According to the Environmental Monitoring—Technical Guideline on Drawing and Revising Analytical Method Standards(HJ 168-2010),the monitoring methods for benzene,toluene,ethylbenzene,p-xylene,m-xylene,isopropylbenzene,o-xylene,and styrene in the Stationary Source Emission—Determination of Benzene and Its Analogies—Bags Sampling/Direct Injection—Gas Chromatography(HJ 1261-2022)are verified,and their linear relationships,detection limits,precision and accuracy are analyzed.
基金supported by the National Natural Science Foundation of China(Nos.41564001 and 41572185)the Natural Science Foundation of Jiangxi Province(No.20151BAB203045)
文摘Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.
文摘The transient plane source(TPS)method is developed recently to measure the thermal conductivity of materials.In the measurement,the heating power is influenced by the heat which is transferred via the probe electrical leads.This fact further influences the measurement accuracy of thermal conductivity.To solve this problem,the influence of heat loss through the electrical leads on the heating power is studied theoretically.The mathematical formula of heat loss is deduced,and the corresponding correction model is presented.A series of measurement experiments on different materials have been conducted by using the hot disk thermal constant analyzer.The results show that the influence of the heat loss on the measurement is sensitive to different test materials and probes with different sizes.When the thermal conductivity of the material is greater than 0.2 W/(m·K),the influence of the heat loss is less than 0.16%,which can be ignored.As to the lower thermal conductivity materials,it is necessary to compensate the heat loss through the electrical leads,and the accuracy of thermal conductivity measurement can be effectively improved.
基金supported by the National Natural Science Foundation of China (Grant No. 50879090)
文摘The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.
基金supported by the National Key Research and Development Program of China(Grants No.2016YFC0402204 and 2016YFC0402207)
文摘The estimation of non-point source pollution loads into the Danjiangkou Reservoir is highly significant to environmental protection in the watershed. In order to overcome the drawbacks of traditional watershed numerical models, a base flow separation method was established coupled with a digital filtering method and a flux method. The digital filtering method has been used to separate the base flows of the Hanjiang,Tianhe, Duhe, Danjiang, Laoguan, and Qihe rivers. Based on daily discharge, base flow, and pollutant concentration data, the flux method was used to calculate the point source pollution load and non-point source pollution load. The results show that:(1) In the year 2013, the total inflow of the six rivers mentioned above accounted for 95.9% of the total inflow to the Danjiangkou Reservoir. The total pollution loads of chemical oxygen demand(CODMn) and total phosphorus(TP) from the six rivers were 58.20 103 t and 1.863 10~3 t, respectively, and the non-point source pollution loads were 39.82 10~3 t and 1.544 10~3 t, respectively, indicating that the non-point source pollution is a major factor(with a contribution rate of 68.4% for CODMnand 82.9% for TP).(2) The Hanjiang River is the most significant contributor of pollution loads to the Danjiangkou Reservoir, and its CODMnand TP contribution rates reached 79.3% and 83.2%, respectively. The Duhe River took the second place.(3) Non-point source pollution mainly occurred in the wet season in 2013, accounting for 80.8% and 90.9% of the total pollution loads of CODMnand TP, respectively. It is concluded that the emphasis of pollution control should be placed on non-point source pollution.
文摘The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatibility can be examined in terms of the distribution of station residuals.For an ideal distribution,the input error is held at the station where it takes place as the station residual and the error is not permitted to spread to other stations.A comparison study of two optimization methods,namely the least squares method and the absolute value method,shows that the distribution with this character constrains the input errors and minimizes their impact,which explains the much more robust performance by the absolute value method in dealing with large and isolated input errors.When the errors in the input data are systematic and/or extreme in that the basic data structure is altered by these errors,none of the optimization methods are able to function.The only means to resolve this problem is the early detection and correction of these errors through a data screening process.An efficient data screening process is of primary importance for AE/MS source location.In addition to its critical role in dealing with those systematic and extreme errors,data screening creates a favorable environment for applying optimization methods.
基金The authors wish to acknowledge financial support from the National Natural Science Foundation of China(51822407 and 51774327)Natural Science Foundation of Hunan Province in China(2018JJ1037)Innovation Driven project of Central South University(2020CX014).
文摘Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.
基金the National Natural Science Foundation of China(No.51009087)
文摘The hull form optimization concerns one of the most important applications of wave-making resistance theories.In recent years,scholars can determine the hull form by using the optimization method based on the computational fluid dynamics(CFD)and other mathematical techniques.In this paper,the hull form optimization method based on the Rankine source method and nonlinear programming(NLP)is discussed;in the optimization process,a hull form modification function is introduced to represent an improved hull surface and to generate a new smooth hull surface by changing its frame lines and bow stem profiles under the prescribed design constraints. Numerical example is given for a practical container hull form.Finally,shape optimization of bow bulls is shown for non-protruding and protruding bow bulls.This study presents a simplified and practical design method to the select frame lines of bow bulls.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
基金the support by the National Natural Science Foundation of China(No.11372184)
文摘Finite water depth effect for wave-body problems are studied by continuous Rankine source method and non-desingularized technique.Free surface and seabed surface profiles are represented by continuous panels rather than a discretization by isolated points.These panels are positioned exactly on the fluid boundary surfaces and therefore no desingularization technique is required.Space increment method is applied for both free surface source and seabed source arrangements to reduce computational cost and improve numerical efficiency.Fourth order Runge-Kutta iteration scheme is adopted on the free surface updating at every time step.The finite water depth effect is studied quantitatively for a series of cylinders with different B/T ratios.The accuracy and efficiency of the proposed model are validated by comparison with published numerical results and experimental data.Numerical results show that hydrodynamic coefficients vary for cylinder bodies with different ratios of B/T.For certain set of B/T ratios the effect of finite water depth increases quickly with the increase of motion frequency and becomes stable when frequency is relatively large.It also shows that water depths have larger hydrodynamic effects on cylinder with larger breadth to draft ratios.Both the heave added mass and damping coefficients increase across the frequency range with the water depths decrease for forced heave motion.The water depths have smaller effects on sway motion response than on heave motion response.
基金Project supported by the National Natural Science Foundation of China （No. 50879090）, and the Key Research Program of Hydrody- namics of China （No. 9140A 14030712JB 11044）
文摘The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method （VSM）, and steepest descent integration method （SDIM） were used to evaluate this type of Green＇s function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves.
文摘Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.
基金supported by National Natural Science Foundation of China (No. 41174085)Chinese Academy of Sciences (KZZD-EW-TZ-19)China Geological Survey (12120113101400)
文摘Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.
基金supported by the National Natural Science Foundation of China(Nos.11775126,11545013,11605101)the Young Elite Scientists Sponsorship Program by CAST(No.2016QNRC001)+1 种基金Science Challenge Project by MIIT of China(No.TZ2018001)Tsinghua University,Initiative Scientific Research Program。
文摘In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.
基金Supported by the National Natural Science Foundation of China under Grant No 11174235
文摘A method of source depth estimation based on the multi-path time delay difference is proposed. When the minimum time arrivals in all receiver depths are snapped to a certain time on time delay-depth plane, time delay arrivals of surface-bottom reflection and bottom-surface reflection intersect at the source depth. Two hydrophones deployed vertically with a certain interval are required at least. If the receiver depths are known, the pair of time delays can be used to estimate the source depth. With the proposed method the source depth can be estimated successfully in a moderate range in the deep ocean without complicated matched-field calculations in the simulations and experiments.
文摘A panel method is described for calculating potential flow around near-surface submarines. The method uses Havelock sources which automatically satisfy the linearized free-surface boundary condition. Outputs from the method include pressure field, pressure drag, wave resistance, vertical force, trim moment and wave pattern. Comparisons are made with model tests for wave resistance of Series 58 and DARPA SUBOFF hulls, as well as with wave resistance, lift force and trim moment of three length-to-diameter variants of the DSTO Joubert submarine hull. It is found that the Havelock source panel method is capable of determining with reasonable accuracy wave resistance, vertical force and trim moment for submarine hulls. Further experimental data are required in order to assess the accuracy of the method for pressure field and wave pattern prediction. The method is implemented in the computer code“HullWave”and offers potential advantages over RANS-CFD codes in terms of speed, simplicity and robustness.
文摘Field experiments were carried out at the National Research Center farm, Nubaria area, Behura Governorate, Egypt, to study the effect of nitrogen source, flushing schedule and irrigation method on emitter clogging. Peanut Giza 5 variety (Arachishy pogaea L.) was planted in sandy soil during two successive growing seasons (2010-2011) in the 1st week of May and harvested after 130 days. Treatments used are: 1) two irrigation methods: surface drip irrigation and sub-surface drip irrigation (SDI;SSDI), 2) nitrogen source (NS):NH4NO3, (NH4)2SO4 and Ca(NO3)2 (NS1, NS2 and NS3) and 3) flushing number (FL) 0, 1 and 4 (FL1, FL2;FL3). The experiments design was split-split plot and three replicates were used. Data obtained were subjected to statistical analysis. The main effects of treatments used on clogging per cent could be written the following ascending orders: SDI 3 2 1, NS1 2 3. Concerning the 1st interaction the following ascending orders denote their effects on clogging percent: SDI × FL3 2 1, SDI × NS1 2 3, SSDI × FL3 2 1, SSDI × NS1 2 3, FL1 × NS1 1 × NS2 1 × NS3, FL2 × NS1 2 × NS2 2 × NS3 and FL3 × NS1 3 × NS2 3 × NS3. The differences between any two treatments and/or any two interactions in clogging percent were significant at the 5% level. The effect of the 2nd interaction on clogging percent was significant at the 5% level. The maximum value of clogging (20.18%) and the lowest one (3.9%) were archived in the interactions: SSDI × FL1 × NS3 and SDI × FL3 × NS1, respectively.
基金The National Natural Science Foundation of China under contract No.11704225the Shandong Provincial Natural Science Foundation under contract No.ZR2016AQ23+1 种基金the State Key Laboratory of Acoustics of Chinese Academy of Sciences under contract No.SKLA201704the National Programe on Global Change and Air-Sea Interaction
文摘Environmental uncertainty represents the limiting factor in matched-field localization. Within a Bayesian framework, both the environmental parameters, and the source parameters are considered to be unknown variables. However, including environmental parameters in multiple-source localization greatly increases the complexity and computational demands of the inverse problem. In the paper, the closed-form maximumlikelihood expressions for source strengths and noise variance at each frequency allow these parameters to be sampled implicitly, substantially reducing the dimensionality and difficulty of the inversion. This paper compares two Bayesian-point-estimation methods： the maximum a posteriori（MAP） approach and the marginal posterior probability density（PPD） approach to source localization. The MAP approach determines the sources locations by maximizing the PPD over all source and environmental parameters. The marginal PPD approach integrates the PPD over the unknowns to obtain a sequence of marginal probability distribution over source range or depth.Monte Carlo analysis of the two approaches for a test case involving both geoacoustic and water-column uncertainties indicates that：（1） For sensitive parameters such as source range, water depth and water sound speed, the MAP solution is better than the marginal PPD solution.（2） For the less sensitive parameters, such as,bottom sound speed, bottom density, bottom attenuation and water sound speed, when the SNR is low, the marginal PPD solution can better smooth the noise, which leads to better performance than the MAP solution.Since the source range and depth are sensitive parameters, the research shows that the MAP approach provides a slightly more reliable method to locate multiple sources in an unknown environment.
