引言

气溶胶与云、雾、降水的相互耦合作用一直是大气科学领域的核心问题之一，而起、放电过程作为雷暴云动力-微物理过程的产物之一，其与气溶胶的关系日益受到重视。可溶气溶胶液滴的同质核化和冰核异质核化是冰晶粒子形成的两种来源，而雷暴云中冰晶粒子微物理特征对电荷的产生及其分层状况有一定影响，是导致闪电发生的重要因素之一。

可溶性气溶胶液滴的同质核化或是在大气冰核帮助下完成的异质核化，是云中冰晶粒子形成的重要来源(Pruppacher and Klett, 1997；Koop et al，2000；Liu and Penner, 2005；DeMott et al，2010)，因此，不同核化过程可能通过改变冰晶粒子特征影响雷暴电过程。由于目前对云微物理发展过程的观测手段比较有限，人们主要基于云模式进行冰晶核化过程对雷暴云电过程影响的研究，而且主要集中在研究大气冰核对雷暴动力、微物理及起电过程的作用上。例如，李泽宇等(2016)基于一维半云和气溶胶分档云模式，对不同初始沙尘浓度的非感应起电过程进行了敏感性试验，模拟结果表明：随着沙尘粒子数浓度的增多，云中冰晶粒子与霰粒子的数浓度都有所增加，初始起电时间提前，空间电荷密度大小增加。苏正军等(2010)基于催化剂检测试验的时间特征数据，对核化速率进行了深入研究，发现同一样品在不同温度下的动态特性呈现出显著差异：在低温区间(低于-16℃)，冰晶的形成是一个快速过程，其主要特征包括凝结冻结核化与凝华核化；相比之下，在较高温度区间(不低于-12℃)，冰晶的形成则主要由慢速过程主导，这一过程涉及接触核化与浸入核化机制。李璐滢等(2021)基于二维雷暴云模式，引入一种新型冰晶异质核化方案，揭示了不同的异质核化过程对雷暴电过程的影响，结果表明：浸润核化是冰晶生成的最重要异质核化过程，较高数浓度的冰晶消耗雷暴云内液态水含量，抑制凇附过程，导致霰粒子比含水量低，表现为较强的负极性非感应起电率；接触核化生成的冰晶量最少，仅对雷暴云中下层3~5 km处的冰晶有贡献，同时霰粒子数浓度较低，导致该方案下的起电过程最弱；沉积核化主要影响云砧处的冰晶，有利于提高霰收集云滴的效率，表现为极高的霰粒子比含水量，促进低温区非感应起电过程的发生。师正等(2023)研究了大气冰核浓度对雷暴云动力、微物理及空间电荷结构的影响，结果表明：大气冰核浓度提升有利于异质核化过程增强，冰晶在高温区大量生成，而同质核化过程被抑制，因此冰晶整体含量降低，引起低温区中霰粒子含量降低和高温区中霰粒子尺度降低。不同大气冰核浓度下的雷暴云空间电荷结构在雷暴云发展初期均呈现负的偶极性电荷结构。在雷暴云旺盛期，随着冰核浓度增加，空间电荷结构由三极性转变为复杂四极性。在雷暴云消散阶段，电荷结构均呈现偶极性，且电荷密度值随着冰核浓度的增加逐渐减小。

此外，科技进步促使多种观测技术被应用于雷暴活动过程的深入研究中。王婷波等(2020)依据国家闪电监测网络提供的云地闪电定位数据，对我国中东部地区的雷暴发生及其发展特性进行了统计分析。张萍萍等(2021)借助FY-4A闪电成像仪，有效地捕捉到了对流系统的发展演变特征。阮悦等(2022)在研究中整合了VLF/LF三维闪电监测数据、S波段双偏振雷达观测以及地面观测资料，揭示出在冰雹云的成熟阶段，云闪频率达到最高，而在其发展阶段则相对较低。大量的实际观测结果以及模式研究工作揭示雷暴云闪电活动特征与空间电荷结构密切相关。例如，通过对日本冬季雷暴的观测，Takeuti et al(1978)和Brook et al(1982)发现正地闪始发时，空间电荷结构均为偶极型。Carey and Rutledge(1998)发现云底部分出现足够大的负电荷堆有利于正地闪的形成。而冰雹云之所以能够生成大量的正地闪，可能与云中的反三极电荷结构有关(张义军等，2006)；对于负地闪而言，底部的次正电荷中心则是个关键(Clarence and Malan, 1957；Williams，1989；Tan et al，2014a)。底部次正电荷区足够大时，将容易产生反极性云闪(张义军等，2002；郄秀书等，2005；Tan et al，2014a)。闪电发生类别及地闪接地点的选择都受雷暴云电环境特征影响(Tan et al，2012；Tan et al，2014b)。雷暴云内不同核化过程对起电过程和空间电荷结构产生重要影响，而空间电荷结构又与放电过程息息相关，因此不同核化过程可能对雷暴云闪电活动特征产生重要影响。但目前仅有Gonçalves et al(2012)基于区域模式研究了细菌性冰核浓度对雷暴云微物理、降水及闪电发生率的作用，发现细菌性冰核会使对流云中的冰晶离子浓度增加，从而使闪电频率增加，总降水量减少。但其放电结果仅基于简单计算的经验公式，并未揭示其中详细的微物理和起电过程。

综上所述，冰晶核化过程对雷暴云起电过程存在影响，但不同核化过程和放电过程关系的研究很少。因此，为了全面了解冰晶核化过程影响雷暴云放电过程的物理原因，有必要开展相应的模式研究。本研究拟在已有的雷暴云三维随机起、放电模式基础上，将同质核化与异质核化参数化方案植入模式中，探讨不同核化方式对雷暴云放电特征的影响，并给出合理性解释。

1 模式介绍及初始场设置

本文所采用的云模式是中国气象科学研究院开发的云模式框架(胡志晋和何观芳，1987；王谦和胡志晋，1990；于达维等，2001)。在此基础上，谭涌波(2006)对Mansell et al(2002)提出的随机放电参数化方案进行了改进，建立了包含雷暴云起、放电参数化方案的三维雷暴云-电耦合模式。该模型视大气为无黏性的可压缩流体，采用非静力平衡的假设。模型涵盖云滴、雨滴、冰晶、霰和雹共5种水成物粒子，包含了22个预测变量和5个诊断变量，其涉及水成物的物理特性、空间总电荷密度、电位及电场强度等电过程的关键参数。下文将对模式采用的主要参数化方案进行简要介绍。

1.1 冰晶核化方案

迄今为止，雷暴云起、放电的模式中往往采用相对简单的经验公式来计算冰晶核化过程。例如，计算冰晶的核化率时往往采用一种与温度或相对于水的过饱和度相关的经验公式(Fletcher 1962；Cooper 1980；Meyers et al，1992)，其具有计算方便的优势，但由于不考虑冰核类型和具体核化物理过程，无法表征真实自然界中大气冰核浓度、时空分布和类型上的变化。想要合理地再现雷暴云起、放电过程必须以正确描述冰晶微物理特征为基础，因此，在雷暴云起、放电模式中完善冰晶核化过程十分重要。

本文引入了Koop and Murray(2016)改进的同质核化方案，通过研究特定温度下界面自由能的变化，确定相应的同质核化速率。直径为10 μm的液滴在-35℃的环境温度下需经过数小时才能冻结，而当低至-38℃后冻结仅需数秒时间，因此-38℃被认定为同质核化阈值(Liu et al，2005)。但环境温度低于阈值时，同质核化并不是必然发生，因为临界分子簇的形成具有随机性，冻结概率的高低取决于液滴体积和冻结持续时间(Knopf et al，2020)。模式中同质核化的计算基于经典成核理论，过冷液滴的同质核化速率取决于特定温度下界面自由能的变化。同时，为减少误差，低于-45℃时，所有液滴被处理为完全冻结，起始冻结温度定为-32℃。具体参数化方案见师正等(2023)。

本文在模式中植入了Diehl and Grützun(2018)提出的冰晶异质核化方案，将异质核分为接触核、浸润核和沉积核三种，综合考虑了各种不溶性气溶胶的核化特性，其中接触冻结核化通常是指干燥的气溶胶颗粒与液滴碰并冻结的过程。方案重点考虑：(1)环境温度低于冻结温度；(2)冰核与过冷液滴相互碰并；(3)作为冰核的气溶胶粒子为干燥状态。浸润冻结核化的温度与液滴体积存在正比的关系。根据随机假设，不可溶气溶胶表面成核时更低的能量势垒加速了液滴内分子簇的形成，提高了核化率，更大体积的液滴能捕获浸润核的概率也会增加。沉积核化与冰面过饱和度密切相关，冰面过饱和度的增加可有效提高冰核活化率。三种异质核化方案详见李璐滢等(2021)。

1.2 电过程参数化方案

模式中综合考虑了非感应起电过程和感应起电过程。霰(雹)与冰晶(雪晶)碰并回弹产生的非感应起电参数化方案采用GZ(Gardiner/Ziegler)方案(Gardiner et al，1985；Ziegler et al，1991；Shi et al，2015)；感应起电过程考虑霰(雹)和云滴之间的碰撞分离，其方案基于Ziegler et al(1991)的方案。

本文采用的随机放电参数化方案基于谭涌波(2006)的方案。当空间电场值达到逃逸击穿区域内电场阈值(随高度变化)或模式中任何位置的电场强度超过常规的空气击穿阈值(160 kV·m^-1)时，闪电启动；通道的随机演变过程采用了Mansell et al(2005)提出的方法，这种方法强调闪电通道的逐步发展，即每次正、负通道各自仅扩展至一个后继的通道点。地闪发生时，通道会把地面电位引到云内，从而导致强烈的后续放电过程。为了更加合理地模拟地闪放电过程以及云内后续放电，在模式中限制了后续放电的次数(小于5次)。为了限制雷暴云中电荷的过量累积，假定每次放电结束后空间电荷减小至原来的70%。

1.3 模式启动条件

在对流云数值模式中，通常采用两种方法在水平均匀的环境场中触发对流过程。一种方法适用于由动力抬升机制引发的对流过程，称为冷空气出流辐合启动机制；另一种为热力扰动启动机制。本研究采取后者，并在其基础上考虑湿度扰动因素，在模式域的下层施加一个范围有限的湿热扰动场，其湿度和温度均高于周边环境，并利用垂直运动方程中的浮力项来初始化对流。模式中采用轴对称形式的椭球体扰动函数：

$ \theta=\theta_0+\Delta \theta \cdot \cos ^2\left(\frac{\pi}{2} \beta\right) $

(1)

$ \begin{gathered} q_{\mathrm{v}}=q_{\mathrm{v} 0}+\left(q_{\mathrm{vs}}-q_{\mathrm{v}}\right) \cos ^2\left(\frac{\pi}{2} \beta\right) \end{gathered} $

(2)

式中: θ为位温，θ₀为初始位温扰动(设定为2.5 K)，Δθ为中心最大扰动位温，q_v为比湿，q_v0为初始湿度扰动(设定为80%)，q_vs为考虑了位温扰动后的饱和比湿。

$ \beta=\left(\frac{x-x_{\mathrm{c}}}{x_{\mathrm{r}}}\right)^2+\left(\frac{y-y_{\mathrm{c}}}{y_{\mathrm{r}}}\right)^2 $

(3)

式中: x，y为所求目标格点的空间位置; x_c, y_c为扰动场的中心坐标; x_r, y_r分别为扰动场在两个坐标方向的半径。

1.4 初始场

本文采用1999年7月23日美国新墨西哥州发生的一次山地雷暴云电过程作为模拟个例(王宁宁，2012；Shi et al，2015；2016；Tan et al，2017)。该个例的探空数据如图 1所示，大气探测数据源于NCEP再分析资料。本次模拟的雷暴为山地弱雷暴，初始地面温度约为12℃，初始露点温度约为8℃。0℃层位于约550 hPa的高度。模拟总时长为90 min，时间步长为2 s。为了探讨大气冰核不同核化方案对放电过程造成的影响，本文进行了三组敏感性模拟试验，分别为：仅考虑同质核化过程(Ho方案)、仅考虑异质核化过程(He方案)以及同时考虑同质核化和异质核化过程(Co方案)。

2 模拟结果

表 1给出了三种方案中闪电(包括正、负地闪和云闪)频次，可见，Ho方案发生闪电的频次最少，Co方案发生闪电频次最多。

图 2给出了不同方案下闪电放电频次随时间的变化。Ho方案的首次云闪出现在第24 min，唯一一次地闪为正地闪，发生在第25 min。云闪活动在第27 min达到峰值(1 min内发生了21次)，最后一次闪电活动记录出现在第51 min。He方案的首次云闪出现在第19 min，闪电活动最频繁的时间点为第23 min，1 min内记录了46次云闪，最后一次云闪发生在第31 min。地闪活动多于第22 min时发生，1 min内集中发生了6次正地闪，3次负地闪。Co方案中，云闪在第19 min首次出现，于第23 min和27 min分别记录了41次和42次云闪。此外，从第22 min开始观测到负地闪，在第30 min和35 min，负地闪活动尤为频繁，分别达到了48次和58次。通过对比分析可知：He和Co方案的闪电活动发生时间相比Ho方案要早，因此异质核化可能有助于闪电的提前发生。在Ho和He方案中，云闪的发生频次显著高于地闪，而在Co方案中，地闪的次数则显著增加，且均为负地闪。这表明当同质核化与异质核化同时发生时，会促使地闪频率上升。

为了揭示不同核化过程对闪电始发点高度的影响，图 3展示了三种方案不同时间的闪电始发点高度。由图 3可见，Ho方案中闪电始发点高度最高，主要分布在6~8 km高度范围内。其中，第25 min时正地闪的始发点位于6.5 km的高度；而云闪则主要发生在第24~51 min内，其始发点最高和最低高度分别位于8 km和6 km，在7.5 km的高度上记录了24次云闪事件。He方案的闪电始发点主要集中在4~6 km的高度，其中，5 km高度上闪电最频繁，在此高度上共发生了5次正地闪、2次负地闪和高达72次的云闪。10次地闪活动均发生在第22~27 min内，其始发点主要在5~6 km；而在第19~31 min内共发生了148次云闪，其始发点高度则在4~6 km。Co方案共发生了107次负地闪，其始发点最高和最低高度分别位于7 km和4.5 km，时间范围为第22~35 min。此外，该方案中还记录了185次云闪，主要集中在4.5~5.5 km的高度，发生时间在第19~35 min。Co方案有3个地闪的频发高度：在4.5、5.5、6.5 km处分别发生了29、39、14次；而其云闪最频发的高度则是5 km，总计发生了147次。此外，Co方案闪电始发高度的垂直分布最宽广。

综上所述，仅考虑同质核化，雷暴云中产生大量云闪，云闪始发点高度相对高(6~8 km)，放电持续时间长。当仅考虑异质核化时，闪电发生的时间提前，闪电始发点高度分布在4~6 km。当同时考虑两种核化过程，雷暴云中放电提前，同时产生大量地闪，闪电始发点高度在4.5~7 km。

3 讨论分析

不同核化过程决定了初始冰晶的含量与分布，并且影响雷暴云微物理特征。图 4为冰晶最大混合比和数浓度的时间-高度分布，从图中可以看出Ho方案中冰晶从第20 min左右开始产生，混合比大于0.5 g· kg^-1的冰晶主要集中在6~9 km，为云中的低温区(低于-20℃)，与Alizadeh-Choobari and Gharaylou(2017)的研究结果类似，且冰晶的分布较为稳定，整体来看，随着时间的推移Ho方案的冰晶受重力作用逐步下沉。He方案冰晶产生的时间要早于Ho方案，冰晶的混合比明显小于Ho方案，并且中心位置明显低于Ho方案，这与Shi and Tan(2014)的研究相吻合：同质核化主要发生于云顶低温区，而异质核化主要发生在高温区。Co方案的冰晶最大混合比的分布情况整体与He方案相似，但是其冰晶混合比大于4 g·kg^-1的范围更广，高温区冰晶是异质核化产生，而低温区冰晶源于同质核化过程，但上层冰晶存在的时间明显短于Ho方案，可以证明异质核化的存在明显抑制了同质核化。

雷暴闪电放电特征与空间雷暴云电荷结构息息相关。因此，为了解释不同核化过程对闪电放电特征影响的物理原因，本文首先给出了不同方案的起电过程的演变特征。非感应起电目前被认为是雷暴云内最重要的起电过程，其产生原理为霰粒子和冰晶粒子碰撞引起的电荷分离。图 5给出了三种方案下非感应起电率随时间的变化，这代表了冰晶荷电情况。从图中可以看出Ho方案的正非感应起电率最大(67.1987 pC·m^-3·s^-1)，这是因为云顶高度附近(7~9 km)高处同质核化生成大量的小尺度冰晶，这些小冰晶与霰粒子产生碰撞分离形成非感应起电过程。这与Shi et al(2024)的模拟结果是一致的。由于温湿条件在雷暴云不同高度上存在差异，冰粒子在不同温湿条件下所带的电荷极性不同。当温度低于反转温度(-15℃)时，霰粒子荷为负电，而冰晶粒子荷为正电。随着大量冰晶低于反转温度，且云水含量大于0.1 g·m^-3时，粒子荷电情况发生了反转，变成霰粒子带正电，而冰晶粒子带负电。此外，随着雷暴进一步发展，云水含量减少至0.1 g· m^-3以下后，电荷的转移不再满足反转条件，冰晶在-15℃以下依旧带正电(图 4b, 4c)。He方案的正非感应起电率的最大值为14.0989 pC· m^-3·s^-1，负非感应起电率的最大绝对值为108.1587 pC·m^-3·s^-1，Co方案与He方案相似：正非感应起电率的最大值为16.2216 pC·m^-3·s^-1，负非感应起电率的最大绝对值为108.151 pC·m^-3·s^-1。这表明异质核化在高温区产生的冰晶对负非感应起电贡献作用明显, 结果与Shi et al(2024)相同。

本文所使用的三维雷暴云模式涉及到多种水成物粒子的物理变化过程及起、放电过程，其中雷暴空间电荷结构对放电特征作用明显。本文假定闪电下行先导通道发展到距地面1 km处即发生地闪(谭涌波，2006)。图 6给出了雷暴云内空间电荷结构的演变特征，取正、负电荷密度各自的极大值。

Ho方案整体上是偶极性电荷结构，持续时间最久，正电荷区的位置是三个方案中最高的(8 km左右)，并且主正电荷区的电荷量更大，在雷暴发展旺盛的时期负电荷区表现出下沉的趋势，主负电荷区的中心位于5~8 km处(图 6a)。这种现象可能是由于雷暴云发生了贝吉龙过程，即当云层伸展到过冷区，其中存在的大量过冷水滴将通过同质核化产生冰晶粒子并抑制云滴形成。因此Ho方案中的闪电始发点普遍较高。如图 7a所示，当雷暴发展至第24 min时，偶极性电荷结构背景下，雷暴云中首次发生云闪，其高度为7.5 km。此外，Ho方案中雷暴在第25 min时偶极性电荷结构背景下发生了一次正地闪，正先导穿过雷暴中下部的负电荷堆接地(图 7b)。第52 min时，雷暴处于消散阶段，此时雷暴电荷结构仍然为偶极性电荷结构，但其高度有所下降，此时云闪起始高度为7 km高度附近(图 7c)。

He方案中，异质核化导致更早的起电过程，雷暴云在前40 min主要呈现出偶极性和三极性电荷结构(图 6b)。此方案首次放电的时间为第19 min，此时云闪起始于主正电荷堆与主负电荷堆之间，约为5 km高度处(图 7d)。在偶极性电荷结构背景下，正电荷堆和负电荷堆之间产生正地闪，其正先导穿过负电荷堆后接地(图 7e)。此外，在第22 min时电荷结构演变为三极型电荷结构，负先导穿过底部次正电荷堆接地(图 7f)。随后，随着雷暴进一步发展，雷暴空间电荷结构逐渐演变为偶极性，闪电以云闪为主。

Co方案中，异质核化促进起电过程提前出现，同时出现两种核化过程导致雷暴起电量更强(相对于Ho和He方案)(图 6c)，Co方案同样以云闪开始，此时雷暴云电荷结构以偶极性为主，电荷量强，云闪的正负通道更长(图 7g)。在第30 min时雷暴电荷结构演变为三极性，负地闪起始于主负电荷堆与次正电荷堆之间(图 7h)。随后，雷暴电荷结构变为偶极性，闪电以云闪为主(图 7i)。如前文所述，Ho和He方案中的正地闪是在偶极性电荷结构下产生的。由图 6b和图 7e可见，在He方案中，下行的正先导穿过底部的主负电荷堆后接地。而在Co方案中，电荷结构以三极性电荷结构为主。由图 6c可以看到，此时中部的主负电荷堆最大电荷量明显高于上部的主正电荷堆。结合图 3c可以发现，部分负地闪起始于主正电荷堆与主负电荷堆之间，此时正极性先导仅限于在负电荷堆内部移动，而负极性先导不仅横向穿越正电荷堆传播到远处，还能抵达地面形成负地闪。因此，Co方案中负地闪活动的次数显著增加。

4 结论与讨论

本文利用三维雷暴云电耦合模式，选取一次山地弱雷暴个例，研究了三种冰晶核化过程对雷暴云放电特征的影响，探讨了不同核化过程对雷暴云的动力学、微物理过程及电过程的作用。分析结果如下：

(1) 同质核化主要发生于云顶低温区，雷暴云通过发生贝吉龙过程，使其中存在的过冷水滴产生大量冰晶粒子并抑制云滴形成。当温度低于反转温度(-15℃)且云水含量低于0.1 g·m^-3时，冰晶粒子为正电荷。因此，该方案下正电荷区的位置最高，且主正电荷区的电荷量更大，雷暴云易呈现出偶极性电荷结构，并且持续时间较长。故仅考虑同质核化时，雷暴云中产生大量云闪，云闪始发点高度相对高，且放电持续时间长。

(2) 由于异质核化主要发生在雷暴云中高温区，导致更早的起电过程，且产生的冰晶中心位置也较低。随着大量冰晶低于反转温度，且云水含量大于0.1 g·m^-3时，粒子荷电极性发生反转，霰粒子带正电，而冰晶粒子带负电。随着雷暴进一步发展，云水含量减少至0.1 g·m^-3以下后，电荷的转移不再满足反转条件，冰晶在-15℃以下依旧带正电，导致雷暴云易呈现出三极性到偶极性电荷结构的转变。此外，仅考虑异质核化时，闪电发生的时间提前，闪电始发点高度相对较低，同时偶极性电荷结构背景下促进了正地闪的发生。

(3) 当两种核化过程同时发生时，雷暴冰晶空间分布与仅考虑异质核化过程相似，但冰晶含量更高，其导致雷暴起电量更强(相对于仅发生同质核化和仅发生异质核化)。同时考虑两种核化过程时，雷暴云中放电提前，闪电频次最大。同时，雷暴三极型电荷结构背景下主负电荷堆和次正电荷堆之间产生了大量负地闪。该方案中闪电不仅可触发于高温区的主负电荷堆和次正电荷堆之间，而且可以触发于低温区的主正电荷堆和主负电荷堆之间，因此闪电始发点高度范围更大。

近年来，随着对对流云中冰晶的观测和理论研究得到发展，人们对云中冰晶生成过程的认识得到了深化，从而有效揭示了冰晶在云微物理和降水过程的作用。次生冰晶产生于冰晶繁生过程中，例如霰粒子碰撞小水滴后产生冰晶碎片，这些冰晶碎片随后和其他过冷水滴相互作用，使得云内冰晶数浓度不断增加。冰晶繁生过程在对流云中尤为重要，是导致对流云内冰晶数浓度快速增加的最重要原因。那么，对流云中的冰晶繁生是否也会影响雷暴云的起、放电过程？不同的冰晶繁生过程对雷暴云起电过程、空间电荷结构影响如何？目前，这类问题还缺乏深入的探讨，而相应的研究工作也鲜有报道。未来将致力于研究冰晶繁生过程对雷暴云起电过程和空间电荷结构的影响，为雷电预警技术的改善提供重要的理论基础。

此外，本文为了对比同质核化和异质核化过程对雷暴放电特征的作用，综合考虑了接触核化、浸润核化和沉积核化三种异质核化过程，并且与异质核化相关的气溶胶浓度设置也相对高，这可能会导致异质核化作用加强。然而，需指出的是，不同学者通过野外观测以及室内试验的方法得到不同条件下冰核浓度过冷却谱经验公式，其往往表征气溶胶中一种机制或者两种成核机制下冰核浓度，这些经验公式也被许多模式所使用(Fletcher，1962；Cooper，1980；Meyers et al，1992)。这主要由于不同地区的大气冰核特征具有显著差异，其大气环境与核化方式也有所差异，本文综合考虑三种异质核化过程，为一种理想状态的数值模拟研究。着重探讨不同核化过程对放电过程产生的不同影响，虽然模拟结果可能不能完全表征实际雷暴中的闪电活动特征，但不同核化过程对雷暴放电过程影响机制可能具有重要研究价值。同时，我们也认识到，为了更加准确地反映真实情况，未来研究应进一步优化模拟参数化方案，并尽可能采用观测事实对模拟结果进行验证。