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.展开更多
We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topograp...We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.展开更多
Using the software ANSYS-19.2/Explicit Dynamics,this study performedfinite-element modeling of the large-diameter steel pipeline cross-section for the Beineu-Bozoy-Shymkent gas pipeline with a non-through straight crac...Using the software ANSYS-19.2/Explicit Dynamics,this study performedfinite-element modeling of the large-diameter steel pipeline cross-section for the Beineu-Bozoy-Shymkent gas pipeline with a non-through straight crack,strengthened by steel wire wrapping.The effects of the thread tensile force of the steel winding in the form of single rings at the crack edges and the wires with different winding diameters and pitches were also studied.The results showed that the strengthening was preferably executed at a minimum value of the thread tensile force,which was 6.4%more effective than that at its maximum value.The analysis of the influence of the winding dia-meters showed that the equivalent stresses increased by 32%from the beginning of the crack growth until the wire broke.The increment in winding diameter decelerated the disclosure of the edge crack and reduced its length by 8.2%.The analysis of the influence of the winding pitch showed that decreasing the distance between the winding turns also led to a 33.6%reduction in the length of the straight crack and a 7.9%reduction in the maximum stres-ses on the strengthened pipeline cross-section.The analysis of the temperature effect on the pipeline material,within a range from-40℃to+50℃,resulted in a crack length change of up to 5.8%.As the temperature dropped,the crack length decreased.Within such a temperature range,the maximum stresses were observed along the cen-tral area of the crack,which were equal to 413 MPa at+50℃and 440 MPa at-40℃.The results also showed that the presence of the steel winding in the pipeline significantly reduced the length of crack propagation up to 8.4 times,depending on the temperature effect and design parameters of prestressing.This work integrated the existing methods for crack localization along steel gas pipelines.展开更多
基金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.
基金the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA14020102)the National Natural Science Foundation of China(Nos.41774125,41530320,41904104)+1 种基金the Key National Research Project of China(No.2018YFC0603300)the S&T Program of Beijing(No.Z181100005718001).
文摘We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.
基金funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan(Grant No.AP19680589).
文摘Using the software ANSYS-19.2/Explicit Dynamics,this study performedfinite-element modeling of the large-diameter steel pipeline cross-section for the Beineu-Bozoy-Shymkent gas pipeline with a non-through straight crack,strengthened by steel wire wrapping.The effects of the thread tensile force of the steel winding in the form of single rings at the crack edges and the wires with different winding diameters and pitches were also studied.The results showed that the strengthening was preferably executed at a minimum value of the thread tensile force,which was 6.4%more effective than that at its maximum value.The analysis of the influence of the winding dia-meters showed that the equivalent stresses increased by 32%from the beginning of the crack growth until the wire broke.The increment in winding diameter decelerated the disclosure of the edge crack and reduced its length by 8.2%.The analysis of the influence of the winding pitch showed that decreasing the distance between the winding turns also led to a 33.6%reduction in the length of the straight crack and a 7.9%reduction in the maximum stres-ses on the strengthened pipeline cross-section.The analysis of the temperature effect on the pipeline material,within a range from-40℃to+50℃,resulted in a crack length change of up to 5.8%.As the temperature dropped,the crack length decreased.Within such a temperature range,the maximum stresses were observed along the cen-tral area of the crack,which were equal to 413 MPa at+50℃and 440 MPa at-40℃.The results also showed that the presence of the steel winding in the pipeline significantly reduced the length of crack propagation up to 8.4 times,depending on the temperature effect and design parameters of prestressing.This work integrated the existing methods for crack localization along steel gas pipelines.
