<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">EE</journal-id>
<journal-id journal-id-type="nlm-ta">EE</journal-id>
<journal-id journal-id-type="publisher-id">EE</journal-id>
<journal-title-group>
<journal-title>Energy Engineering</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-0118</issn>
<issn pub-type="ppub">0199-8595</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">43497</article-id>
<article-id pub-id-type="doi">10.32604/ee.2023.043497</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Research on Carbon Emission for Preventive Maintenance of Wind Turbine Gearbox Based on Stochastic Differential Equation</article-title>
<alt-title alt-title-type="left-running-head">Research on Carbon Emission for Preventive Maintenance of Wind Turbine Gearbox Based on Stochastic Differential Equation</alt-title>
<alt-title alt-title-type="right-running-head">Research on Carbon Emission for Preventive Maintenance of Wind Turbine Gearbox Based on Stochastic Differential Equation</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Su</surname><given-names>Hongsheng</given-names></name></contrib>
<contrib id="author-2" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Dong</surname><given-names>Lixia</given-names></name><email>donglixia0425@163.com</email></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Yu</surname><given-names>Xiaoying</given-names></name></contrib>
<contrib id="author-4" contrib-type="author">
<name name-style="western"><surname>Liu</surname><given-names>Kai</given-names></name></contrib>
<aff><institution>School of Automation and Electrical Engineering, Lanzhou Jiaotong University</institution>, <addr-line>Lanzhou, 730070</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Lixia Dong. Email: <email>donglixia0425@163.com</email></corresp>
</author-notes>
<pub-date date-type="collection" publication-format="electronic"><year>2024</year></pub-date>
<pub-date date-type="pub" publication-format="electronic"><day>26</day><month>3</month><year>2024</year></pub-date>
<volume>121</volume>
<issue>4</issue>
<fpage>973</fpage>
<lpage>986</lpage>
<history>
<date date-type="received">
<day>04</day><month>7</month><year>2023</year></date>
<date date-type="accepted">
<day>09</day><month>10</month><year>2023</year></date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2023 Su et al.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Su et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_EE_43497.pdf"></self-uri>
<abstract>
<p>Time based maintenance (TBM) and condition based maintenance (CBM) are widely applied in many large wind farms to optimize the maintenance issues of wind turbine gearboxes, however, these maintenance strategies do not take into account environmental benefits during full life cycle such as carbon emissions issues. Hence, this article proposes a carbon emissions computing model for preventive maintenance activities of wind turbine gearboxes to solve the issue. Based on the change of the gearbox state during operation and the influence of external random factors on the gearbox state, a stochastic differential equation model (SDE) and corresponding carbon emission model are established, wherein SDE is applied to model the evolution of the device state, whereas carbon emission is used to implement carbon emissions computing. The simulation results indicate that the proposed preventive maintenance cannot ensure reliable operation of wind turbine gearboxes but reduce carbon emissions during their lifespan. Compared with TBM, CBM minimizes unit carbon emissions without influencing reliable operation, making it an effective maintenance method.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Stochastic differential equation (SDE)</kwd>
<kwd>condition-based maintenance (CBM)</kwd>
<kwd>carbon emissions</kwd>
</kwd-group>
<funding-group>
<award-group id="awg1">
<funding-source>Basic Science Research Program through the National Natural Science Foundation of China</funding-source>
<award-id>61867003</award-id>
</award-group>
<award-group id="awg2">
<funding-source>Key Project of Science and Technology Research and Development Plan of China Railway Co., Ltd.</funding-source>
<award-id>N2022X009</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<p>Wind power generation is one of the most mature technologies in renewable energy. However, wind farms are often situated in remote locations characterized by harsh working environments and elevated equipment heights, which significantly challenge the maintenance of wind turbine equipment. Wind turbines, as critical components of the wind power system, possess a complex structure and operate under demanding conditions. They are subject to various loads and environmental influences, leading to component degradation, faults, and failures. Particularly, gearboxes, as the components with the highest failure rate in wind turbines, frequently experience malfunctions and shutdowns. Consequently, reliability analysis is imperative to optimize maintenance strategies. Preventive maintenance is a pivotal method for ensuring equipment operates in a stable and efficient condition. It involves proactive equipment monitoring to identify early signs of failure and implementing maintenance activities to prevent equipment degradation, thus reducing downtime and maintenance losses. Traditional preventive maintenance strategies, such as time-based maintenance (TBM), often result in over or under-maintenance scenarios [<xref ref-type="bibr" rid="ref-1">1</xref>&#x2013;<xref ref-type="bibr" rid="ref-3">3</xref>].</p>
<p>Condition-based maintenance (CBM) represents a more advanced approach, focusing on monitoring and diagnosing equipment conditions to detect abnormalities and tailor maintenance plans accordingly. This strategy has gained widespread application in wind farms. In preventive maintenance models for wind turbines, the Weibull model is frequently employed to depict the reliability changes of wind turbines. Service age regression and failure rate increment factors are used to represent the failure recovery rate and change rate post-maintenance, respectively [<xref ref-type="bibr" rid="ref-4">4</xref>]. This approach is particularly advantageous in the maintenance strategy of units under unique offshore environmental conditions, optimizing maintenance scope primarily based on system performance. The impact of maintenance on reliability, including the correlation between lifecycle reliability and the timing of maintenance and failure, has been studied through simulation in reference [<xref ref-type="bibr" rid="ref-5">5</xref>]. A proportional hazard model, describing the relationship between the operating status of a gearbox and its failure rate, has been utilized to predict the maintenance interval of gearboxes, with reliability as the core maintenance objective [<xref ref-type="bibr" rid="ref-6">6</xref>]. Therefore, the operating status of a gearbox is crucial in determining the optimal inspection time. Su et al. explored the operation and reliability of gearboxes under preventive maintenance, considering both TBM and CBM, by establishing a unified maintenance model based on stochastic process theory [<xref ref-type="bibr" rid="ref-7">7</xref>&#x2013;<xref ref-type="bibr" rid="ref-9">9</xref>].</p>
<p>Previous literature has focused on predicting gearbox maintenance intervals with reliability as the primary objective. However, it often overlooks the environmental impacts on gearbox condition and the associated carbon emissions throughout the gearbox&#x2019;s lifecycle. Prolonged operation of wind equipment not only increases preventive maintenance costs but also leads to performance deterioration. The operation of wind turbines is inherently linked with escalated energy consumption, operating costs, and carbon emissions. Thus, implementing appropriate maintenance decisions can significantly reduce both carbon emissions and operational costs. Maintenance decisions, as opposed to technological advancements, present a more feasible approach for reducing carbon emissions.</p>
<p>Franciosi et al. integrated the concept of circular economy into the regular preventive maintenance model, highlighting the importance of considering sustainability factors in daily maintenance activities for sustainable development [<xref ref-type="bibr" rid="ref-10">10</xref>]. A joint optimization model for maintenance and production plans, considering carbon emissions, was developed to determine the optimal maintenance and production strategy, factoring in energy consumption per unit of equipment and resultant carbon emissions [<xref ref-type="bibr" rid="ref-11">11</xref>]. Afrinaldi et al. proposed a model for determining optimal preventive replacements to minimize economic and environmental impacts [<xref ref-type="bibr" rid="ref-12">12</xref>]. A maintenance optimization model that considers carbon emissions was introduced in reference [<xref ref-type="bibr" rid="ref-13">13</xref>], incorporating improvement factors into the equipment failure rate and carbon emissions model to highlight maintenance imperfections. Subsequently, a multi-objective decision model was developed to minimize both carbon emissions and cost rates. Liu et al. proposed an opportunistic maintenance strategy for wind turbines, accounting for structural and random correlations, as well as carbon emissions. This strategy involves introducing service age regression factors and failure rate increasing factors into the failure rate and carbon emissions models to depict incomplete maintenance. Through case analysis, the optimal opportunistic maintenance timing was determined, elucidating the correlation coefficient and sensitivity of carbon emissions [<xref ref-type="bibr" rid="ref-14">14</xref>]. Although these studies consider carbon emissions during equipment maintenance and introduce factors to describe maintenance imperfections, they often neglect the impact of random minor factors on equipment condition.</p>
<p>There is a discernible relationship between equipment maintenance activities and carbon emissions. The assessment of equipment condition directly influences the choice of maintenance strategies. In conducting preventive maintenance, timely adjustment of maintenance schedules is essential to safeguard environmental benefits and minimize carbon emissions. Therefore, this article focuses on the wind turbine gearbox, considering the impact of external random disturbances on equipment state and the reduction in equipment failure and emission rates due to maintenance activities. Additionally, the role of equipment recycling in emission reduction is examined. By optimizing the preventive maintenance plan for gearboxes, the environmental benefits of the equipment are balanced, and the model&#x2019;s effectiveness is demonstrated through numerical examples.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Mathematical Model and Solution</title>
<sec id="s2_1">
<label>2.1</label>
<title>Description of Gearbox Model</title>
<p>Current equipment state models primarily utilize ordinary differential equations, which are based on fault rates but often neglect the impact of external disturbances during operation. Given that wind turbine equipment operation is influenced by random factors such as environmental conditions, weather changes, and routine inspections, it is essential to consider these aspects in equipment state modeling. Accordingly, to account for non-homogeneous random disturbances acting on the equipment, the Ito-type stochastic differential equation (SDE) is employed to construct the state change model. In this model, a random function X(t) represents the health state of the gearbox at time t. The evolution of X(t) reflects the state changes over the lifespan of the wind turbine gearbox. The external random interference experienced by the gearbox during operation is denoted by {B(t)}, a standard Brownian motion, and X(t) undergoes a random differential over the interval [0, t].</p>
<p>The state transition model for the wind power gearbox is represented by:</p>
<disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>dt</mml:mtext></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<p>wherein, <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents the equipment&#x2019;s failure rate function, <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> denotes the random disturbance coefficient, and <italic>k</italic> is the state fluctuation parameter, referred to as the equipment state fluctuation rate. This model is based on the following assumptions:</p>
<p>Assumption 1: The maintenance of the gearbox is instantaneously completed, with the expected random disturbance value being zero.</p>
<p>Assumption 2: <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> are <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:mi>B</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> measurable functions defined on <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:mrow><mml:mo>[</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, if <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:mi>&#x03BC;</mml:mi><mml:mo>&#x2264;</mml:mo><mml:mi>N</mml:mi></mml:math></inline-formula> and all <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:mi>t</mml:mi><mml:mo>&#x2264;</mml:mo><mml:mi>T</mml:mi></mml:math></inline-formula> meet the following conditions:
<disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mrow><mml:mo>|</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x003C;</mml:mo><mml:mi>K</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:mi>x</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>y</mml:mi><mml:mo>|</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:mrow><mml:mo>|</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x003C;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mrow><mml:mtext>T</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mi>x</mml:mi><mml:mo>|</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Under these assumptions, the solution to the stochastic differential equation possesses a unique strong solution.</p>
<p>In life data analysis, the Weibull distribution is extensively used due to its ability to fit a wide range of sample data by adjusting its shape, scale, and position parameters. The variability in its shape parameters enables the description of diverse data trends, including decline, stability, and growth. The degradation of gearboxes aligns with the Weibull distribution, effectively characterizing the life distribution of wind turbine gearboxes. Based on this, the gearbox model assumes failure through component fatigue, with its shape parameters offering considerable flexibility in data fitting. Therefore, the fundamental failure rate model for gearboxes can be formulated as:
<disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:msub><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mi>&#x03B2;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula>wherein, <italic>&#x03B2;</italic> is a shape parameter, and <italic>&#x03B7;</italic> is a proportion parameter.</p>
<p>The changes in the state of the wind turbine gearbox are influenced by a combination of the inherent degradation of the gearbox and external factors. The failure rate of the wind turbine gearbox is a continuous function, representing the basic failure rate. Consequently, the comprehensive failure rate of the wind power gearbox can be determined by incorporating the state fluctuation rate [<xref ref-type="bibr" rid="ref-15">15</xref>].
<disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mi>&#x03B2;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>The reliability function of the gearbox is denoted by <xref ref-type="disp-formula" rid="eqn-7">Eq. (7)</xref>:
<disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mo>&#x2212;</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mo>&#x2212;</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>The <italic>&#x03B6;</italic>-th preventive maintenance interval of the gearbox is denoted by <italic>T</italic><sub><italic>&#x03B6;</italic></sub>. Therefore, <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><mml:mo stretchy="false">[</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula> is used to represent each preventive maintenance interval within the lifecycle of the gearbox.</p>
<p>Consequently, the number of preventive maintenance interventions can be expressed as:
<disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03C2;</mml:mi></mml:mrow></mml:munderover><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></disp-formula></p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Model Analysis</title>
<p>Based on the failure rate model for the wind power gearbox, the likelihood function is formulated using <xref ref-type="disp-formula" rid="eqn-7">Eq. (7)</xref> in conjunction with the failure density function <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> as follows:
<disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B2;</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03B7;</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>=</mml:mo><mml:mspace width="negativethinmathspace" /><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:munderover><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>=</mml:mo><mml:mspace width="negativethinmathspace" /><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:munderover><mml:mfrac><mml:mi>&#x03B2;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>&#x22C5;</mml:mo><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mspace width="negativethinmathspace" /><mml:mo>&#x2212;</mml:mo><mml:mspace width="negativethinmathspace" /><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mspace width="negativethinmathspace" /><mml:mo>&#x2212;</mml:mo><mml:mspace width="negativethinmathspace" /><mml:mi>&#x03B1;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>The logarithmic Likelihood function is then derived, as illustrated in <xref ref-type="disp-formula" rid="eqn-10">Eq. (10)</xref>:
<disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B2;</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03B7;</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mfrac><mml:mi>&#x03B2;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>+</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B2;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:mo>[</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>The Newton-Raphson iteration method is employed to develop the following iteration formula:
<disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:mi>M</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B1;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mtable columnalign="left left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:msup><mml:mi>&#x03B2;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B2;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B2;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B1;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B7;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:msup><mml:mi>&#x03B7;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B7;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B1;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B1;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B2;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B1;</mml:mi><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>ln</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:msup><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>This yields:
<disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mtable rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mi>&#x03B2;</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B7;</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B1;</mml:mi></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mi>g</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mtable rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mi>&#x03B2;</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B7;</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B1;</mml:mi></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>By solving <xref ref-type="disp-formula" rid="eqn-13">Eq. (13)</xref>, the parameter estimation values for the above fault rate model can be determined. The maintenance timing for CBM is deduced based on a reliability threshold. It is assumed that at time t, the reliability of the gearbox equipment either reaches or falls below a threshold, denoted as <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, When this occurs, maintenance is required to ensure that the condition of the gearbox remains within the normal operating range. Herein, <italic>R</italic><sub><italic>ref</italic></sub> represents the reliability threshold of the device.</p>
<sec id="s2_2_1">
<label>2.2.1</label>
<title>TBM Analysis</title>
<p>TBM assumes a consistent expected sample in each cycle [<xref ref-type="bibr" rid="ref-16">16</xref>]. An ordinary differential equation (ODE) is utilized to construct the state transition model for the wind power gearbox, as delineated below:
<disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:mrow><mml:mtext>dE</mml:mtext></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>t</mml:mi><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>t</mml:mi></mml:math></disp-formula></p>
<p>The state transition model, developed using ODE, is:
<disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mspace width="negativethinmathspace" /><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mi>&#x03B2;</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mi>&#x03B7;</mml:mi></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B2;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B4;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>t</mml:mi></mml:math></disp-formula>wherein, <italic>&#x03B4;</italic> represents the regression coefficient. Since <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><mml:mrow><mml:mover><mml:mrow><mml:mtext>X</mml:mtext></mml:mrow><mml:mo stretchy="false">&#x005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the average expected value of X(t), the maximum likelihood estimation method is then applied to solve <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref>.</p>
</sec>
</sec>
</sec>
<sec id="s3">
<label>3</label>
<title>Carbon Emission Model</title>
<p>Throughout the lifecycle of wind turbines, environmental impacts primarily involve carbon emissions resulting from energy consumption during use and recycling. Recycling is recognized as a means to reduce carbon emissions. The need for maintenance is primarily due to the continuous degradation and performance deterioration of wind turbines during operation [<xref ref-type="bibr" rid="ref-17">17</xref>].</p>
<p>In operation, the gearbox transmits power from the wind to the generator. The rotor&#x2019;s low speed in the wind turbine must be accelerated by the gearbox to achieve the necessary speed for power generation. To ensure braking capacity, braking devices are installed at both the input and output ends of the gearbox. These work in conjunction with the braking system to halt the wind turbine drive system. Consequently, during operation, forced lubrication of gears and bearings is critical for safety. This lubrication system comprises an oil supply device, a filtration system, an oil/air cooling device, and connecting pipelines. The energy consumption of this lubrication system is a contributor to carbon emissions. The primary objectives of lubricating the gearbox are to reduce friction and provide cooling, thereby ensuring efficient operation and extending the service life of the gearbox [<xref ref-type="bibr" rid="ref-18">18</xref>].</p>
<p>As such, the gearbox consumes energy during operation, leading to carbon emissions. For every 1 kWh of electricity consumed, carbon emissions amount to 0.272 kg, and carbon dioxide emissions are 0.997 kg.</p>
<p>When sudden malfunctions occur in a wind turbine&#x2019;s gearbox, minor repairs are undertaken to maintain normal operation. However, these repairs neither alter the gearbox&#x2019;s failure rate nor its carbon emissions; they simply restore operational health. Increasing the frequency of preventive maintenance enhances the duration of reliable gearbox operation, reducing the likelihood of failures and the need for minor maintenance.</p>
<p>Prolonging the operation of the wind gearbox escalates energy consumption and environmental impact, accompanied by operational performance degradation. Throughout the gearbox&#x2019;s lifecycle, increasing the frequency of preventive maintenance for wind turbine equipment reduces the probability of unexpected gearbox failures and accordingly diminishes the frequency of fault maintenance. However, excessive preventive maintenance within a cycle can lead to unnecessary resource waste and an increase in carbon emissions during that cycle, as illustrated in <xref ref-type="fig" rid="fig-1">Fig. 1</xref> [<xref ref-type="bibr" rid="ref-14">14</xref>].</p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>Variation in carbon emissions relative to the number of preventive maintenance interventions and preventive maintenance intervals</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-1.tif"/>
</fig>
<p>During the operational phase of a wind power project, which is the most prolonged phase of the entire project lifecycle, activities predominantly involve replacing equipment and materials and consuming energy resources. The carbon emissions during the gearbox&#x2019;s lifespan are mainly due to energy used in operation, maintenance, construction, installation, and retirement of the wind turbine equipment [<xref ref-type="bibr" rid="ref-14">14</xref>].</p>
<p>CO<sub>2</sub> emissions and energy consumption of wind farms are largely concentrated in the production and transportation stage of wind turbines, accounting for 89.62% of the total CO<sub>2</sub> emissions and 113.95% of the total energy consumption. CO<sub>2</sub> emissions and energy consumption in the infrastructure and operation and maintenance stages are the second and third highest, respectively. Recycling in the waste disposal stage results in CO<sub>2</sub> emissions and energy consumption accounting for 46.24% and &#x2212;88.72% of the total in the gearbox retirement stage, respectively [<xref ref-type="bibr" rid="ref-19">19</xref>]. <xref ref-type="fig" rid="fig-2">Fig. 2</xref> shows the distribution of CO<sub>2</sub> emissions and energy consumption across these four stages.</p>
<fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>Ratio of CO<sub>2</sub> emissions and energy consumption in each stage of the wind farm</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-2.tif"/>
</fig>
<p>The carbon emissions over the lifespan of gearboxes are divided into three parts: emissions generated during the operation phase, emissions during the maintenance phase, emissions during the production phase, and emissions recovered during the retirement phase.</p>
<p>Therefore, the carbon emissions <italic>GT</italic> over the lifespan of the gearbox are represented by:
<disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:mi>G</mml:mi><mml:mi>T</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:math></disp-formula>wherein <italic>G</italic><sub><italic>1</italic></sub> denotes the carbon emissions generated during the operation phase of the gearbox, equivalent to the total energy consumption during operation. This is expressed as:
<disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mi>H</mml:mi></mml:math></disp-formula></p>
<p><italic>G</italic><sub><italic>use</italic></sub>-Carbon emissions per unit of energy consumed;</p>
<p><italic>H</italic>-Total energy consumption during equipment operation.</p>
<p><italic>G</italic><sub><italic>2</italic></sub> represents the carbon emissions produced during the gearbox&#x2019;s maintenance phase, calculated as the total number of maintenance occurrences multiplied by the carbon emissions from a single maintenance event.
<disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>&#x03B6;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><italic>&#x03B6;</italic>-The total number of maintenance times within the gearbox&#x2019;s lifecycle;</p>
<p><italic>G</italic><sub><italic>m</italic></sub>-The carbon emissions generated by a single gearbox maintenance activity.</p>
<p>Furthermore, <italic>G</italic><sub><italic>3</italic></sub> refers to the carbon emissions generated during the production phase and the carbon emissions recovered during the retirement phase. As the carbon emissions recovered during the retirement phase are considered a reduction in total emissions, the carbon emissions for the production and retirement phases are calculated as the emissions generated during production minus the emissions recovered during retirement.
<disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B5;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><mml:mi>&#x03B5;</mml:mi></mml:math></inline-formula>-Recovery coefficient;</p>
<p><italic>G</italic><sub><italic>p</italic></sub>-Carbon emissions generated by manufacturing the device.</p>
<p>From the analysis above, it can be concluded that:
<disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mi>G</mml:mi><mml:mi>T</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mi>H</mml:mi><mml:mo>+</mml:mo><mml:mi>&#x03B6;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B5;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>The energy consumption function of the gearbox during operation is represented by:
<disp-formula id="eqn-21"><label>(21)</label><mml:math id="mml-eqn-21" display="block"><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>b</mml:mi></mml:math></disp-formula>wherein, <italic>a</italic> and <italic>b</italic> are parameters of the device energy consumption function. Maintenance activities not only restore the equipment&#x2019;s performance but also impact its carbon emissions. The relationship between the energy consumption functions of the gearbox, before and after the <italic>n</italic>-th preventive maintenance, is expressed as:</p>
<disp-formula id="eqn-22"><label>(22)</label><mml:math id="mml-eqn-22" display="block"><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>e</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x220F;</mml:mo><mml:mrow><mml:mi>e</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>b</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<p>where, <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> represents the rate of change in the gearbox failure rate following the <italic>e</italic>-th preventive maintenance.</p>
<p>The total energy consumption of the gearbox during operation is calculated as:
<disp-formula id="eqn-23"><label>(23)</label><mml:math id="mml-eqn-23" display="block"><mml:mi>H</mml:mi><mml:mo>=</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mo>&#x22EF;</mml:mo><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></disp-formula></p>
<p>The unit time carbon emissions over the gearbox&#x2019;s lifecycle can be represented as:
<disp-formula id="eqn-24"><label>(24)</label><mml:math id="mml-eqn-24" display="block"><mml:mi>G</mml:mi><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>G</mml:mi><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mi>H</mml:mi><mml:mo>+</mml:mo><mml:mi>&#x03B6;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B5;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03B6;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B5;</mml:mi><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
</sec>
<sec id="s4">
<label>4</label>
<title>Example Analysis</title>
<sec id="s4_1">
<label>4.1</label>
<title>Instance Validation</title>
<p>The state monitoring methods for wind power gearboxes typically include vibration monitoring, acoustic monitoring, oil grinding monitoring, and temperature monitoring, with both offline and online monitoring approaches. Under fault conditions, the vibration signal of a wind power gearbox is generally nonlinear and non-stationary [<xref ref-type="bibr" rid="ref-20">20</xref>]. Therefore, amplitude data is selected as the condition monitoring data and fault data for the wind power gearbox. This paper utilizes the actual fault monitoring data from a wind farm unit&#x2019;s gearbox for simulation. The specified service life for this type of wind turbine is 15 years, with the vibration monitoring data presented in <xref ref-type="table" rid="table-1">Table 1</xref> and some fault data shown in <xref ref-type="table" rid="table-2">Table 2</xref>.</p>
<table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>Gearbox status monitoring value</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Time/h</th>
<th>Amplitude/mm</th>
<th>Time/h</th>
<th>Amplitude/mm</th>
</tr>
</thead>
<tbody>
<tr>
<td>78</td>
<td>3.025</td>
<td>2037</td>
<td>6.744</td>
</tr>
<tr>
<td>816</td>
<td>3.081</td>
<td>2279</td>
<td>6.829</td>
</tr>
<tr>
<td>352</td>
<td>4.179</td>
<td>2493</td>
<td>7.032</td>
</tr>
<tr>
<td>609</td>
<td>4.376</td>
<td>2658</td>
<td>8.177</td>
</tr>
<tr>
<td>789</td>
<td>4.592</td>
<td>2899</td>
<td>8.853</td>
</tr>
<tr>
<td>963</td>
<td>4.802</td>
<td>3107</td>
<td>8.989</td>
</tr>
<tr>
<td>1147</td>
<td>5.032</td>
<td>3284</td>
<td>9.295</td>
</tr>
<tr>
<td>1225</td>
<td>6.272</td>
<td>3476</td>
<td>11.784</td>
</tr>
<tr>
<td>1396</td>
<td>6.458</td>
<td>3658</td>
<td>13.538</td>
</tr>
<tr>
<td>1579</td>
<td>6.597</td>
<td>3821</td>
<td>15.752</td>
</tr>
<tr>
<td>1725</td>
<td>6.621</td>
<td>3897</td>
<td>18.651</td>
</tr>
<tr>
<td>1913</td>
<td>6.685</td>
<td>4017</td>
<td>21.93</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-2">
<label>Table 2</label>
<caption>
<title>Fault data sample of gearbox</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Number</th>
<th>1</th>
<th>2</th>
<th>3</th>
<th>4</th>
<th>5</th>
<th>6</th>
<th>7</th>
<th>8</th>
<th>9</th>
<th>10</th>
</tr>
</thead>
<tbody>
<tr>
<td>Life/h</td>
<td>5327</td>
<td>6955</td>
<td>6281</td>
<td>4212</td>
<td>4406</td>
<td>6335</td>
<td>4528</td>
<td>7395</td>
<td>6382</td>
<td>4954</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Based on the state monitoring data and the Newton-Raphson model parameter solving method, the following can be deduced: <inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><mml:mi>&#x03B2;</mml:mi><mml:mo>=</mml:mo><mml:mn>2.46</mml:mn><mml:mo>,</mml:mo><mml:mi>&#x03B7;</mml:mi><mml:mo>=</mml:mo><mml:mn>733.68</mml:mn><mml:mo>,</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mn>0.335</mml:mn></mml:math></inline-formula>.</p>
<p>Consequently, the state transition model of the gearbox is formulated as follows:
<disp-formula id="eqn-25"><label>(25)</label><mml:math id="mml-eqn-25" display="block"><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mfrac><mml:mn>2.46</mml:mn><mml:mn>733.68</mml:mn></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mn>733.68</mml:mn></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>1.46</mml:mn></mml:mrow></mml:msup><mml:mo>&#x22C5;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mn>0.335</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>0.00127</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-26"><label>(26)</label><mml:math id="mml-eqn-26" display="block"><mml:mi>R</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mi>X</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mo>&#x2212;</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mi>t</mml:mi><mml:mn>733.68</mml:mn></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2.46</mml:mn></mml:mrow></mml:msup><mml:mo>&#x22C5;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mn>0.335</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mi>x</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>The state transition model of a wind power gearbox under Time-Based Maintenance (TBM) and Condition-Based Maintenance (CBM) is simulated and validated using MATLAB. Changes in transmission status and reliability over time are depicted in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>Gearbox reliability variation diagram under TBM and CBM</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-3.tif"/>
</fig>
<p>It is observed that both the operating status and reliability of the gearbox diminish over time. According to the reliability requirements for normal operation of the wind power gearbox, set at 0.9, maintenance is initiated when the reliability falls below this threshold. State-Based Maintenance (SBM) is contingent on the current health status of the wind power gearbox, with maintenance timing being variable. Upon reaching the maintenance threshold, the gearbox enters a new maintenance cycle.</p>
<p>In contrast, Time-Based Maintenance is executed after a predetermined number of operational hours. It is noted that during the first maintenance cycle, the gearbox&#x2019;s reliability did not reach the threshold before maintenance, leading to over-maintenance. Conversely, during the second cycle, the gearbox&#x2019;s reliability had already fallen below the threshold, indicative of insufficient maintenance.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Carbon Emission Analysis</title>
<sec id="s4_2_1">
<label>4.2.1</label>
<title>Carbon Emission Analysis during TBM Maintenance</title>
<p>Carbon emission analysis is performed on the gearbox during Time-Based Maintenance (TBM). As derived from <xref ref-type="disp-formula" rid="eqn-18">Eq. (18)</xref>, <italic>&#x03B6;</italic> represents the number of preventive maintenance interventions for the gearbox during TBM.</p>
<p>Namely <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><mml:mi>&#x03B6;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. Given that TBM schedules and intervals are fixed, the number of preventive maintenance cycles required over the gearbox&#x2019;s lifecycle is constant. With a TBM interval of 3500 h and a gearbox lifespan of 15 years, operating 300 days annually for an average of 13 h per day, the number of TBM maintenance cycles for the gearbox totals 17. Namely <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>17</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>s</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2720</mml:mn><mml:mo>,</mml:mo><mml:mi>H</mml:mi><mml:mo>=</mml:mo><mml:mn>108155</mml:mn></mml:math></inline-formula> <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>11450</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:mrow><mml:mtext>p</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>8</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mn>7</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>.</p>
<p>From <xref ref-type="disp-formula" rid="eqn-24">Eq. (24)</xref>, it can be inferred that:
<disp-formula id="ueqn-27"><mml:math id="mml-ueqn-27" display="block"><mml:mi>G</mml:mi><mml:mi>W</mml:mi><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>G</mml:mi><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>2720</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>108155</mml:mn><mml:mo>+</mml:mo><mml:mn>17</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>11450</mml:mn><mml:mo>+</mml:mo><mml:mn>0.2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>8</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mn>7</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mn>5216.4</mml:mn></mml:math></disp-formula></p>
<p>Hence, during TBM maintenance of the gearbox, the Global Warming Potential (GWP) of the gearbox is calculated to be 5216.4 g.</p>
</sec>
<sec id="s4_2_2">
<label>4.2.2</label>
<title>Carbon Emission Analysis during CBM Maintenance</title>
<p>For CBM maintenance of the gearbox, the number of preventive repairs <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:mi>&#x03B6;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>C</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, during the gearbox&#x2019;s lifecycle, can be determined using <xref ref-type="disp-formula" rid="eqn-8">Eqs. (8)</xref> and <xref ref-type="disp-formula" rid="eqn-24">(24)</xref>. It is found that the <italic>GWP</italic> per unit energy consumption of the gearbox is 4236.8 g, with <italic>N</italic><sub><italic>C</italic></sub> &#x003D; 11. As the number of preventive repairs varies, the GWP by the gearbox also changes correspondingly. The overall trend observed is an initial decrease followed by an increase.</p>
<p><xref ref-type="fig" rid="fig-4">Fig. 4</xref> illustrates the trend of carbon emissions <italic>GT</italic> and <italic>GWP</italic> of the gearbox over its lifespan as a function of the number of preventive repairs.</p>
<fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Changes in GWP and GT under CBM</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-4.tif"/>
</fig>
<p>The carbon emission model results for the wind turbine gearbox reveal that the minimum <italic>GWP</italic> value, 4236.8 g, is achieved when <italic>N</italic><sub><italic>C</italic></sub> &#x003D; 11.</p>
</sec>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Sensitivity Analysis</title>
<p>Given the complexity of factors influencing carbon emissions and the inherent uncertainty in the CBM process, it is crucial to understand how various impact factors may vary throughout the entire lifecycle of the equipment. This paper conducts a sensitivity analysis of three key parameters: carbon emissions per unit of energy consumption of gearboxes <italic>G</italic><sub><italic>use</italic></sub>, carbon emissions generated by maintenance activities of wind power gearboxes <italic>G</italic><sub><italic>m</italic></sub>, and carbon emissions generated during the manufacturing of a gearbox <italic>G</italic><sub><italic>p</italic></sub>.</p>
<p>(1) Carbon emissions per unit of energy consumed, <italic>G</italic><sub><italic>use</italic></sub>.</p>
<p>Variations in <italic>G</italic><sub><italic>use</italic></sub> directly influence the <italic>GWP</italic> of the gearbox.</p>
<p>As depicted in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>, when <italic>G</italic><sub><italic>use</italic></sub> is set at 3720, the minimum <italic>GWP</italic> value is achieved at <italic>N</italic><sub><italic>C</italic></sub> &#x003D; 10; conversely, when <italic>G</italic><sub><italic>use</italic></sub> is 1720, the optimal number of preventive maintenance interventions for minimizing <italic>GWP</italic> is <italic>N</italic><sub><italic>C</italic></sub> &#x003D; 14.</p>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title><italic>GWP</italic> varies with changes in <italic>G</italic><sub><italic>use</italic></sub></title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-5.tif"/>
</fig>
<p>An increase in <italic>G</italic><sub><italic>use</italic></sub> escalates the <italic>GWP</italic> of the gearbox, leading to higher carbon emissions during normal operation. Hence, adjusting the number of preventive maintenance sessions for the gearbox is a viable strategy for reducing its carbon emissions. As shown in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>, a decrease in <italic>G</italic><sub><italic>use</italic></sub> value makes the number of preventive maintenance times more sensitive to changes.</p>
<p>(2) Carbon emissions generated by a single maintenance activity <italic>G</italic><sub><italic>m</italic></sub>.</p>
<p>When the values of the parameter <italic>G</italic><sub><italic>m</italic></sub> are set at 7450, 11450, and 15450, respectively, the <italic>GWP</italic> value remains constant. Notably, when the minimum preventive maintenance frequency for optima<italic>l GWP</italic> is <italic>N</italic><sub><italic>C</italic></sub> &#x003D; 11, the variation in carbon emissions per unit operation of the gearbox is depicted in <xref ref-type="fig" rid="fig-6">Fig. 6</xref>. It is evident that the carbon emissions associated with the gearbox&#x2019;s unit energy consumption are not significantly affected by variations in <italic>G</italic><sub><italic>m</italic></sub>.</p>
<fig id="fig-6">
<label>Figure 6</label>
<caption>
<title><italic>GWP</italic> variation with changes in <italic>G</italic><sub><italic>m</italic></sub></title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-6.tif"/>
</fig>
<p>(3) The carbon emissions <italic>G</italic><sub><italic>p</italic></sub> generated by manufacturing a device.</p>
<p>The <italic>GWP</italic> responds to changes in the parameter <italic>G</italic><sub><italic>p</italic></sub> value. When <italic>G</italic><sub><italic>p</italic></sub> values are set at 1.6 &#x00D7; 10<sup>8</sup>, 8 &#x00D7; 10<sup>7</sup> and 7 &#x00D7; 10<sup>7</sup>, the preventive maintenance frequencies <italic>N</italic><sub><italic>C</italic></sub> required for the gearbox to achieve the minimum <italic>GWP</italic> value are 13, 11, and 10, respectively. As illustrated in <xref ref-type="fig" rid="fig-7">Fig. 7</xref>, an increase in the <italic>G</italic><sub><italic>p</italic></sub> value leads to a more pronounced sensitivity in GWP fluctuations.</p>
<fig id="fig-7">
<label>Figure 7</label>
<caption>
<title><italic>GWP</italic> variation with changes in <italic>G</italic><sub><italic>p</italic></sub></title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_43497-fig-7.tif"/>
</fig>
<p>In the carbon emission model of the wind turbine gearbox, adjusting the number of preventive maintenance sessions proves to be an effective strategy for reducing carbon emissions per unit of operation. By modulating the frequency of preventive maintenance, either by increasing or decreasing it, the gearbox&#x2019;s carbon footprint can be optimized. This approach facilitates a balance between operational efficiency and environmental impact, highlighting the importance of strategic maintenance planning in reducing ecological footprints.</p>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Conclusions</title>
<p>(1) This research demonstrates that compared to traditional TBM utilizing ODE, CBM employing SDE can more accurately track and simulate the state changes of wind power gearboxes, offering enhanced real-time accuracy.</p>
<p>(2) The carbon emissions model, when applied to both CBM with SDE-based and TBM with ODE-based approaches, reveals that carbon emissions under CBM are generally lower than those under TBM. However, over the long term, the carbon emissions from both methodologies are expected to converge, as TBM represents the expected outcome of CBM.</p>
<p>(3) The research findings indicate that the stochastic differential equation model is adept at incorporating the effects of daily inspections, maintenance, and environmental factors on the gearbox, thereby better ensuring the reliable operation of wind turbine equipment. The carbon emission models are capable of calculating the carbon emissions under both TBM and CBM. Utilizing these models for carbon emission assessments allows for timely adjustments to the preventive maintenance schedules of gearboxes, thereby reducing carbon emissions while ensuring their reliable operation.</p>
<p>(4) The results show that carbon emissions under CBM are generally lower than those under TBM. Sensitivity analysis of carbon emissions arising from energy consumption factors, such as <italic>G</italic><sub><italic>use</italic></sub> (carbon emissions per unit of energy consumed), <italic>G</italic><sub><italic>m</italic></sub> (carbon emissions from a maintenance activity), and <italic>G</italic><sub><italic>p</italic></sub> (carbon emissions from manufacturing a device), suggests that a decrease in the <italic>G</italic><sub><italic>use</italic></sub> value and an increase in the <italic>G</italic><sub><italic>p</italic></sub> value make the carbon emissions per unit of energy consumption in the gearbox more responsive to changes.</p>
</sec>
</body>
<back>
<ack>
<p>The authors would like to express our sincere appreciation to the anonymous referees for providing valuable suggestions and comments that have significantly contributed to the improvement of our manuscript.</p>
</ack>
<sec><title>Funding Statement</title>
<p>This research was supported by Basic Science Research Program through the National Natural Science Foundation of China (Grant No. 61867003) and Key Project of Science and Technology Research and Development Plan of China Railway Co., Ltd. (N2022X009).</p>
</sec>
<sec><title>Author Contributions</title>
<p>The authors confirm contribution to the paper as follows: study conception and design: Hongsheng Su, Lixia Dong, Xiaoying Yu; data collection: Lixia Dong, Kai Liu; analysis and interpretation of results: Lixia Dong, Hongsheng Su, Xiaoying Yu, Kai Liu; draft manuscript preparation: Hongsheng Su, Lixia Dong, Xiaoying Yu, Kai Liu. All authors reviewed the results and approved the final version of the manuscript.</p>
</sec>
<sec sec-type="data-availability"><title>Availability of Data and Materials</title>
<p>The authors confirm that the data supporting the findings of this study are available within the article.</p>
</sec>
<sec sec-type="COI-statement"><title>Conflicts of Interest</title>
<p>The authors declare that they have no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>1.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>M&#x00E1;rquez</surname>, <given-names>G. P. F.</given-names></string-name>, <string-name><surname>Tobias</surname>, <given-names>M. A.</given-names></string-name>, <string-name><surname>P&#x00E9;rez</surname>, <given-names>P. M. J.</given-names></string-name>, <string-name><surname>Papaelias</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2012</year>). <article-title>Condition monitoring of wind turbines: Techniques and methods</article-title>. <source>Renewable Energy</source><italic>,</italic> <volume>46</volume><italic>,</italic> <fpage>169</fpage>&#x2013;<lpage>178</lpage>.</mixed-citation></ref>
<ref id="ref-2"><label>2.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tchakoua</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Wamkeue</surname>, <given-names>R.</given-names></string-name>, <string-name><surname>Ouhrouche</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Hasnaoui</surname>, <given-names>S. F.</given-names></string-name>, <string-name><surname>Tameghe</surname>, <given-names>A. T.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2014</year>). <article-title>Wind turbine condition monitoring: State-of-the-art review, new trends, and future challenges</article-title>. <source>Energies</source><italic>,</italic> <volume>7</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>2595</fpage>&#x2013;<lpage>2630</lpage>.</mixed-citation></ref>
<ref id="ref-3"><label>3.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Song</surname>, <given-names>T. X.</given-names></string-name>, <string-name><surname>Tan</surname>, <given-names>T. Y.</given-names></string-name>, <string-name><surname>Han</surname>, <given-names>G. C.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Research on preventive maintenance strategies and systems for in-service ship equipment</article-title>. <source>Polish Maritime Research</source><italic>,</italic> <volume>29</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>85</fpage>&#x2013;<lpage>96</lpage>.</mixed-citation></ref>
<ref id="ref-4"><label>4.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Joel</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Kazem</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Keld</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Christopher</surname>, <given-names>D.</given-names></string-name></person-group> (<year>2015</year>). <article-title>Effect of preventive maintenance intervals on reliability and maintenance costs of wind turbine gearboxes</article-title>. <source>Wind Energy</source><italic>,</italic> <volume>18</volume><italic>(</italic><issue>11</issue><italic>),</italic> <fpage>2013</fpage>&#x2013;<lpage>2024</lpage>.</mixed-citation></ref>
<ref id="ref-5"><label>5.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chan</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Mo</surname>, <given-names>J.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Life cycle reliability and maintenance analyses of wind turbines</article-title>. <source>Energy Procedia</source><italic>,</italic> <volume>110</volume><italic>,</italic> <fpage>328</fpage>&#x2013;<lpage>333</lpage>.</mixed-citation></ref>
<ref id="ref-6"><label>6.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ling</surname>, <given-names>M. H.</given-names></string-name>, <string-name><surname>So</surname>, <given-names>H. Y.</given-names></string-name>, <string-name><surname>Balakrishnan</surname>, <given-names>N.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Likelihood inference under proportional hazards model for one-shot device testing</article-title>. <source>IEEE Transactions on Reliability</source><italic>,</italic> <volume>65</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>446</fpage>&#x2013;<lpage>458</lpage>.</mixed-citation></ref>
<ref id="ref-7"><label>7.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Kang</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Su</surname>, <given-names>H. S.</given-names></string-name></person-group> (<year>2015</year>). <article-title>Reliability analysis and control strategy design on preventive maintenance of repairable equipment</article-title>. <source>International Journal Control Automation and Systems</source><italic>,</italic> <volume>8</volume><italic>,</italic> <fpage>59</fpage>&#x2013;<lpage>80</lpage>.</mixed-citation></ref>
<ref id="ref-8"><label>8.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Su</surname>, <given-names>H. S.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Stochastic model analysis and control strategy design on preventive maintenance based on condition</article-title>. <source>International Journal Control Automation and Systems</source><italic>,</italic> <volume>9</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>194</fpage>&#x2013;<lpage>214</lpage>.</mixed-citation></ref>
<ref id="ref-9"><label>9.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Su</surname>, <given-names>H. S.</given-names></string-name></person-group> (<year>2014</year>). <article-title>Preventive maintenance model analysis based on condition</article-title>. <source>International Journal of Security Applications</source><italic>,</italic> <volume>8</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>353</fpage>&#x2013;<lpage>366</lpage>.</mixed-citation></ref>
<ref id="ref-10"><label>10.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Franciosi</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Lambiase</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Miranda</surname>, <given-names>S.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Sustainable maintenance: A periodic preventive maintenance model with sustainable spare parts management</article-title>. <source>IFAC-PapersOnLine</source>. <ext-link ext-link-type="uri" xlink:href="https://www.sciencedirect.com/science/article/pii/S2405896317334584">https://www.sciencedirect.com/science/article/pii/S2405896317334584</ext-link> <comment>(accessed on 15/11/2023)</comment>.</mixed-citation></ref>
<ref id="ref-11"><label>11.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zied</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Nidhal</surname>, <given-names>R.</given-names></string-name>, <string-name><surname>Ali</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Maintenance on leasing sales strategies for manufacturing/remanufacturing system with increasing failure rate and carbon emission</article-title>. <source>International Journal of Production Research</source><italic>,</italic> <volume>58</volume><italic>(</italic><issue>21</issue><italic>),</italic> <fpage>6616</fpage>&#x2013;<lpage>6637</lpage>.</mixed-citation></ref>
<ref id="ref-12"><label>12.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Afrinaldi</surname>, <given-names>F.</given-names></string-name>, <string-name><surname>Taufik.</surname></string-name>, <string-name><surname>Tasman</surname>, <given-names>M. A.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>H. C.</given-names></string-name>, <string-name><surname>Hasan</surname>, <given-names>A.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Minimizing economic and environmental impacts through an optimal preventive replacement schedule: Model and application</article-title>. <source>Journal of Cleaner Production</source><italic>,</italic> <volume>143</volume><italic>,</italic> <fpage>882</fpage>&#x2013;<lpage>893</lpage>.</mixed-citation></ref>
<ref id="ref-13"><label>13.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>Y. T.</given-names></string-name>, <string-name><surname>Liu</surname>, <given-names>Q. M.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>B. S.</given-names></string-name></person-group> (<year>2023</year>). <article-title>Research on optimization of preventive maintenance of equipment considering carbon emissions</article-title>. <source>Industrial Engineering and Management</source><italic>,</italic> <volume>28</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>43</fpage>&#x2013;<lpage>51</lpage> <comment>(In Chinese)</comment>.</mixed-citation></ref>
<ref id="ref-14"><label>14.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname>, <given-names>Q. M.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>Z. N.</given-names></string-name>, <string-name><surname>Xia</surname>, <given-names>T. B.</given-names></string-name>, <string-name><surname>Minchih</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>J. X.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Integrated structural dependence and stochastic dependence for opportunistic maintenance of wind turbines by considering carbon emissions</article-title>. <source>Energies</source><italic>,</italic> <volume>39</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>127</fpage>&#x2013;<lpage>136</lpage>.</mixed-citation></ref>
<ref id="ref-15"><label>15.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Su</surname>, <given-names>H. S.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>D. T.</given-names></string-name>, <string-name><surname>Duan</surname>, <given-names>X. P.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Condition maintenance decision of wind turbine gearbox based on stochastic differential equation</article-title>. <source>Energies</source><italic>,</italic> <volume>13</volume><italic>(</italic><issue>17</issue><italic>),</italic> <fpage>138</fpage>&#x2013;<lpage>159</lpage>.</mixed-citation></ref>
<ref id="ref-16"><label>16.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Su</surname>, <given-names>H. S.</given-names></string-name>, <string-name><surname>Huangfu</surname>, <given-names>L. L.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Preventive maintenance model analysis on wind-turbine gearbox under stochastic disturbance</article-title>. <source>Energy Reports</source><italic>,</italic> <volume>8</volume><italic>(</italic><issue>Supplement 1</issue><italic>),</italic> <fpage>224</fpage>&#x2013;<lpage>231</lpage>.</mixed-citation></ref>
<ref id="ref-17"><label>17.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Bi</surname>, <given-names>X. T.</given-names></string-name>, <string-name><surname>Yang</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Yang</surname>, <given-names>S. Y.</given-names></string-name></person-group> (<year>2021</year>). <article-title>LCA-based regional distribution and transference of carbon emissions from wind farms in China</article-title>. <source>Energies</source><italic>,</italic> <volume>15</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>198</fpage>.</mixed-citation></ref>
<ref id="ref-18"><label>18.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Silvia</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Samuel</surname>, <given-names>S.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Practical example of modification of a gearbox lubrication system</article-title>. <source>Lubricants</source><italic>,</italic> <volume>10</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>110</fpage>.</mixed-citation></ref>
<ref id="ref-19"><label>19.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Xu</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Yu</surname>, <given-names>J. H.</given-names></string-name>, <string-name><surname>Chen</surname>, <given-names>Y.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Carbon emission and energy consumption accounting analysis of wind farm operation</article-title>. <source>Journal of Dalian Polytechnic University</source><italic>,</italic> <volume>39</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>281</fpage>&#x2013;<lpage>288</lpage> <comment>(In Chinese)</comment>.</mixed-citation></ref>
<ref id="ref-20"><label>20.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Qian</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Ma</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Cross</surname>, <given-names>P.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Integrated data-driven model-based approach to condition monitoring of the wind turbine gearbox</article-title>. <source>IET Renewable Power Generation</source><italic>,</italic> <volume>11</volume><italic>(</italic><issue>9</issue><italic>),</italic> <fpage>1177</fpage>&#x2013;<lpage>1185</lpage>.</mixed-citation></ref>
</ref-list>
</back></article>