Drilling and blasting are the two most significant operations in open pit mines that play a crucial role in downstream stages. While previous research has focused on optimizing these operations as two separate parts o...Drilling and blasting are the two most significant operations in open pit mines that play a crucial role in downstream stages. While previous research has focused on optimizing these operations as two separate parts or merely in a specific parameter, this paper proposes a system dynamic model(SDM) for drilling and blasting operations as an interactive system. In addition, some technical and economic uncertainties such as rock density, uniaxial compressive strength, bit life and operating costs are considered in this system to evaluate the different optimization results. For this purpose, Vensim simulation software is utilized as a powerful dynamic tool for both modelling and optimizing under deterministic and uncertain conditions. It is concluded that an integrated optimization as opposed to the deterministic approach can be efficiently achieved. This however is dependent on the parameters that are considered as uncertainties.展开更多
It is difficult to determine the optimal weight on bit (WOB) of the double-driven bottomhole assembly (DD-BHA, with double stabilizers and a bent housing positive displacement motor (PDM)) which is employed during ver...It is difficult to determine the optimal weight on bit (WOB) of the double-driven bottomhole assembly (DD-BHA, with double stabilizers and a bent housing positive displacement motor (PDM)) which is employed during vertical and fast drilling. High WOB leads to well deviation out of control, and low WOB leads to low rate of penetration (ROP). So considering the rock physical properties, the anisotropy index function of polycrystalline diamond compact (PDC) bit was derived with the structure and cutting performance parameters of the bit, and the effect of natural hole deviation tendencies on the performance of DD-BHA resisting deviation was represented. The concept of elliptic deformation ratio was used to characterize the performance of DD-BHA resisting deviation. Eventually, a model calculating the critical WOB was established. By comparing the model predictions with the measured hole angle changes in the field, the results show that the model predictions are accurate with error less than 5.8%, which can meet the operational requirements in the projects. Furthermore, the model was adopted to justify and guide the operating conditions and parameters during drilling, which shows that the optimum WOB predicted by the model can not only control deviation but also improve ROP effectively. The model is independent on the formation characteristics of blocks, so it can be expanded widely to other oilfields.展开更多
In drilling operation,a large saving in time and money would be achieved by reducing the drilling time,since some of the costs are time-dependent.Drilling time could be minimized by raising the penetration rate.In the...In drilling operation,a large saving in time and money would be achieved by reducing the drilling time,since some of the costs are time-dependent.Drilling time could be minimized by raising the penetration rate.In the comparative optimization method,by using the records of the first drilled wells and comparing the criteria like penetration rate,cost per foot and specific energy,the drilling parameters of the next wells being drilled can be optimized in each depth interval.In the mathematical optimization technique,some numerical equations to model the penetration rate,bit wear rate and hydraulics would be used to minimize the drilling cost and time as much as possible and improve the results of the primary comparative optimization.In this research,as a case study the Iranian Khangiran gas field has been evaluated to optimize the drilling costs.A combination of the mentioned optimization techniques resulted in an optimal well which reduced the drilling time and cost considerably in comparison with the wells already drilled.展开更多
To optimize cutting control parameters and provide scientific evidence for controlling cutting forces,cutting force modeling and cutting control parameter optimization are researched with one tool adopted to orbital d...To optimize cutting control parameters and provide scientific evidence for controlling cutting forces,cutting force modeling and cutting control parameter optimization are researched with one tool adopted to orbital drill holes in aluminum alloy 6061.Firstly,four cutting control parameters(tool rotation speed,tool revolution speed,axial feeding pitch and tool revolution radius)and affecting cutting forces are identified after orbital drilling kinematics analysis.Secondly,hybrid level orthogonal experiment method is utilized in modeling experiment.By nonlinear regression analysis,two quadratic prediction models for axial and radial forces are established,where the above four control parameters are used as input variables.Then,model accuracy and cutting control parameters are analyzed.Upon axial and radial forces models,two optimal combinations of cutting control parameters are obtained for processing a13mm hole,corresponding to the minimum axial force and the radial force respectively.Finally,each optimal combination is applied in verification experiment.The verification experiment results of cutting force are in good agreement with prediction model,which confirms accracy of the research method in practical production.展开更多
During the last three decades,evolutionary algorithms(EAs)have shown superiority in solving complex optimization problems,especially those with multiple objectives and non-differentiable landscapes.However,due to the ...During the last three decades,evolutionary algorithms(EAs)have shown superiority in solving complex optimization problems,especially those with multiple objectives and non-differentiable landscapes.However,due to the stochastic search strategies,the performance of most EAs deteriorates drastically when handling a large number of decision variables.To tackle the curse of dimensionality,this work proposes an efficient EA for solving super-large-scale multi-objective optimization problems with sparse optimal solutions.The proposed algorithm estimates the sparse distribution of optimal solutions by optimizing a binary vector for each solution,and provides a fast clustering method to highly reduce the dimensionality of the search space.More importantly,all the operations related to the decision variables only contain several matrix calculations,which can be directly accelerated by GPUs.While existing EAs are capable of handling fewer than 10000 real variables,the proposed algorithm is verified to be effective in handling 1000000 real variables.Furthermore,since the proposed algorithm handles the large number of variables via accelerated matrix calculations,its runtime can be reduced to less than 10%of the runtime of existing EAs.展开更多
The purpose of this study was to grasp current potential problems of dose error in intensity-modulated proton therapy (IMPT) plans. We were interested in dose differences of the Varian Eclipse treatment planning syste...The purpose of this study was to grasp current potential problems of dose error in intensity-modulated proton therapy (IMPT) plans. We were interested in dose differences of the Varian Eclipse treatment planning system (TPS) and the fast dose calculation method (FDC) for single-field optimization (SFO) and multi-field optimization (MFO) IMPT plans. In addition, because some authors have reported dosimetric benefit of a proton arc therapy with ultimate multi-fields in recent years, we wanted to evaluate how the number of fields and beam angles affect the differences for IMPT plans. Therefore, for one brain cancer patient with a large heterogeneity, SFO and MFO IMPT plans with various multi-angle beams were planned by the TPS. Dose distributions for each IMPT plan were calculated by both the TPS’s conventional pencil beam algorithm and the FDC. The dosimetric parameters were compared between the two algorithms. The TPS overestimated 400 - 500 cGy (RBE) for minimum dose to the CTV relative to the dose calculated by the FDC. These differences indicate clinically relevant effect on clinical results. In addition, we observed that the maximum difference in dose calculated between the TPS and the FDC was about 900 cGy (RBE) for the right optic nerve, and this quantity also has a possibility to have a clinical effect. The major difference was not seen in calculations for SFO IMPT planning and those for MFO IMPT planning. Differences between the TPS and the FDC in SFO and MFO IMPT plans depend strongly on beam arrangement and the presence of a heterogeneous body. We advocate use of a Monte Carlo method in proton treatment planning to deliver the most precise proton dose in IMPT.展开更多
Diamond turning based on a fast tool servo(FTS)is widely used in freeform optics fabrication due to its high accuracy and machining efficiency.As a new trend,recently developed high-frequency and long-stroke FTS units...Diamond turning based on a fast tool servo(FTS)is widely used in freeform optics fabrication due to its high accuracy and machining efficiency.As a new trend,recently developed high-frequency and long-stroke FTS units are independently driven by a separate control system from the machine tool controller.However,the tool path generation strategy for the independently controlled FTS is far from complete.This study aims to establish methods for optimizing tool path for the independent control FTS to reduce form errors in a single step of machining.Different from the conventional integrated FTS control system,where control points are distributed in a spiral pattern,in this study,the tool path for the independent FTS controller is generated by the ring method and the mesh method,respectively.The machined surface profile is predicted by simulation and the parameters for the control point generation are optimized by minimizing the deviation between the predicted and the designed surfaces.To demonstrate the feasibility of the proposed tool path generation strategies,cutting tests of a two-dimensional sinewave and a micro-lens array were conducted and the results were compared.As a result,after tool path optimization,the peak-to-valley form error of the machined surface was reduced from 429 nm to 56 nm for the two-dimensional sinewave by using the ring method,and from 191 nm to 103 nm for the micro-lens array by using the mesh method,respectively.展开更多
Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex str...Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.展开更多
文摘Drilling and blasting are the two most significant operations in open pit mines that play a crucial role in downstream stages. While previous research has focused on optimizing these operations as two separate parts or merely in a specific parameter, this paper proposes a system dynamic model(SDM) for drilling and blasting operations as an interactive system. In addition, some technical and economic uncertainties such as rock density, uniaxial compressive strength, bit life and operating costs are considered in this system to evaluate the different optimization results. For this purpose, Vensim simulation software is utilized as a powerful dynamic tool for both modelling and optimizing under deterministic and uncertain conditions. It is concluded that an integrated optimization as opposed to the deterministic approach can be efficiently achieved. This however is dependent on the parameters that are considered as uncertainties.
基金Projects(50804040,51004082)supported by the National Natural Science Foundation of China
文摘It is difficult to determine the optimal weight on bit (WOB) of the double-driven bottomhole assembly (DD-BHA, with double stabilizers and a bent housing positive displacement motor (PDM)) which is employed during vertical and fast drilling. High WOB leads to well deviation out of control, and low WOB leads to low rate of penetration (ROP). So considering the rock physical properties, the anisotropy index function of polycrystalline diamond compact (PDC) bit was derived with the structure and cutting performance parameters of the bit, and the effect of natural hole deviation tendencies on the performance of DD-BHA resisting deviation was represented. The concept of elliptic deformation ratio was used to characterize the performance of DD-BHA resisting deviation. Eventually, a model calculating the critical WOB was established. By comparing the model predictions with the measured hole angle changes in the field, the results show that the model predictions are accurate with error less than 5.8%, which can meet the operational requirements in the projects. Furthermore, the model was adopted to justify and guide the operating conditions and parameters during drilling, which shows that the optimum WOB predicted by the model can not only control deviation but also improve ROP effectively. The model is independent on the formation characteristics of blocks, so it can be expanded widely to other oilfields.
文摘In drilling operation,a large saving in time and money would be achieved by reducing the drilling time,since some of the costs are time-dependent.Drilling time could be minimized by raising the penetration rate.In the comparative optimization method,by using the records of the first drilled wells and comparing the criteria like penetration rate,cost per foot and specific energy,the drilling parameters of the next wells being drilled can be optimized in each depth interval.In the mathematical optimization technique,some numerical equations to model the penetration rate,bit wear rate and hydraulics would be used to minimize the drilling cost and time as much as possible and improve the results of the primary comparative optimization.In this research,as a case study the Iranian Khangiran gas field has been evaluated to optimize the drilling costs.A combination of the mentioned optimization techniques resulted in an optimal well which reduced the drilling time and cost considerably in comparison with the wells already drilled.
基金Supported by the National Natural Science Foundation of China(50975141)the Aviation Science Fund(20091652018,2010352005)the National Science and Technology Major Project of the Ministry of Science and Technology of China(2012ZX04003031-4)
文摘To optimize cutting control parameters and provide scientific evidence for controlling cutting forces,cutting force modeling and cutting control parameter optimization are researched with one tool adopted to orbital drill holes in aluminum alloy 6061.Firstly,four cutting control parameters(tool rotation speed,tool revolution speed,axial feeding pitch and tool revolution radius)and affecting cutting forces are identified after orbital drilling kinematics analysis.Secondly,hybrid level orthogonal experiment method is utilized in modeling experiment.By nonlinear regression analysis,two quadratic prediction models for axial and radial forces are established,where the above four control parameters are used as input variables.Then,model accuracy and cutting control parameters are analyzed.Upon axial and radial forces models,two optimal combinations of cutting control parameters are obtained for processing a13mm hole,corresponding to the minimum axial force and the radial force respectively.Finally,each optimal combination is applied in verification experiment.The verification experiment results of cutting force are in good agreement with prediction model,which confirms accracy of the research method in practical production.
基金This work was supported in part by the National Key Research and Development Program of China(2018AAA0100100)the National Natural Science Foundation of China(61822301,61876123,61906001)+2 种基金the Collaborative Innovation Program of Universities in Anhui Province(GXXT-2020-051)the Hong Kong Scholars Program(XJ2019035)Anhui Provincial Natural Science Foundation(1908085QF271).
文摘During the last three decades,evolutionary algorithms(EAs)have shown superiority in solving complex optimization problems,especially those with multiple objectives and non-differentiable landscapes.However,due to the stochastic search strategies,the performance of most EAs deteriorates drastically when handling a large number of decision variables.To tackle the curse of dimensionality,this work proposes an efficient EA for solving super-large-scale multi-objective optimization problems with sparse optimal solutions.The proposed algorithm estimates the sparse distribution of optimal solutions by optimizing a binary vector for each solution,and provides a fast clustering method to highly reduce the dimensionality of the search space.More importantly,all the operations related to the decision variables only contain several matrix calculations,which can be directly accelerated by GPUs.While existing EAs are capable of handling fewer than 10000 real variables,the proposed algorithm is verified to be effective in handling 1000000 real variables.Furthermore,since the proposed algorithm handles the large number of variables via accelerated matrix calculations,its runtime can be reduced to less than 10%of the runtime of existing EAs.
文摘The purpose of this study was to grasp current potential problems of dose error in intensity-modulated proton therapy (IMPT) plans. We were interested in dose differences of the Varian Eclipse treatment planning system (TPS) and the fast dose calculation method (FDC) for single-field optimization (SFO) and multi-field optimization (MFO) IMPT plans. In addition, because some authors have reported dosimetric benefit of a proton arc therapy with ultimate multi-fields in recent years, we wanted to evaluate how the number of fields and beam angles affect the differences for IMPT plans. Therefore, for one brain cancer patient with a large heterogeneity, SFO and MFO IMPT plans with various multi-angle beams were planned by the TPS. Dose distributions for each IMPT plan were calculated by both the TPS’s conventional pencil beam algorithm and the FDC. The dosimetric parameters were compared between the two algorithms. The TPS overestimated 400 - 500 cGy (RBE) for minimum dose to the CTV relative to the dose calculated by the FDC. These differences indicate clinically relevant effect on clinical results. In addition, we observed that the maximum difference in dose calculated between the TPS and the FDC was about 900 cGy (RBE) for the right optic nerve, and this quantity also has a possibility to have a clinical effect. The major difference was not seen in calculations for SFO IMPT planning and those for MFO IMPT planning. Differences between the TPS and the FDC in SFO and MFO IMPT plans depend strongly on beam arrangement and the presence of a heterogeneous body. We advocate use of a Monte Carlo method in proton treatment planning to deliver the most precise proton dose in IMPT.
基金supported by Japan Society for the Promotion of Science,Grant-in-Aid for Scientific Research(B),Project Number 21H01230.
文摘Diamond turning based on a fast tool servo(FTS)is widely used in freeform optics fabrication due to its high accuracy and machining efficiency.As a new trend,recently developed high-frequency and long-stroke FTS units are independently driven by a separate control system from the machine tool controller.However,the tool path generation strategy for the independently controlled FTS is far from complete.This study aims to establish methods for optimizing tool path for the independent control FTS to reduce form errors in a single step of machining.Different from the conventional integrated FTS control system,where control points are distributed in a spiral pattern,in this study,the tool path for the independent FTS controller is generated by the ring method and the mesh method,respectively.The machined surface profile is predicted by simulation and the parameters for the control point generation are optimized by minimizing the deviation between the predicted and the designed surfaces.To demonstrate the feasibility of the proposed tool path generation strategies,cutting tests of a two-dimensional sinewave and a micro-lens array were conducted and the results were compared.As a result,after tool path optimization,the peak-to-valley form error of the machined surface was reduced from 429 nm to 56 nm for the two-dimensional sinewave by using the ring method,and from 191 nm to 103 nm for the micro-lens array by using the mesh method,respectively.
基金We would like to acknowledge the sponsorship of the National Natural Science Foundation of China(42004092,42030103,41974119)Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(Grant No.2021QNLM020001-6)Young Elite Scientists Sponsorship Program by CAST(2021QNRC001).
文摘Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.
