一种用于计算多区域建筑房间温湿度的数学模型

1.本发明涉及一种用于计算多区域建筑房间温湿度的数学模型。


背景技术：

2.hvac系统是现代建筑中广泛使用的系统之一，占建筑物能耗的60％。如何在确保居住者热舒适的情况下，提高hvac系统的能效对建筑节能意义重大。而且计算hvac系统的能耗是hvac系统进行节能优化控制研究的基础。因此，建立一个能准确计算多区域建筑房间温湿度的数学模型是非常有意义的。
3.研究人员们为了建立建筑数学模型，利用物理领域的类比原理，用电阻-电容(rc)网络方法来表示建筑围护结构。rc网络方法能将建筑围护结构复杂的热传导过程简化为简单的数学计算问题。如王胜伟等利用rc网络方法建立了简化的建筑能耗模型，该模型能可靠地预测建筑的动态热工性能，并获得了令人满意的准确性和鲁棒性。如金旭东等利用rc网络方法将多区域建筑的围护结构划分为一层结构，建立了多区域建筑模型，该模型虽然能够准确计算各个房间的能耗，但是没有考虑到多区域建筑空气的耦合问题，会导致模型计算的室温与实际的室温存在一定的偏差。如雷蕾等利用rc网络方法将建筑围护结构划分为两层结构，建立了多区域建筑数学模型，该研究考虑了室内空气对流和上下楼层之间的换热，该模型虽然能够准确计算各个房间的温湿度，但是没有明确考虑通过门窗渗透进入室内的能量以及人体、照明和设备的散热，这会导致模型的计算值与实际值存在一定的偏差。因此，本发明在充分考虑建筑围护结构热传导过程、室内空气对流、上下楼层之间的换热、通过门窗渗透进入室内的能量、人体散热量、照明散热量、设备散热量的情况下，建立一个能准确计算多区域建筑房间温湿度的数学模型。


技术实现要素：

4.本发明所解决的技术问题是提供一种能准确计算多区域建筑房间温湿度的数学模型，为房间温湿度的控制模型开发设计提供理论依据。本发明根据大量试验数据，利用rc网络方法分析和描述围护结构的热传导过程，得到了一种满足研究和工程要求的数学模型。经过多次计算检验，该数学模型能准确计算多区域建筑中房间的温湿度。
5.为实现以上技术目的，本发明采用以下技术方案：
6.一种用于计算多区域建筑房间温湿度的数学模型，建模过程包括以下步骤：
7.步骤一：将建筑围护结构划分为两层结构；
8.步骤二：将空调房间划分为送风区、工作区、回风区；
9.步骤三：建立建筑围护结构各节点和内外表面的能量平衡方程；
10.步骤四：建立送风区、工作区、回风区的能量平衡方程和湿度平衡方程；
11.步骤五：对公式进行化简，并消除公式中的变量t
wa,out
,t
wa,in
,t
fl,out
和t
fl,in
；
12.步骤六：将单房间数学模型扩展到整个建筑的多区域数学模型。
13.进一步地，步骤一中将建筑围护结构划分为两层结构，每层结构的热容量和热阻
对应于每层结构质点上的热容和热阻，用c表示建筑围护结构的热容，用表示建筑围护结构的热阻，其中，建筑围护结构的外表面与第二节点之间的厚度是墙厚的1/4，建筑围护结构内两个节点之间的厚度是墙厚的1/2，建筑围护结构的第一个节点与内表面之间的厚度是墙厚的1/4。
14.更进一步地，利用rc网络方法将建筑围护结构划分为两层结构，所述rc网络方法是利用电阻和电容的网络来描述和分析建筑围护结构的热传导过程。
15.进一步地，步骤三中通过分析和描述分别得到墙体内表面、第二节点、第一节点和外表面的能量平衡方程：
[0016][0017][0018][0019][0020]
式中，λ
wa
、δ
wa
、ρ
wa
、c
wa
、a
wa
分别表示墙体的导热系数、厚度、密度、热容、面积；t表示温度，k表示换热系数，下标wa2、wa,out、wa1、wa,in分别表示墙体的第二节点、外表面、第一个节点、墙体的内表面；t
out
表示室外空气温度；t
inr
表示室内的回风温度；et
sunfu
表示被外墙吸收的太阳辐射能量；et
sunen
表示被内墙吸收的太阳辐射能量；et
inner
表示室内热源通过辐射传递到内墙表面的能量。
[0021]
进一步地，步骤三中通过分析和描述分别得到地板/天花板内表面、第二节点、第一节点和外表面的能量平衡方程：
[0022][0023][0024][0025]
[0026]
式中，λ
fl
、δ
fl
、ρ
fl
、c
fl
、a
fl
分别表示地板/天花板的导热系数、厚度、密度、热容、面积；t
fl,do
表示楼下的温度；t表示温度，k表示换热系数，下标fl2、fl,out、fl1、fl,in分别表示地板/天花板的第二节点、外表面、第一个节点、内表面。
[0027]
进一步地，步骤四中通过分析分别得到送风区、工作区、回风区的能量平衡方程：
[0028][0029][0030][0031]
式中，ρa表示空气密度；v
ins
表示送风区的体积；ca表示空气的热容；t
ins
表示送风区的空气温度；a
wa,s
表示墙体与送风区的接触面积；a
fl,s
表示地板/天花板与送风区的接触面积；g
sa1
表示送风空气质流量；t
sa1
表示送风温度；g
sa,sw
表示送风区向工作区输送的送风质流量；v
inw
表示工作区的体积；t
inw
表示工作区的空气温度；a
wa,w
表示墙体与工作区的接触面积；a
fl,w
表示地板/天花板与工作区的接触面积；k
win
表示窗的换热系数；a
win
表示窗的面积；k
door
表示门的换热系数；a
door
表示门的面积；g
sa,wr
表示工作区向回风区输送的空气质流量；t
inr
表示回风区的空气温度；v
adj,i
表示相邻房间的体积；r表示换气次数；t
adj,i
表示相邻房间的空气温度；l表示房间某朝向上的门窗缝隙长度；l表示每米门窗缝隙的基准渗风量；m表示门窗缝隙的渗风量综合修正系数；c
p
表示干空气的定压比热；v表示渗透空气的体积流量；n表示房间内最大班次人数；g表示每个人的平均散热量；φ表示考虑不同工作场所的性质，取φ＝1.0；μ1表示同时使用系数；μ2表示整流器散热系数；μ3表示安装系数；n
il
表示灯具安装功率；η1表示电机容量利用系数；η2表示电动机负荷系数；η3表示同时使用系数；η表示电动机效率；n
eq
表示电动机安装功率；v
inr
表示回风区的体积；a
wa,r
表示墙体与回风区的接触面积；a
fl,r
表示地板/天花板与回风区的接触面积；g
sa2
表示回风空气质流量；t
sa2
表示回风温度。
[0032]
进一步地，步骤四通过分析分别得到中送风区、工作区、回风区的湿度平衡方程：
[0033][0034][0035][0036]
式中，m
ins
表示送风区的空气湿度；m
sa
表示送风湿度；m
inw
表示工作区的空气湿度；madj,i
表示相邻房间的空气湿度；mr表示每个人的散湿量；m
inr
表示回风区的空气湿度。
[0037]
进一步地，步骤五中利用向前差分对公式进行化简：
[0038][0039]
更进一步地，简化后的公式整理成矩阵形式：
[0040]
x
′
＝ax+bi+cu
[0041]
式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为预测值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量。
[0042]
进一步地，步骤六中所述的多区域建筑数学模型为：
[0043][0044]
式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为计算值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量。
[0045]
本发明具有以下优点：
[0046]
(1)本发明涉及一种用于计算多区域建筑房间温湿度的数学模型。多区域建筑房间温湿度的变化跟房间内外环境和空调系统运行情况息息相关。为产品研发和实际工程应用提供更多的理论依据，本发明利用rc网络方法分析和描述围护结构的热传导过程，提出了一种用于计算多区域建筑房间温湿度的数学模型。本发明解决的技术问题是提供一种能准确计算房间温湿度的数学模型，为房间温湿度的控制模型开发设计提供理论依据。本发明在考虑空调系统运行和房间内外环境的情况下，利用rc网络方法分析和描述建筑围护结构的热传导过程，得到了一种满足研究和工程要求的数学模型。经过多次计算检验，该数学模型能准确计算多区域建筑中房间的温湿度。
[0047]
(2)本发明将建筑围护结构的热传导和房间空气能量交换等复杂问题转化为简单的数学问题。能够实时计算出房间温湿度随着房间内外环境和空调系统运行情况的实时变化值。填补了目前尚无分析和描述多区域建筑房间温湿度数学模型的空白，为准确计算hvac系统能耗和hvac系统节能优化控制提供了一种新的思路。
附图说明
[0048]
图1为利用rc网络方法对建筑围护结构进行划分的示意图；
[0049]
图2为空调房间分区示意图；
[0050]
图3为多区域建筑中某一房间在夏季时室内温湿度的计算值和实际值实时对比曲线图；
[0051]
图4为多区域建筑中某一房间在冬季时室内温湿度的计算值和实际值实时对比曲线图。
具体实施方式
[0052]
为了更好地理解本发明，现采用以下实施例加以说明，以下实施例属于本发明的保护范围，但不限制本发明的保护范围。
[0053]
一种用于计算多区域建筑房间温湿度的数学模型，建模过程包括以下步骤：
[0054]
步骤一：利用rc网络方法将建筑围护结构划分为两层结构；
[0055]
步骤一中rc网络方法是利用电阻和电容的网络来描述和分析建筑围护结构的热传导过程，将建筑围护结构划分为两层结构，每层结构的热容量和热阻对应于每层结构质点上的热容和热阻，用c表示建筑围护结构的热容，用表示建筑围护结构的热阻，其中，建筑围护结构的外表面与第二节点之间的厚度是墙厚的1/4，建筑围护结构内两个节点之间的厚度是墙厚的1/2，建筑围护结构的第一个节点与内表面之间的厚度是墙厚的1/4，利用rc网络方法对建筑围护结构进行划分的示意图如图1所示；
[0056]
步骤二：将空调房间划分为送风区、工作区、回风区；
[0057]
步骤二中将空调房间划分为送风区、工作区、回风区的示意图如图2所示；
[0058]
步骤三：建立建筑围护结构各节点和内外表面的能量平衡方程；
[0059]
步骤三中通过分析和描述分别得到墙体内表面、第二节点、第一节点和外表面的能量平衡方程：
[0060][0061][0062][0063][0064]
式中，λ
wa
、δ
wa
、ρ
wa
、c
wa
、a
wa
分别表示墙体的导热系数、厚度、密度、热容、面积；t表示温度，k表示换热系数，下标wa2、wa,out、wa1、wa,in分别表示墙体的第二节点、外表面、第一个节点、墙体的内表面；t
out
表示室外空气温度；t
inr
表示室内的回风温度；et
sunfu
表示被外墙吸收的太阳辐射能量；et
sunen
表示被内墙吸收的太阳辐射能量；et
inner
表示室内热源通过辐射传递到内墙表面的能量；
[0065]
步骤三中通过分析和描述分别得到地板/天花板内表面、第二节点、第一节点和外表面的能量平衡方程：
[0066]
[0067][0068][0069][0070]
式中，λ
fl
、δ
fl
、ρ
fl
、c
fl
、a
fl
分别表示地板/天花板的导热系数、厚度、密度、热容、面积；t
fl,do
表示楼下的温度；t表示温度，k表示换热系数，下标fl2、fl,out、fl1、fl,in分别表示地板/天花板的第二节点、外表面、第一个节点、内表面；
[0071]
步骤四：建立送风区、工作区、回风区的能量平衡方程和湿度平衡方程；
[0072]
步骤四中通过分析分别得到送风区、工作区、回风区的能量平衡方程：
[0073][0074][0075][0076]
式中，ρa表示空气密度；v
ins
表示送风区的体积；ca表示空气的热容；t
ins
表示送风区的空气温度；a
wa,s
表示墙体与送风区的接触面积；a
fl,s
表示地板/天花板与送风区的接触面积；g
sa1
表示送风空气质流量；t
sa1
表示送风温度；g
sa,sw
表示送风区向工作区输送的送风质流量；v
inw
表示工作区的体积；t
inw
表示工作区的空气温度；a
wa,w
表示墙体与工作区的接触面积；a
fl,w
表示地板/天花板与工作区的接触面积；k
win
表示窗的换热系数；a
win
表示窗的面积；k
door
表示门的换热系数；a
door
表示门的面积；g
sa,wr
表示工作区向回风区输送的空气质流量；t
inr
表示回风区的空气温度；v
adj,i
表示相邻房间的体积；r表示换气次数；t
adj,i
表示相邻房间的空气温度；l表示房间某朝向上的门窗缝隙长度；l表示每米门窗缝隙的基准渗风量；m表示门窗缝隙的渗风量综合修正系数；c
p
表示干空气的定压比热；v表示渗透空气的体积流量；n表示房间内最大班次人数；g表示每个人的平均散热量；φ表示考虑不同工作场所的性质，取φ＝1.0；μ1表示同时使用系数；μ2表示整流器散热系数；μ3表示安装系数；n
il
表示灯具安装功率；η1表示电机容量利用系数；η2表示电动机负荷系数；η3表示同时使用系数；η表示电动机效率；n
eq
表示电动机安装功率；v
inr
表示回风区的体积；a
wa,r
表示墙体与回风区的接
触面积；a
fl,r
表示地板/天花板与回风区的接触面积；g
sa2
表示回风空气质流量；t
sa2
表示回风温度；
[0077]
步骤四中通过分析分别得到送风区、工作区、回风区的湿度平衡方程：
[0078][0079][0080][0081]
式中，m
ins
表示送风区的空气湿度；m
sa
表示送风湿度；m
inw
表示工作区的空气湿度；m
adj,i
表示相邻房间的空气湿度；mr表示每个人的散湿量；m
inr
表示回风区的空气湿度；
[0082]
步骤五：对公式进行化简，并消除公式中的变量t
wa,out
,t
wa,in
,t
fl,out
和t
fl,in
；
[0083]
步骤五中利用向前差分对公式进行化简：
[0084][0085]
消除公式中的变量t
wa,out
,t
wa,in
,t
fl,out
和t
fl,in
，得：
[0086][0087][0088]
[0089][0090][0091]
[0092][0093][0094][0095][0096]
将简化后的公式整理成矩阵形式：
[0097]
x
′
＝ax+bi+cu
[0098]
式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为预测值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量；
[0099]
x'、x、i、u的表达式如下所示：
[0100][0101]
a的表达式如下所示：
[0102]
[0103][0104]
其中：
[0105]
[0106][0107][0108][0109][0110]
b的表达式如下所示：
[0111]
[0112][0113]
c的表达式如下所示：
[0114][0115]
其中：
[0116][0117]
[0118][0119]
步骤六：将单房间数学模型扩展到整个建筑的多区域数学模型；
[0120]
步骤六中所述的多区域建筑数学模型为：
[0121][0122]
式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为计算值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量。
[0123]
为了验证该数学模型能够准确计算多区域建筑中房间的温湿度，对桂林某大型商业建筑办公区中某个房间在冬夏两季的温湿度进行计算。建筑围护结构、室内外环境的参数如表1所示。房间温度和相对湿度采用si7021-a20-im数字相对湿度和温度传感器进行测量，其中，湿度比是从相对湿度和温度传感器的测量值转换而来。并把该传感器安装在房间中间距离地面1.5米处，每分钟记录一次房间的温度和相对湿度，该测量持续了2400s。在采集房间的温度和相对湿度的同时，室外温度和相对湿度也采用si7021-a20-im数字相对湿度和温度传感器进行测量。室外辐射强度采用sm206太阳能电表测量。
[0124]
表1建筑围护结构和室内外环境的参数
[0125]
[0126][0127]
图3和图4分别为多区域建筑中某个房间在夏季和冬季时室内温湿度的计算值和实际值实时对比曲线。
[0128]
从图3可以看出，房间在0-1200s期间，温度、湿度比的计算值和实测值变化曲线急剧下降，而且温度和湿度比的计算值比实测值略大，但两者相差均小于0.5℃。在1200s-2400s期间，温度、湿度比的计算值和实测值变化曲线趋于平缓，并逐渐重合。从图4可以看出，在0-1200s期间，房间的实测温度和湿度比的变化曲线始终在计算温度和湿度比的变化曲线上方。房间的温度和湿度比在1200s-2400s期间达到了稳定状态，两者的变化曲线逐渐重合。这表明，本发明建立的数学模型能够准确地、实时地计算出房间温度和湿度比。本发明提供了一种能准确计算多区域建筑房间温湿度的数学模型，为房间温湿度的控制模型开发设计提供理论依据。
[0129]
尽管已经描述和叙述了被看作本发明的示范实施例，本领域技术人员将会明白，可以对其作出各种改变和替换，而不会脱离本发明的精神。另外，可以做出许多修改以将特定情况适配到本发明的教义，而不会脱离在此描述的本发明中心概念。所以，本发明不受限于在此披露的特定实施例，但本发明可能还包括属于本发明范围的所有实施例及其等同物。

技术特征：
1.一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，建模过程包括以下步骤：步骤一：将建筑围护结构划分为两层结构；步骤二：将空调房间划分为送风区、工作区、回风区；步骤三：建立建筑围护结构各节点和内外表面的能量平衡方程；步骤四：建立送风区、工作区、回风区的能量平衡方程和湿度平衡方程；步骤五：对公式进行化简，并消除公式中的变量t
wa,out
,t
wa,in
,t
fl,out
和t
fl,in
；步骤六：将单房间数学模型扩展到整个建筑的多区域数学模型。2.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤一中将建筑围护结构划分为两层结构，每层结构的热容量和热阻对应于每层结构质点上的热容和热阻，用c表示建筑围护结构的热容，用表示建筑围护结构的热阻，其中，建筑围护结构的外表面与第二节点之间的厚度是墙厚的1/4，建筑围护结构内两个节点之间的厚度是墙厚的1/2，建筑围护结构的第一个节点与内表面之间的厚度是墙厚的1/4。3.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，利用rc网络方法将建筑围护结构划分为两层结构，所述rc网络方法是利用电阻和电容的网络来描述和分析建筑围护结构的热传导过程。4.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤三中通过分析和描述分别得到墙体内表面、第二节点、第一节点和外表面的能量平衡方程：衡方程：衡方程：衡方程：式中，λ
wa
、δ
wa
、ρ
wa
、c
wa
、a
wa
分别表示墙体的导热系数、厚度、密度、热容、面积；t表示温度，k表示换热系数，下标wa2、wa,out、wa1、wa,in分别表示墙体的第二节点、外表面、第一个节点、墙体的内表面；t
out
表示室外空气温度；t
inr
表示室内的回风温度；et
sunfu
表示被外墙吸收的太阳辐射能量；et
sunen
表示被内墙吸收的太阳辐射能量；et
inner
表示室内热源通过辐射传递到内墙表面的能量。5.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤三中通过分析和描述分别得到地板/天花板内表面、第二节点、第一节点和外表面的能量平衡方程：
式中，λ
fl
、δ
fl
、ρ
fl
、c
fl
、a
fl
分别表示地板/天花板的导热系数、厚度、密度、热容、面积；t
fl,do
表示楼下的温度；t表示温度，k表示换热系数，下标fl2、fl,out、fl1、fl,in分别表示地板/天花板的第二节点、外表面、第一个节点、内表面。6.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤四中通过分析分别得到送风区、工作区、回风区的能量平衡方程：于，步骤四中通过分析分别得到送风区、工作区、回风区的能量平衡方程：于，步骤四中通过分析分别得到送风区、工作区、回风区的能量平衡方程：式中，ρ
a
表示空气密度；v
ins
表示送风区的体积；c
a
表示空气的热容；t
ins
表示送风区的空气温度；a
wa,s
表示墙体与送风区的接触面积；a
fl,s
表示地板/天花板与送风区的接触面积；g
sa1
表示送风空气质流量；t
sa1
表示送风温度；g
sa,sw
表示送风区向工作区输送的送风质流量；v
inw
表示工作区的体积；t
inw
表示工作区的空气温度；a
wa,w
表示墙体与工作区的接触面积；a
fl,w
表示地板/天花板与工作区的接触面积；k
win
表示窗的换热系数；a
win
表示窗的面积；k
door
表示门的换热系数；a
door
表示门的面积；g
sa,wr
表示工作区向回风区输送的空气质流量；t
inr
表示回风区的空气温度；v
adj,i
表示相邻房间的体积；r表示换气次数；t
adj,i
表示相邻房间的空气温度；l表示房间某朝向上的门窗缝隙长度；l表示每米门窗缝隙的基准渗风量；m表示门窗缝隙的渗风量综合修正系数；c
p
表示干空气的定压比热；v表示渗透空气的体积流量；n表示房间内最大班次人数；g表示每个人的平均散热量；φ表示考虑不同工作场所的性质，
取φ＝1.0；μ1表示同时使用系数；μ2表示整流器散热系数；μ3表示安装系数；n
il
表示灯具安装功率；η1表示电机容量利用系数；η2表示电动机负荷系数；η3表示同时使用系数；η表示电动机效率；n
eq
表示电动机安装功率；v
inr
表示回风区的体积；a
wa,r
表示墙体与回风区的接触面积；a
fl,r
表示地板/天花板与回风区的接触面积；g
sa2
表示回风空气质流量；t
sa2
表示回风温度。7.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤四中通过分析分别得到送风区、工作区、回风区的湿度平衡方程：于，步骤四中通过分析分别得到送风区、工作区、回风区的湿度平衡方程：于，步骤四中通过分析分别得到送风区、工作区、回风区的湿度平衡方程：式中，m
ins
表示送风区的空气湿度；m
sa
表示送风湿度；m
inw
表示工作区的空气湿度；m
adj,i
表示相邻房间的空气湿度；m
r
表示每个人的散湿量；m
inr
表示回风区的空气湿度。8.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤五中利用向前差分对公式进行化简：9.根据权利要求8所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，简化后的公式整理成矩阵形式：x
′
＝ax+bi+cu式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为预测值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量。10.根据权利要求1所述的一种用于计算多区域建筑房间温湿度的数学模型，其特征在于，步骤六中所述的多区域建筑数学模型为：式中，a、b、c分别为房间的系统矩阵、输入矩阵和扰动矩阵，x'为计算值的向量，x为实际值的向量，i为输入值的向量，u为扰动量的向量。

技术总结
本发明公开了一种用于计算多区域建筑房间温湿度的数学模型，其建模过程包括以下步骤：(1)将建筑围护结构划分为两层结构；(2)将空调房间划分为送风区、工作区、回风区；(3)建立建筑围护结构各节点和内外表面的能量平衡方程；(4)建立送风区、工作区、回风区的能量平衡方程和湿度平衡方程；(5)对公式进行化简；(6)将单房间数学模型扩展到整个建筑的多区域数学模型。本发明在考虑空调系统运行和房间内外环境的情况下，利用RC网络方法分析和描述建筑围护结构的热传导过程，得到了一种满足研究和工程要求的数学模型，经过多次计算检验，本发明的数学模型能准确计算多区域建筑中房间的温湿度。的温湿度。的温湿度。


技术研发人员：雷蕾 梁利霞 胡佳敏
受保护的技术使用者：浙江理工大学
技术研发日：2023.03.27
技术公布日：2023/7/25