<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1347789</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1347789</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Decentralized control strategy of air-conditioning loads for primary frequency regulation based on environment information</article-title>
<alt-title alt-title-type="left-running-head">Jiang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2023.1347789">10.3389/fenrg.2023.1347789</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Baofeng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Zheng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xiong</surname>
<given-names>Wei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2525044/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Lingyu</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Pengyu</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2592697/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liao</surname>
<given-names>Siyang</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Dispatching and Control Center</institution>, <institution>Central China Branch of State Grid Corporation of China</institution>, <addr-line>Wuhan</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Electrical Engineering and Automation</institution>, <institution>Wuhan University</institution>, <addr-line>Wuhan</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2337497/overview">Zhengmao Li</ext-link>, Aalto University, Finland</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2274284/overview">Xiao Hu</ext-link>, Northeast Electric Power University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2244063/overview">Yuchen Fang</ext-link>, Tsinghua University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Pengyu Wang, <email>2023102070046@whu.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>15</day>
<month>01</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1347789</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>18</day>
<month>12</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Jiang, Yang, Xiong, Chen, Wang and Liao.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Jiang, Yang, Xiong, Chen, Wang and Liao</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The existing research on orderly participation of air conditioning load in frequency regulation needs a relatively complete communication architecture, but it cannot be achieved on a large scale in practice. This paper proposes a decentralized control method for the air-conditioning load control to participate in the power system primary frequency regulation, so that the air-conditioning loads can respond to system frequency deviation in an orderly manner without communication. Firstly, the triggering frequency is set to prevent the air-conditioning load from being over- or under-adjusted based on the frequency response interval, which is determined according to the outdoor temperature. Secondly, the response priority is arranged according to the holding time to ensure an orderly response of the air-conditioning load without communication. Finally, a simulation example of a thousand air-conditioning loads under different control strategy proves that the proposed method can achieve better frequency regulation effect.</p>
</abstract>
<kwd-group>
<kwd>decentralized control</kwd>
<kwd>primary frequency regulation</kwd>
<kwd>demand response</kwd>
<kwd>thermostatically controlled loads</kwd>
<kwd>environmental information</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Sustainable Energy Systems</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>Maintaining frequency stability is a prerequisite for ensuring the stable operation of the power system and electricity supply for users. The real-time power balance between the generator side and the user side plays a decisive role in the stability of the system frequency. The traditional frequency regulation means of power system mainly relies on the generator side to track load fluctuation and stabilize the accident disturbance (<xref ref-type="bibr" rid="B5">Kundur et al., 1994</xref>). However, as a large number of renewable energy sources (RESs) are connected to the system, the randomness of RESs power output increases the probability and amplitude of system frequency fluctuations, and greatly increases the frequency regulation capacity requirements of the power system (<xref ref-type="bibr" rid="B4">Joos et al., 2000</xref>; <xref ref-type="bibr" rid="B11">Obaid et al., 2019</xref>). On the other hand, with the increasing of RESs proportion and the decreasing of conventional units proportion, the primary and secondary frequency regulation capacity of the traditional units and the system inertia will decline correspondingly, which brings further challenges to the frequency stability of the power system (<xref ref-type="bibr" rid="B9">Magdy et al., 2021</xref>).</p>
<p>Demand response (DR) participating in power system frequency regulation can greatly improve the frequency regulation capability and dynamic characteristics of the power system, which has become a research hotspot in the field of frequency regulation. The air-conditioning load is huge, accounting for about one-third of the total summer electricity load of Spain (<xref ref-type="bibr" rid="B13">Science Daily, 2012</xref>). In China, Jiangsu&#x2019;s air-conditioning load power consumption accounts for about 45% of the total summer load power consumption (<xref ref-type="bibr" rid="B15">State Grid and State Grid Jiangsu Electric Power Co Ltd, 2019</xref>).</p>
<p>The basic idea of the air-conditioning load participating in the frequency regulation in the existing study is to make the air-conditioning load respond to the system frequency deviation automatically according to the local measurement of system frequency. In terms of the overall effect, the droop response of the load power of the entire air-conditioning cluster to the system frequency can be realized. However, the quantity of the air conditioning loads is huge. How to coordinate multiple air-conditioning loads within the cluster and orderly exert the regulating ability of different air-conditioning loads is the main problem to be solved (<xref ref-type="bibr" rid="B7">Li et al., 2020</xref>).</p>
<p>Previous studies mainly focus on the coordinated control of multiple air-conditioning loads based on full communication conditions. That is, the master station collects the status of all air-conditioning loads and then according to a certain sorting rule makes centralized decision to issue control commands to different air-conditioning loads when the frequency deviation occurs in the system. Centralized control uses the unified analysis and calculation of the control center to determine which air-conditioner loads participate in frequency response under different grid states, and then directly issues switch commands to each air-conditioning load. In (<xref ref-type="bibr" rid="B8">Lu and Zhang, 2013</xref>), a centralized load controller and its control algorithm are designed to generate a priority list and issue scheduling commands, which can optimize the response sequence of air-conditioning loads to obtain better continuous adjustment reserves. Some scholars found that using &#x201c;transition time&#x201d; instead of indoor temperature as the basis for priority ranking can improve the regulation effect of system frequency (<xref ref-type="bibr" rid="B20">Zhou et al., 2015</xref>). It needs to be pointed out that the centralized control method has high requirements for the completeness and timeliness of communication. However, it is actually difficult to establish complete and timely communication between the control center and each air-conditioning load under the existing conditions. Therefore, this type of control mode has great limitations in practical applications.</p>
<p>Another air-conditioning load control mode is &#x201c;centralized-distributed control,&#x201d; which adopts the basic mode of &#x201c;upward aggregation-centralized analysis-unified setting-decentralized response.&#x201d; The response mechanism is set for each air-conditioning load in advance, so that they can respond immediately when deviation occurs. In <xref ref-type="bibr" rid="B2">Biegel et al. (2013)</xref>, <xref ref-type="bibr" rid="B19">Zhang et al. (2013)</xref>, <xref ref-type="bibr" rid="B14">Song et al. (2019)</xref>, the control center sends the same control signal &#x3b1; to each air-conditioning load so that it will turn on/off according to the probability &#x3b1;, so as to achieve different adjustments as a whole. In <xref ref-type="bibr" rid="B1">Bao et al. (2015)</xref>, <xref ref-type="bibr" rid="B17">Wu et al. (2018)</xref>, different triggering frequencies are set for different air-conditioning loads through priority sorting, so that air-conditioning loads with higher priority have lower trigger frequencies, so as to achieve orderly control of air-conditioning loads. The centralized-distributed control presents good control characteristics. Because the judgment mechanism of whether the air-conditioning load participates in the response is set in the user&#x2019;s local area, the control mode can meet the frequency modulation time scale requirements. However, the centralized-distributed control still requires complete communication conditions and requires the control center to establish communication with each air-conditioning load to obtain its operating status. Such complete communication has not actually been established in the actual power grid, and its cost is extremely high and difficult to achieve. Under the condition of no communication, if each air-conditioning load can make an autonomous decision on its own control strategy based on its own state to realize global coordination at the cluster level, it is of great practical significance for the air-conditioning loads to participate in the system frequency modulation.</p>
<p>The decentralized control without communication has no upward aggregation and centralized analysis links. In most existing studies, the air-conditioning load adjusts the temperature setting or switch state by itself based on the monitored system frequency information. In <xref ref-type="bibr" rid="B18">Xu et al. (2011)</xref>, the set temperature of the thermostatically controlled load is changed according to the pre-set fixed proportional coefficient and the system frequency deviation. However, because the proportional coefficient is fixed, the response amount will fluctuate greatly under different air-conditioning load totals, increasing the uncertainty of frequency modulation effect. In <xref ref-type="bibr" rid="B10">Molina-Garc&#xed;a et al. (2011)</xref>, different triggering frequency and reconnecting frequency are set for different air-conditioning loads as the basis for participating in response and reconnect to system. Besides, the action delay is set to reduce measurement noise. However, because the triggering frequency and reconnecting frequency are both fixed and cannot be adjusted according to its own working status, the response time of the air-conditioning load is difficult to be guaranteed. The decentralized control does not rely on communication facilities, has low cost and high reliability, and is easy to implement. However, the uncertainty caused by the inability to unified dispatch and control makes it difficult for frequency regulation.</p>
<p>In this case that the triggering frequency is fixed, the response capacity has not been fully utilized and it is difficult to meet various environmental conditions. In response to the above problems, this paper proposes a decentralized air-conditioning demand response strategy by autonomous local decision-making based on local temperature information. The local temperature information and system frequency information measured by the air-conditioning load are calculated locally to set the appropriate triggering frequency. The main contributions are as follows:<list list-type="simple">
<list-item>
<p>1) Decentralized control strategy to optimize the frequency response interval of the air-conditioning cluster is proposed by monitoring the environmental temperature. According to the relationship between environmental temperature and total air-conditioning load, frequency response interval is adjusted. Thus, the over-adjustment or under-adjustment of the air-conditioning loads can be avoided in the incomplete communication condition.</p>
</list-item>
<list-item>
<p>2) The response priority of the air-conditioning load is sorted based on the &#x201c;holding time,&#x201d; which is defined according to the room temperature and the set temperature. Thus, the relatively orderly control of the air-conditioning loads can be realized to provide a longer frequency adjustment time without compromising the service life of the air-conditioning.</p>
</list-item>
</list>
</p>
<p>The rest is organized as follows: <xref ref-type="sec" rid="s2">Section 2</xref> introduces basement such as air-conditioning load droop control, frequency response action mechanism and equivalent model. <xref ref-type="sec" rid="s3">Section 3</xref> explains the air-conditioning triggering frequency setting mechanism. <xref ref-type="sec" rid="s4">Section 4</xref> introduces the working mechanism of the air-conditioning local controller. <xref ref-type="sec" rid="s5">Section 5</xref> is simulations.</p>
</sec>
<sec id="s2">
<title>2 Problem formulation for demand response</title>
<sec id="s2-1">
<title>2.1 Architecture</title>
<p>Air-conditioning load participating in primary frequency regulation usually adopts proportional control based on the principle of analog generator primary frequency regulation. Air-conditioning active power adjustment is proportionally controlled by the system frequency deviation calculated in Eq. <xref ref-type="disp-formula" rid="e1">1</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>ref</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>ref</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the nominal frequency, <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the system frequency. When the system frequency is low <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, the total air-conditioning load power should be lowered, and some air-conditioning loads should be switched off. When the system frequency is high (<inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:mfenced open="" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, the total air-conditioning load power should be adjusted up, and some air-conditioning loads should be switched on. Considering that the two situations are similar, for the convenience of description, this paper discusses scenarios where the system frequency is low and the demand response is a load reduction. In an ideal state, the response amount should be calculated in Eq. <xref ref-type="disp-formula" rid="e2">2</xref> and Eq. <xref ref-type="disp-formula" rid="e3">3</xref>:<disp-formula id="e2">
<mml:math id="m6">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>db</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mtext>air_con</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>db</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>db</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m7">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mtext>air</mml:mtext>
<mml:mtext>con</mml:mtext>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>db</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<inline-formula id="inf5">
<mml:math id="m8">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum frequency deviation allowed by the system. <inline-formula id="inf6">
<mml:math id="m9">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the dead zone for the primary frequency regulation. <inline-formula id="inf7">
<mml:math id="m10">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum response that the air-conditioning loads can provide. <inline-formula id="inf8">
<mml:math id="m11">
<mml:mrow>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>con</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is a proportional coefficient.</p>
<p>In the actual response process, the power of the fixed-frequency air-conditioning cannot be continuously adjusted. Thus, the continuous adjustment of the total power of the air-conditioning cluster can be achieved by switching large-scale air-conditioning loads.</p>
<p>For a single air-conditioning load, the triggering frequency is pre-set, and the system frequency is monitored at the same time. By comparing the triggering frequency and the system frequency deviation, it is determined whether and how the air-conditioning load will respond. The air-conditioning response mechanism is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, and the ordinate represents the system frequency deviation. When the system frequency drops to <inline-formula id="inf9">
<mml:math id="m12">
<mml:mrow>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>ref</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the air-conditioning load switches off; when the system frequency rises to <inline-formula id="inf10">
<mml:math id="m13">
<mml:mrow>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>ref</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the air-conditioning load switches on. Since this paper only discusses the case where the system frequency is low <inline-formula id="inf11">
<mml:math id="m14">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> mentioned in this paper refers to <inline-formula id="inf12">
<mml:math id="m15">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>. Besides, for the convenience of narrative, <inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> mentioned in this paper refers to <inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>. The core problem to be done in this paper is to make dynamic decisions based on the air-conditioning status and environmental information.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Air-conditioning load response mechanism.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g001.tif"/>
</fig>
<p>For air-conditioning clusters, different triggering frequencies are set for each air-conditioning load to achieve different responses under different frequency deviations and the correlation in Eq. <xref ref-type="disp-formula" rid="e2">2</xref>.</p>
</sec>
<sec id="s2-2">
<title>2.2 Equivalent model of air-conditioning</title>
<p>This paper uses the classic ETP model as the air-conditioning model. All details about ETP model can be found in <xref ref-type="bibr" rid="B16">Wen et al. (2021)</xref>.</p>
<p>Indoor temperature can be expressed by Eq. <xref ref-type="disp-formula" rid="e4">4</xref>.<disp-formula id="e4">
<mml:math id="m18">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi mathvariant="italic">PR</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi mathvariant="italic">PR</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:msup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2a;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2a;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <inline-formula id="inf15">
<mml:math id="m19">
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the rated power of the air-conditioning. <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the equivalent heat capacity. <inline-formula id="inf17">
<mml:math id="m21">
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the equivalent thermal resistance. <inline-formula id="inf18">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf19">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the indoor temperature and the outdoor temperature respectively. The time interval between time <inline-formula id="inf20">
<mml:math id="m24">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and time <inline-formula id="inf21">
<mml:math id="m25">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> is <inline-formula id="inf22">
<mml:math id="m26">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the air-conditioning energy efficiency ratio.</p>
<p>The indoor temperature fluctuates between the maximum allowed value <inline-formula id="inf23">
<mml:math id="m27">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and minimum allowed value <inline-formula id="inf24">
<mml:math id="m28">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. <inline-formula id="inf25">
<mml:math id="m29">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf26">
<mml:math id="m30">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are determined by the user set value <inline-formula id="inf27">
<mml:math id="m31">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>set</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and threshold <inline-formula id="inf28">
<mml:math id="m32">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.<disp-formula id="e5">
<mml:math id="m33">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>set</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>&#x2001;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>set</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m34">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>if&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3c;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>&#x2009;if&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mtext>others</mml:mtext>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m35">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>if&#x2009;</mml:mtext>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>if&#x2009;</mml:mtext>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Air-conditioning load switch status change.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g002.tif"/>
</fig>
</sec>
</sec>
<sec id="s3">
<title>3 Design of decentralized control for air-conditioning loads based on local temperature</title>
<sec id="s3-1">
<title>3.1 Frequency response interval setting based on environmental information</title>
<p>Compared with the centralized control mode with communication, one of the characteristics of the decentralized control is uncertainty (<xref ref-type="bibr" rid="B6">Li et al., 2021</xref>). For the same air-conditioning cluster, there must be load peaks and load valleys. If the triggering frequencies of all the air-conditioning loads are between the fixed interval <inline-formula id="inf29">
<mml:math id="m36">
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> under any circumstances, the number of air-conditioning loads which response is far different when the same frequency deviation occurs in different scenarios. Whether it is over-adjustment during peak air-conditioning load or under-adjustment during low air-conditioning load, it will adversely affect the response results, and the uncertainty of response will also bring inconvenience to the planning of other frequency adjustment facilities. Therefore, it is necessary to set different frequency response intervals for different load conditions, so that the response amount of the air-conditioning load is basically stable when the same frequency deviation occurs in different scenarios.</p>
<p>The main factor that affects residents&#x2019; switching on and off of air-conditioning is weather, which is prominently manifested in temperature (<xref ref-type="bibr" rid="B3">De Felice et al., 2013</xref>). For the analytical expression of the correlation between the total air-conditioning load and the outdoor temperature, many documents have been studied (<xref ref-type="bibr" rid="B12">Pourmousavi and Nehrir, 2014</xref>). Because it is not included in the research content of this paper, it will not be repeated.</p>
<p>For the same system frequency deviation, compared to the larger total power of the air-conditioning cluster, when the total power of the air-conditioning cluster is small, in order to keep the air-conditioning cluster a relatively stable response so that its participation in demand response effect is not affected, the air-conditioning must be triggered between a smaller frequency interval, that is, the frequency response span <inline-formula id="inf30">
<mml:math id="m37">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>interval</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is smaller.</p>
<p>The frequency response interval can be given:<disp-formula id="e8">
<mml:math id="m38">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
<disp-formula id="e9">
<mml:math id="m39">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="normal">k</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mrow>
<mml:mtext>air_con</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>all</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf31">
<mml:math id="m40">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the dead zone of response, and its function is to prevent the air-conditioning loads from malfunctioning and avoid noise. The relationship between the total power of the air-conditioning cluster and the setting of the frequency response interval reflected in Eqs <xref ref-type="disp-formula" rid="e8">8</xref>, <xref ref-type="disp-formula" rid="e9">9</xref> is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The frequency response interval setting.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g003.tif"/>
</fig>
<p>Combining <xref ref-type="fig" rid="F3">Figure 3</xref>, Eqs <xref ref-type="disp-formula" rid="e8">8</xref>, <xref ref-type="disp-formula" rid="e9">9</xref> and the correlation between the total air-conditioning load and the outdoor temperature, the corresponding relationship between outdoor temperature and frequency response interval can be obtained.</p>
</sec>
<sec id="s3-2">
<title>3.2 Response priority ranking based on hold time</title>
<p>In order to ensure users&#x2019; comfort, when the indoor temperature rises to the upper limit <inline-formula id="inf32">
<mml:math id="m41">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> allowed by the user, the air-conditioning must be switched on again for cooling. Since different air-conditionings are in different temperature states at the moment when they participate in demand response, the time to switch on after they participate in demand response is also different. This interval is defined as &#x201c;the holding time&#x201d;. As shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, when the air-conditioning is off (point A), the hold time is <inline-formula id="inf33">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>; when the air-conditioning is on (point B), the hold time is <inline-formula id="inf34">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. In addition, it is easy to know that when the room temperature is <inline-formula id="inf35">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the holding time is the longest, and it is defined as &#x201c;the maximum holding time&#x201d;, which is <inline-formula id="inf36">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The holding time of air-conditioning load.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g004.tif"/>
</fig>
<p>According to Eqs <xref ref-type="disp-formula" rid="e4">4</xref>&#x2013;<xref ref-type="disp-formula" rid="e7">7</xref>, the calculation formulas for holding time and maximum holding time can be obtained:<disp-formula id="e10">
<mml:math id="m46">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3c4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mtext>man</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>ciot</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<disp-formula id="e11">
<mml:math id="m47">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>mat</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>eut</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>According to <xref ref-type="fig" rid="F4">Figure 4</xref> and Eq. <xref ref-type="disp-formula" rid="e11">11</xref>, it can be seen that the farther the indoor temperature is from the upper limit set by the user at the time of demand response, the longer the hold time. For system frequency modulation, the longer the hold time, the longer the air conditioning participates in the frequency modulation, which has a positive effect on the frequency modulation effect. For the air-conditioning itself, the holding time is too short and frequent switching will affect its service life. Therefore, the air-conditioning with a long holding time should be turned off first, that is, the air-conditioning has a higher priority to participate in the response.</p>
<p>It needs to be pointed out that for the same air-conditioning, as the indoor temperature drops, the holding time gradually increases, but the relationship between the two is not linear. According to <xref ref-type="fig" rid="F2">Figure 2</xref>, it can be seen that as the indoor temperature rises, its increasing speed gradually slows down. For different air-conditionings, it can be seen from Eq. <xref ref-type="disp-formula" rid="e10">10</xref> that the holding time is not only related to the indoor temperature, but also related to the heat capacity and thermal resistance. Therefore, for different air-conditioning loads, even if the indoor temperature is the same, the holding time is different, such as <italic>A</italic> and <italic>A</italic>&#x2019;in <xref ref-type="fig" rid="F4">Figure 4</xref>. However, for demand response, the grid pays more attention to the holding time of the air-conditionings participation in the response, so the holding time is used here instead of temperature.</p>
</sec>
<sec id="s3-3">
<title>3.3 Autonomous decision of triggering frequency based on response priority</title>
<p>The higher the priority of the air-conditioning, the sooner it should participate in system frequency regulation, that is, the lower the triggering frequency <inline-formula id="inf37">
<mml:math id="m48">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> should be set. Combined with the aforementioned priority ranking rule of &#x201c;the longer the holding time, the higher the priority,&#x201d; it can be seen that the air conditioning triggering frequency setting with the longer holding time should be lower.</p>
<p>Considering that the decentralized control does not have communication function, the triggering frequency is calculated locally according to the environmental information and air-conditioning parameters. The calculation method is as Eq. <xref ref-type="disp-formula" rid="e12">12</xref> and the complete flow chart for triggering frequency calculation is as <xref ref-type="fig" rid="F5">Figure 5</xref>.<disp-formula id="e12">
<mml:math id="m49">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>where <inline-formula id="inf38">
<mml:math id="m50">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the triggering frequency of air-conditioner <inline-formula id="inf39">
<mml:math id="m51">
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> at time, <inline-formula id="inf40">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>max</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the maximum holding time of air conditioner <inline-formula id="inf41">
<mml:math id="m53">
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf42">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the holding time of air-conditioner <inline-formula id="inf43">
<mml:math id="m55">
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> at time <inline-formula id="inf44">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>The flow chart for triggering frequency&#x2019;s calculation.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g005.tif"/>
</fig>
<p>Combining Eqs <xref ref-type="disp-formula" rid="e10">10</xref>&#x2013;<xref ref-type="disp-formula" rid="e11">11</xref>, Eq. <xref ref-type="disp-formula" rid="e12">12</xref> can be transformed into Eq. <xref ref-type="disp-formula" rid="e13">13</xref>:<disp-formula id="e13">
<mml:math id="m57">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>th</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>dh</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>th</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mtext>th</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>In addition, in order to avoid frequent switching on and off causing damage to the service life of the air-conditionings, a lock time is set. When the hold time is less than the lock time, the air-conditioning is locked and the triggering frequency is not set, that is, it does not participate in demand response.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Local controller mechanism and method application</title>
<p>The air-conditioning load local controller has five modules: information acquisition module, calculation module, monitoring module, control module and locking unit, which are respectively responsible for multi-party information acquisition, triggering frequency and parameter calculation, system frequency monitoring, air-conditioning switch control and air-conditioning lock.</p>
<p>As shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, the working mode of the local controller for air-conditioning load is as follows:</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Air-conditioning load local control mechanism.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g006.tif"/>
</fig>
<p>The information acquisition module transmits real-time updated information to the calculation module. The computing module performs online calculations of the air conditioner&#x2019;s holding time and trigger frequency, and sends this information respectively to the locking unit and monitoring module. By comparing system frequency offset and trigger frequency, the monitoring module makes decisions on air conditioner switching instructions, while the locking unit determines whether the air conditioner participates in the demand response by comparing the hold time and lock time. The control module receives instructions from both the monitoring module and locking unit in order to switch the operating state of the air conditioner.</p>
<sec id="s4-1">
<title>4.1 Information acquisition module</title>
<p>1) Collection of environmental information: sampling the indoor temperature and outdoor temperature where the air-conditioning is located once every minute; 2) Collection of external input information: temperature set by the user, frequency response dead zone, frequency response interval, etc., can be updated annually/quarterly according to actual needs, with extremely low timeliness and communication requirements and various implementation methods; 3) Collection of setting information for air-conditioning manufacturers: temperature setting threshold, lock time, rated power, etc.</p>
</sec>
<sec id="s4-2">
<title>4.2 Calculation module</title>
<p>1) Fit the curve based on historical data and combine Eq. <xref ref-type="disp-formula" rid="e5">5</xref> to calculate the equivalent heat capacity and thermal resistance; 2) Accept the data information from the information acquisition module, and calculate the hold time online according to (12) (13); 3) Calculate the triggering frequency according to (14).</p>
</sec>
<sec id="s4-3">
<title>4.3 Monitoring module</title>
<p>1) Real-time monitoring of the system frequency; 2) Once the frequency deviation exceeds the triggering frequency, a switch command will be issued.</p>
</sec>
<sec id="s4-4">
<title>4.4 Control module</title>
<p>Accept commands issued by the calculation module to directly control the air-conditioning switch.</p>
</sec>
<sec id="s4-5">
<title>4.5 Locking unit</title>
<p>Accept the hold time transmitted by the calculation module, compare the hold time with the lock time, if the hold time is less than the lock time, lock the control module to prevent the air-conditioning from participating in demand response.</p>
</sec>
<sec id="s4-6">
<title>4.6 Method application</title>
<p>The proposed method theoretically enables a large number of decentralized air conditioners to participate in the primary frequency regulation of a real-world scenario, providing auxiliary frequency regulation services for the power system. Firstly, in simulating the participation of air conditioning in primary frequency regulation in real scenarios, a power system frequency response (SFR) model can be constructed to accurately reflect the system frequency response characteristics through model parameter identification method. Secondly, based on the method proposed in this paper, the local controller of the air conditioner will dynamically calculate the triggering frequency and holding time. Different triggering frequencies make different response sequences of the air conditioner. By monitoring the dynamic frequency deviation of the system, the air conditioner can realize local autonomous decision-making under the premise of ensuring the operation life and the comfort requirements of users, which means the air conditioner cluster changes response power according to the variation of the system frequency in an orderly manner, and realize the rapid recovery of the system frequency.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Case study</title>
<p>The simulation assumes that 1,000 air conditioners equipped with local controllers are in the on state at the current outdoor temperature, and the temperature settings of each air-conditioning are the same. To simulate the differentiated features of actual scenarios, the indoor temperature where the air-conditioning is located is regarded as uniformly distributed within the allowable room temperature range, and the rated power of the air-conditioning load, the equivalent heat capacity of the room, and the equivalent thermal resistance obey the normal distribution (<xref ref-type="bibr" rid="B18">Xu et al., 2011</xref>). The outdoor temperature is regarded as the same and constant due to the short frequency modulation time scale and the parameter settings for simulation are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Parameter settings.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Parameter</th>
<th align="center">Value</th>
<th align="center">unit</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<inline-formula id="inf45">
<mml:math id="m58">
<mml:mrow>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">1,000</td>
<td align="center">
<inline-formula id="inf46">
<mml:math id="m59">
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf47">
<mml:math id="m60">
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Normal <inline-formula id="inf48">
<mml:math id="m61">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5,0.2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf49">
<mml:math id="m62">
<mml:mrow>
<mml:mo>&#x2103;</mml:mo>
<mml:mo>/</mml:mo>
<mml:mtext>kW</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf50">
<mml:math id="m63">
<mml:mrow>
<mml:mi mathvariant="normal">C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Normal <inline-formula id="inf51">
<mml:math id="m64">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>0.1,0.02</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf52">
<mml:math id="m65">
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mtext>kWh</mml:mtext>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:mo>&#x2218;</mml:mo>
</mml:msup>
<mml:mi mathvariant="normal">C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf53">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Normal <inline-formula id="inf54">
<mml:math id="m67">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1.5,0.2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf55">
<mml:math id="m68">
<mml:mtext>kW</mml:mtext>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf56">
<mml:math id="m69">
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mo>&#x2005;</mml:mo>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf57">
<mml:math id="m70">
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>24,26</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf58">
<mml:math id="m71">
<mml:mrow>
<mml:mo>&#x2103;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf59">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mrow>
<mml:mtext>out</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">35</td>
<td align="center">
<inline-formula id="inf60">
<mml:math id="m73">
<mml:mrow>
<mml:mo>&#x2103;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf61">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mrow>
<mml:mtext>in</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">Uniform <inline-formula id="inf62">
<mml:math id="m75">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>24,26</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf63">
<mml:math id="m76">
<mml:mrow>
<mml:mo>&#x2103;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In a certain air-conditioning cluster, the smaller the triggering frequency value, the more priority the air-conditioning participates in the demand response. In other words, the triggering frequency represents the order in which the air-conditioning participates in the demand response. If the triggering frequency value is randomly and uniformly assigned to the air-conditioning loads without priority calculation and sorting, the room temperature of the air-conditioning has nothing to do with the triggering frequency value, and it is in a disordered state most of the time, such as <xref ref-type="fig" rid="F7">Figure 7B</xref>. In some special cases, there may also be extreme situations in which the room temperature decreases as the triggering frequency increases, such as <xref ref-type="fig" rid="F7">Figure 7A</xref>. If <xref ref-type="fig" rid="F7">Figures 7A, B</xref> is assigned to the triggering frequency through the priority arrangement of this paper, then <xref ref-type="fig" rid="F7">Figure 7C</xref> can be obtained. The three scenes of <xref ref-type="fig" rid="F7">Figures 7A&#x2013;C</xref> are called scene 1, scene 2 and scene 3 respectively, and simulations are performed in the three cases respectively.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Room temperature and triggering frequency distribution <bold>(A)</bold> The triggering frequency increases as room temperature decreases <bold>(B)</bold> A disordered state <bold>(C)</bold> The triggering frequency increases as room temperature increases (strategy in this paper).</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g007.tif"/>
</fig>
<p>
<statement content-type="case" id="CASE_1">
<label>CASE 1</label>
<p>Step disturbance</p>
<p>Assuming that the system suffers a step disturbance of 0.1 p.u. (account for <inline-formula id="inf64">
<mml:math id="m77">
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> of the total electric load), such as generator loss, the total load of the air-conditioning cluster is set to <inline-formula id="inf65">
<mml:math id="m78">
<mml:mrow>
<mml:mn>0.2</mml:mn>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.u. <xref ref-type="fig" rid="F8">Figure 8</xref> is the simulation results. Observing and comparing the simulation results can draw the following conclusions:</p>
<p>Compared with scene l and scene 2, scene 3 has been prioritized, and the lowest system frequency has been significantly improved, which is about <inline-formula id="inf66">
<mml:math id="m79">
<mml:mrow>
<mml:mn>0.04</mml:mn>
<mml:mo>&#x2005;</mml:mo>
<mml:mtext>Hz</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> higher than scene l. And it also behaves more smoothly in the subsequent frequency fluctuation phase. The RMSE of the frequency stabilization stage in the three scenarios are 0.0146, 0.0070, and 0.0029, and the MSE are 0.0049, 0.0039, and 0.0026, which respectively reflect the size of the system frequency deviation and the degree of system frequency fluctuation. Both of them show a decreasing state from scene l to scene 3, indicating that the overall system frequency deviation and the fluctuation range of scene 3 are both the smallest.</p>
<p>Each point in <xref ref-type="fig" rid="F9">Figure 9</xref> represents one of the air-conditioning loads that participated in the demand response. The horizontal and vertical coordinates are respectively the moment of its participation in the demand response and its holding time (p.u.). Due to the difference of the holding time of the air conditioners, the hold time is equivalent to the identification code of the air conditioner. That is, from left to right, it represents the air-conditioning load&#x2019;s participation in demand response sequence. <xref ref-type="fig" rid="F10">Figure 10</xref> represents the change of the probability distribution of the indoor temperature of the air-conditioning over time. This calculation can be obtained by the kernel density estimation of the ksdensity function in MATLAB according to the sample data. Analyzing <xref ref-type="fig" rid="F9">Figure 9</xref> and <xref ref-type="fig" rid="F10">Figure 10</xref>, we can draw conclusions:</p>
<p>In scene 1, air-conditioning loads basically follow the principle that the shorter the hold time is, the faster the response is, which can also be observed in <xref ref-type="fig" rid="F10">Figure 10</xref>. The probability of air-conditioning at indoor temperature close to <inline-formula id="inf67">
<mml:math id="m80">
<mml:mrow>
<mml:msup>
<mml:mn>26</mml:mn>
<mml:mo>&#x2218;</mml:mo>
</mml:msup>
<mml:mi mathvariant="normal">C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> first decreases and then increases. Therefore, because some air-conditioning loads at indoor temperature close to <inline-formula id="inf68">
<mml:math id="m81">
<mml:mrow>
<mml:msup>
<mml:mn>26</mml:mn>
<mml:mo>&#x2218;</mml:mo>
</mml:msup>
<mml:mi mathvariant="normal">C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> immediately reopened after responding, the frequency modulation effect was not good, and there is still a large amount of air-conditioning loads response after <inline-formula id="inf69">
<mml:math id="m82">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2005;</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>In scene 2, the air-conditioning loads participating in the demand response have a relatively uniform distribution of indoor temperature, and the probability distribution of the indoor temperature varies less. However, because there is still air-conditioning loads with short holding time participating in the response, the rapid load rebound is obvious, resulting in a poor frequency modulation effect.</p>
<p>In scene 3, due to the differences in parameters such as heat capacity and thermal resistance, the air-conditioning response sequence does not strictly follow the principle of &#x201c;air conditioning with a long holding time responds first,&#x201d; but the air-conditioning load that participates in demand response as a whole has a longer holding time, which obtains a longer frequency adjustment time and a better frequency modulation effect.</p>
</statement>
</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Simulation results subjected to step disturbance <bold>(A)</bold> System frequency fluctuation <bold>(B)</bold> Air-conditioning load power reduction.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Response order and hold time subjected to step disturbance <bold>(A)</bold> The triggering frequency increases as room temperature decreases <bold>(B)</bold> A disordered state <bold>(C)</bold> The triggering frequency increases as room temperature increases (strategy in this paper).</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Indoor temperature distribution changes with time subjected to step disturbance <bold>(A)</bold> The triggering frequency increases as room temperature decreases <bold>(B)</bold> A disordered state <bold>(C)</bold> The triggering frequency increases as room temperature increases (strategy in this paper).</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g010.tif"/>
</fig>
<p>
<statement content-type="case" id="CASE_2">
<label>CASE 2</label>
<p>Continuous disturbance</p>
<p>The system sometimes encounters continuous disturbances, such as the continuous sharp drop in wind power output, etc. Set it to <xref ref-type="fig" rid="F11">Figure 11</xref>. The total load of the air-conditioning cluster is still set to 0.2p.u. The simulations were performed in the three scenes of <xref ref-type="fig" rid="F7">Figure 7</xref>. Simulation results as shown in <xref ref-type="fig" rid="F12">Figure 12</xref>.</p>
<p>Observing <xref ref-type="fig" rid="F12">Figure 12</xref>, it can be found that the frequency deviation of scene 1 is the largest and the fluctuation is the most intense, scene 2 is the second, and the frequency deviation of scene 3 is the smallest and the fluctuation is small. In addition, the RMSE of the frequency stabilization stage in the three scenarios are 0.0173, 0.0074, and 0.0053, and the MSE are 0.0089, 0.0051, and 0.0043, which can also be used as evidence.</p>
<p>The meanings of <xref ref-type="fig" rid="F13">Figure 13</xref> and <xref ref-type="fig" rid="F14">Figure 14</xref> are similar to Case 1. The air-conditioning loads involved in demand response in Scene 1 has a shorter holding time, and there are still many air-conditioning loads responding around 10&#xa0;s, which directly leads to the subsequent frequency fluctuations.</p>
</statement>
</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>The continuous disturbance.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Simulation results subjected to continuous disturbance <bold>(A)</bold> System frequency fluctuation <bold>(B)</bold> Air-conditioning load power reduction.</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g012.tif"/>
</fig>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Response order and hold time subjected to continuous disturbance <bold>(A)</bold> The triggering frequency increases as room temperature decreases <bold>(B)</bold> A disordered state <bold>(C)</bold> The triggering frequency increases as room temperature increases (strategy in this paper).</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g013.tif"/>
</fig>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Indoor temperature distribution changes with time subjected to continuous disturbance <bold>(A)</bold> The triggering frequency increases as room temperature decreases <bold>(B)</bold> A disordered state <bold>(C)</bold> The triggering frequency increases as room temperature increases (strategy in this paper).</p>
</caption>
<graphic xlink:href="fenrg-11-1347789-g014.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>This paper proposes an adaptive control strategy for air-conditioning loads based on environmental perception. Decentralized control strategy for coordination of multiple air-conditioning loads is proposed based on the frequency response interval and the response priority in the condition that the communication is incomplete. On the one hand, by optimizing the setting of the frequency response interval, the over- or under-adjustment of the air-conditioning loads is avoided to obtain a better frequency modulation effect, and it will not cause excessive overshoot to bring new impact to the grid. On the other hand, the definition of holding time and prioritization realizes relatively orderly control of air-conditioning loads, which will not damage the service life of air-conditioning while gaining longer frequency adjustment time for the power grid. Based on the method proposed in this paper, the relative order control of air-conditioning cluster under incomplete communication can be realized. However, this paper mainly focuses on the orderly control strategy of decentralized air conditioning participating in primary frequency regulation, without considering the robustness of the control strategy under system uncertainties. We will conduct further research on this issue in the future.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>BJ: Funding acquisition, Project administration, Writing&#x2013;review and editing. ZY: Data curation, Investigation, Writing&#x2013;review and editing. WX: Conceptualization, Methodology, Writing&#x2013;review and editing. LC: Formal Analysis, Methodology, Writing&#x2013;original draft, Validation. PW: Methodology, Software, Visualization, Writing&#x2013;review and editing. SL: Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s9">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the technology project of Central China Branch of State Grid Corporation of China (SGHZ0000DKJS2310392). This paper is supported by the technology project: research on key technologies for inter provincial peak shaving and supply guarantee adapting to multi section constraints (Project No. 521400230001).</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>Authors BJ, ZY, and WX were employed by Central China Branch of State Grid Corporation of China.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Design of a hybrid hierarchical demand response control scheme for the frequency control</article-title>. <source>Generation Transm. Distribution IET</source> <volume>9</volume> (<issue>15</issue>), <fpage>2303</fpage>&#x2013;<lpage>2310</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2015.0628</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Biegel</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Hansen</surname>
<given-names>L. H.</given-names>
</name>
<name>
<surname>Andersen</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Stoustrup</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Primary control by ON/OFF demand-side devices</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>4</volume> (<issue>4</issue>), <fpage>2061</fpage>&#x2013;<lpage>2071</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2013.2257892</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De Felice</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Alessandri</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ruti</surname>
<given-names>P. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Electricity demand forecasting over Italy: potential benefits using numerical weather prediction models</article-title>. <source>Electr. Power Syst. Res.</source> <volume>104</volume>, <fpage>71</fpage>&#x2013;<lpage>79</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2013.06.004</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Joos</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ooi</surname>
<given-names>B. T.</given-names>
</name>
<name>
<surname>Mcgillis</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Galiana</surname>
<given-names>F. D.</given-names>
</name>
<name>
<surname>Marceau</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2000</year>). &#x201c;<article-title>The potential of distributed generation to provide ancillary services</article-title>,&#x201d; in <conf-name>Proceedings of the 2000 Power Engineering Society Summer Meeting (Cat. No.00CH37134)</conf-name>, <conf-loc>Seattle, WA, USA</conf-loc>, <conf-date>July 2000</conf-date>, <fpage>1762</fpage>&#x2013;<lpage>1767</lpage>. <pub-id pub-id-type="doi">10.1109/PESS.2000.868792</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kundur</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Balu</surname>
<given-names>N. J.</given-names>
</name>
<name>
<surname>Lauby</surname>
<given-names>M. G.</given-names>
</name>
</person-group> (<year>1994</year>). &#x201c;<article-title>Control of active power and reactive power</article-title>,&#x201d; in <source>Power system stability and control</source> (<publisher-loc>NY</publisher-loc>: <publisher-name>McGraw-Hill</publisher-name>), <fpage>159</fpage>&#x2013;<lpage>161</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Stochastic-weighted robust optimization based bilayer operation of a multi-energy building microgrid considering practical thermal loads and battery degradation</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>13</volume> (<issue>2</issue>), <fpage>668</fpage>&#x2013;<lpage>682</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2021.3126776</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Robust coordination of a hybrid AC/DC multienergy ship microgrid with flexible voyage and thermal loads</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>11</volume> (<issue>4</issue>), <fpage>2782</fpage>&#x2013;<lpage>2793</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2020.2964831</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Design considerations of a centralized load controller using thermostatically controlled appliances for continuous regulation reserves</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>4</volume> (<issue>2</issue>), <fpage>914</fpage>&#x2013;<lpage>921</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2012.2222944</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Magdy</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Bakeer</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Alhasheem</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Superconducting energy storage technology-based synthetic inertia system control to enhance frequency dynamic performance in microgrids with high renewable penetration</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>6</volume> (<issue>1</issue>), <fpage>36</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-021-00212-z</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Molina-Garc&#xed;a</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bouffard</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Kirschen</surname>
<given-names>D. S.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Decentralized demand-side contribution to primary frequency control</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>26</volume> (<issue>1</issue>), <fpage>411</fpage>&#x2013;<lpage>419</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2010.2048223</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Obaid</surname>
<given-names>Z. A.</given-names>
</name>
<name>
<surname>Cipcigan</surname>
<given-names>L. M.</given-names>
</name>
<name>
<surname>Abrahim</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Muhssin</surname>
<given-names>M. T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Frequency control of future power systems: reviewing and evaluating challenges and new control methods</article-title>. <source>J. Mod. Power Syst. Clean. Energy</source> <volume>7</volume>, <fpage>9</fpage>&#x2013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1007/s40565-018-0441-1</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pourmousavi</surname>
<given-names>S. A.</given-names>
</name>
<name>
<surname>Nehrir</surname>
<given-names>M. H.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Introducing dynamic demand response in the LFC model</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>29</volume> (<issue>4</issue>), <fpage>1562</fpage>&#x2013;<lpage>1572</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2013.2296696</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="web">
<collab>Science Daily</collab> (<year>2012</year>). <article-title>Air conditioning consumers one third of peak electric consumption in the summer</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.sciencedaily.com">https://www.sciencedaily.com</ext-link> (Accessed October 22, 2012)</comment>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Shahidehpour</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>State space modeling and control of aggregated TCLs for regulation services in power grids</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>10</volume> (<issue>4</issue>), <fpage>4095</fpage>&#x2013;<lpage>4106</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2018.2849321</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="web">
<collab>State Grid, State Grid Jiangsu Electric Power Co Ltd</collab> (<year>2019</year>). <article-title>The highest dispatching power load of Jiangsu Power Grid reaches a new peak</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://www.js.sgcc.com.cn">http://www.js.sgcc.com.cn</ext-link> (Accessed July 26, 2019)</comment>.</citation>
</ref>
<ref id="B16">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Wen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Q. H.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>A joint strategy of air-conditioning and electrolytic aluminum participating in primary frequency modulation</article-title>,&#x201d; in <conf-name>Proceedings of the 2021 IEEE 5th Conference on Energy Internet and Energy System Integration (EI2)</conf-name>, <conf-loc>Taiyuan, China</conf-loc>, <conf-date>October 2021</conf-date>, <fpage>1368</fpage>&#x2013;<lpage>1373</lpage>. <pub-id pub-id-type="doi">10.1109/EI252483.2021.9713545</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Hierarchical control of residential HVAC units for primary frequency regulation</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>9</volume> (<issue>4</issue>), <fpage>3844</fpage>&#x2013;<lpage>3856</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2017.2766880</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ostergaard</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Togeby</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Demand as frequency controlled reserve</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>26</volume> (<issue>3</issue>), <fpage>1062</fpage>&#x2013;<lpage>1071</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2010.2080293</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Lian</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Kalsi</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Aggregated modeling and control of air conditioning loads for demand response</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>28</volume> (<issue>4</issue>), <fpage>4655</fpage>&#x2013;<lpage>4664</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2013.2266121</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B. B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>X. Q.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Provision of supplementary load frequency control via aggregation of air conditioning loads</article-title>. <source>Energies</source> <volume>8</volume> (<issue>12</issue>), <fpage>14098</fpage>&#x2013;<lpage>14117</lpage>. <pub-id pub-id-type="doi">10.3390/en81212417</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>