DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 摘 要 金融系统中的多种关系虽然提高了金融服务的便捷性，但也使得风险传导的渠道更为复杂。本文致力于构建包含银行同业拆借关系、银企间借贷关系和企业间担保关系的银企间多层金融网络。基于银行和企业两类主体的资产负债表内及表外业务，阐述各主体间的债务债权关系，明确风险构成要素。综合考虑初始冲击导致的个体风险、信用风险以及市场流动性风险等因素，从分层传导和分步传染两个角度描述风险传导机制，提出网络结构中的传染效应。同时，基于Debtrank算法提出风险传导过程中的测度方法，构建对个体与整体风险状态的测度思路。 关键词 :多层金融网络，银企间网络，风险传导机制，风险测度建模 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 2. 基于银企间多层金融网络的模型构建 2.1. 银企间网络模型的构成要素 本文建立一个以银行和企业两类主体为节点的多层网络，包含 N + M 个节点，分别是N家银行和M家企业。如图1所示，第1层是银行间的同业拆借关系层(实线指代拆借关系)，理论上存在与其他节点无任何关系的节点，故假设至多存在N个持有风险资产的银行b，第2层是银行与企业的借贷关系层(点线指代借贷关系)，同理，至多存在N个持有风险资产的银行b和M个企业f，第3层是企业向银行借款而存在的担保关系层(虚线指代担保关系)，同理，至多存在M个持有风险资产的企业f。定义构建的复杂多层金融网络模型为 G ( E ， V ) ，集合V中的各节点表示银行b、企业f；当借贷关系或者担保关系存在时，银行之间、银行与企业之间或者企业之间存在联系，记作l，属于集合E。 金融风险是在金融市场运营中的不同参与者之间相互作用中产生的，本文建立银企间多层网络模型来描述不同参与者之间的关联性，关键在于确定不同参与者之间的债务债权关系。我们用邻接矩阵 A j i ， γ = { a j i ， γ } N * N ( M ) 来表示各节点之间的关系， γ 表示不同层：银行间同业拆借层(第1层)、银企间借贷层 (第2层)、企业担保层(第3层)； a j i ， γ 表示第 γ 层上的第i个节点拥有第j个节点的债权金额，相当于i的风险敞口，按惯例方向是从节点j指向节点i；统一定义 a i i ， γ = 0 表示不存在自己对自己的债权。 对于银行间的同业拆借层(第1层)，银行j向银行b借款，因为存在无法还款的可能性，则该笔借款 a j b ， 1 对银行b来说是一项风险资产；再者银行间的同业拆借往往期限较短，目的在于平复银行间流动性盈亏，解决临时性的资金短缺，因而银行间同业拆借结构与银行间同业拆借行为同步演化，是一个动态网络层。同理，企业f向银行b借款，因为存在无法还款的可能性，则该笔借款 a f i ， 2 对银行b来说也是一项风险资产；且无论长期借款还是短期借款，均是一个不断演化的过程，因而银企间借贷层(第2层)也是一个动态网络层。另外，如果企业f向银行b借款过程中存在第三方担保企业1k，即存在担保关系，网络中出现企业间担保层(第3层)。担保方所承担的债务属于从属借贷合同的债务，也是一个动态演变的过程，所以企业间担保结构也是一个动态网络层。 图1. 简化的网络结构图 2.2. 银行和企业两类主体的微观结构 在微观层面上，该多层网络中的所有节点都满足资产负债表的特征，且不同节点的资产负债表之间存在一定的联系。因此本文将使用简化的资产负债表代表每个节点的状态，以期能够量化整个网络的动态变化。 银行的资产 V b ( 0 ) 细分为银行间资产 V J B = ∑ j = 1 N a j b 、银企间资产 V F B = ∑ f = 1 M a f b 、其他资产 V O B ( 0 ) = φ 1 V b ( 0 ) ；银行的负债和权益栏分为银行间负债 L J B 、其他负债 L O B 和权益 E b ，如表1所示。对于银行b的资产负债表中每个同业资产 a j b ，在银行j的资产负债表中都存在与之对应的负债，即 L b j = a j b ，故银行间负债 L J B = ∑ j = 1 N a b j ；并设定其他负债 L O B = V b ( 0 ) − L J B − E b ( 0 ) ；权益 E b 作为损失缓冲器， E b ( 0 ) = ϕ 2 V b ( 0 ) ；且 V J B + V F B + V O B = L J B + L O B + E b 关系恒定。 同理，简化的企业资产负债表如表2所示，将企业的资产 V f ( 0 ) 细分为流动性资产 V L F = φ 3 V f ( 0 ) 和非流动性资产(其他资产) V O F = ( 1 − φ 3 ) V f ( 0 ) ；企业的负债和权益栏分为银企间负债 L B F = ∑ b = 1 N a f b 、其他负 债 L O F 和权益 E f ；显然两张资产负债表中，企业的银企间负债与银行的银企间资产是相互对应的；权益 E f 作为资产缓冲器， E f ( 0 ) = ϕ 4 V f ( 0 ) ；且 V L F + V O F = L B F + L O F + E f 关系恒定；另外，由于担保业务属于表外业务，而企业对外提供担保实际上形成了或有负债，因而一旦企业被要求履行担保义务时，体现在资产负债表上为资产受损，这里假定流动性资产受损。 Table 1 Simplified bank balance shee 资产 负债及权益 银行间资产 V J B 银行间负债 L J B 银企间资产 V F B 其他负债 L O B 其他资产 V O B 权益 E b 表1. 简化的银行资产负债表 Table 2 Simplified firm balance shee 资产 负债及权益 流动性资产 V L F 银企间负债 L B F 非流动性资产(其他资产) V O F 其他负债 L O F 权益 E f 表2. 简化的企业资产负债表
system: [SYSTEM]

user: 5. 结论 本文首先分析了金融网络系统中银行间同业拆借、银企间借贷和企业间担保三种关系，基于银行和企业这两类主体间的复杂关系，构建了多层网络，并基于银行和企业两类主体的资产负债表内及表外业务，阐述各主体间的债务债权关系；其次，本文从分层和分步两个角度分析了风险传导过程，对银行间市场、银企间市场和企业担保圈的风险链条式传染以及系统性风险的成因进行了解释；最后，基于Debtrank算法构建了风险测度的思路，分别提出对个体和整体的测度，以期能够使金融系统中的银行与实体经济体系中的企业更加深刻地体会宏观环境从而进行有利的经济行为，更为监管部门制定风险防范措施提供理论依据，从而提高金融风险政策调控效果。
system: [SYSTEM]

user: 摘 要 采用旋转压实仪成型了三种沥青混合料圆柱形试件，并钻芯取得动态模量试验试件。使用UTM-30试验系统在不同试验温度和荷载频率条件下，对三种不同类型沥青混合料进行了动态模量试验。分析了沥青混合料动态模量与试验温度、荷载频率、沥青种类及混合料类型的变化关系。利用西格摩德(Sigmoidal)函数构建了沥青混合料的动态模量主曲线。研究表明：在相同试验条件下，试验温度升高和荷载频率减小使动态模量减小，温度降低和荷载频率增大使动态模量提高。采用高模量沥青对SMA-13沥青混合料动态模量的提高能力优于SBS改性沥青。20℃试验温度下，沥青混凝土的动态模量主要与高模量改性沥青胶结料有关。动态模量主曲线与动态模量实测值有较好的相关性。 关键词 :道路工程，沥青混合料，动态模量，主曲线 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 路面承受的荷载主要是行驶中的车辆对路面施加的动态荷载，并且我国大多数铺面材料为沥青混合料，其为黏弹性组合体。材料在一定加载形式及加载频率下的动态响应称为动态模量，相比于抗压回弹模量，路面结构的力学响应能较好地通过动态模量这一指标进行表征。因此沥青路面设计新规范 [ 1 ] 将沥青混合料的动态模量作为关键指标之一。因此，对动态模量这一设计指标的研究对沥青路面结构设计及路面性能耐久性有较高的价值。 徐志荣 [ 2 ] 等利用静压成型切割试件，通过无侧限抗压测试方法得到试件动态模量数据，并对比了不同成型方式试件的空隙率对动态模量产生的差别进行了分析。结果表明：静压成型更贴合动态模量测定，试验结果离散性较小。文献 [ 3 ] 在规定温度下用超声波直接测定沥青混合料动态模量的试验过程，并将超声透射法测得的动态模量与标准动态模量试验值进行了比较。研究表明在频率为2，5，8和10 Hz的情况下，超声透射计算得到的模量大于标准动态测试得到的模量。王昊鹏 [ 4 ] 等通过简单性能试验机得到了两种高模量改性沥青混合料的动态模量试验值，求得了动态模量主曲线，并通过主曲线预测了试验中两种高模量沥青混合料的高温稳定性及低温抗裂性。文献 [ 5 ] 研究了沥青混合料动态模量的拟合模型，利用Arrhenius公式，基于时间–温度等效原理，得到了沥青混合料的主曲线。李强 [ 6 ] 等采用不同受力方式对3种沥青混合料进行试验，研究沥青混合料动态模量主曲线与受力模式、围压、应变水平三个因素的关系，结果表明：沥青混合料受力状态、围压及应变水平影响其动态模量。动态模量试验值与各因素之间存在较好相关性。大量研究也在动态模量的预估模型取得丰硕的成果 [ 7 ] [ 8 ] [ 9 ] [ 10 ] [ 11 ] 。 结合以上研究，本文主要分析不同沥青混合料在不同试验温度和荷载频率情况下，混合料动态模量变化情况。对沥青混合料的动态模量主曲线进行研究。
system: [SYSTEM]

user: 2. 原材料与试验方法 2.1. 试验原料 试验采用高模量改性沥青(HMB)及成品SBS改性沥青，沥青相应实验检测指标均按照相应规范严格执行其各项技术指标见表1所示。 Table 1 Asphalt related technical indicator 沥青种类 针入度/0.1 mm 延度/cm 动力粘度/(Pa.s) 软化点/℃ HMB 43.4 57.4 4.425 93.1 SBS I-C 75 32.8 2.878 70.5 表1. 沥青相关技术指标 试验采用的粗为优质的玄武岩集料，相应技术指标如表2所示。细集料为优质机制砂，技术指标见表3。矿质填料选用亲水系数小，安定性好的石灰石矿粉，各项性质均符合要求。 Table 2 Technical specifications of coarse aggregate 技术指标 集料规格 10~20 mm 10~15 mm 5~10 mm 3~5 mm 技术要求 毛体积相对密度 2.688 2.662 2.630 2.605 / 表观相对密度 2.727 2.714 2.695 2.737 ≮2.45 吸水率/% 0.530 0.718 0.955 1.786 ≯3.0 针片状含量/% 8.5 10.5 8.0 / / 与沥青的粘附等级/级 5 5 5 5 ≮4 表2. 粗集料技术指标 Table 3 Technical Specifications of Fine aggregat 技术指标 表观相对密度 坚固性 > (0.3 mm部分) 亚甲蓝值 棱角性(流动时间) 2.740 0.4 0.8 33 表3. 细集料技术指标 2.2. 沥青混合料组成设计 试验涉及到的三种沥青混合料分别为(SMA-13级配 + HMB)，(SMA-13级配 + SBS改性沥青)，(AC-20 + SBS改性沥青)。SMA-13级配矿料组成为10~15 mm：5~10 mm：3~5 mm：机制砂 = 37:37:6:10:10，纤维含量按石料3‰添加。AC-20级配矿料组成为10~20 mm：3~10 mm：机制砂：矿粉 = 38:28:28:6，最后确定三种沥青混合料油石比分别为6.1%，7.9%及4%。按最佳油石比成型马歇尔试件，并测定试件相应体积指标。三种沥青混合料技术指标如表4所示。 2.3. 试验方案 试件成型使用美国PINE公司AFG1C型旋转压实仪，旋转压实75次成型高度为170 mm，直径为150 mm地圆柱形试件。待试件固化形成强度后钻芯取样，切割形成高度为150 mm，直径为100 mm的混合料试件。 Table 4 Technical specifications of asphalt mixtur 混合料种类 表观相对密度 空隙比/% 矿料间隙率/% 沥青饱和度/% 稳定度/KN 流值/mm SMA-13 + SBS改性沥青 2.417 3.8 17.8 78.7 10.5 34.0 SMA-13 + HMB 2.338 3.6 17.9 74.1 13.4 30.5 AC-20 + HMB 2.565 4.9 17.23 65.20 12.2 27.8 表4. 沥青混合料技术指标 试件制作完成后，在UTM-30多功能沥青混合料测试机上进行动态模量试验。试验中对试件施加轴向正弦压应力荷载，平均微应变范围控制在85~115 μm。三种沥青混合料在5℃、20℃及40℃温度下分别采用0.1 Hz、0.5 Hz、1 Hz、5 Hz、10 Hz和25 Hz加载频率测定试件的动态模量及相位角，每个试验温度下进行三组平行实验。考虑到沥青混合料的黏弹性及混合料结构内温度和荷载频率对沥青混合料力学响应的影响，本文分析不同试验温度、不同荷载频率及不同混合料级配对沥青混合料动态模量和相位角的影响。
system: [SYSTEM]

user: 3. 试验结果分析 3.1. 不同沥青对试验结果的影响 在20℃试验温度下，(SMA-13级配 + HMB)及(SMA-13级配 + SBS)沥青混合料在不同荷载频率下的动态模量和相位角变化如图1(a)、(b)所示。 图1. 20℃下不同沥青混合料动态模量试验结果对比 由图可知(a)中沥青混合料采用高模量沥青动态模量明显高于SBS改性沥青，(b)中相位角有所减小。说明了在相同温度及相同级配下，试验所用的高模量沥青混合料相比于SBS改性沥青混合料模量有相应提高，抗变形能力有一定程度的提升。高模量改性沥青的掺入相比于SBS改性沥青提高了沥青混合料的动态模量，提高了荷载作用下的不可恢复变形，延长路面服役周期。原因在于高模量改性沥青相比于SBS改性沥青，增强了集料与沥青的粘结能力，从而使混合料抗变形能力提高，动态模量随之增大。高模量沥青混合料有较好的应用前景。 3.2. 不同加载频率对试验结果的影响 在相同试验温度20℃不同加载频率下，三种沥青混合料动态模量试验结果见表5。 Table 5 Dynamic modulus test results at different loading frequencie 加载频率 SMA-13级配 + HMB SMA-13级配 + SBS AC-20级配 + HMB 动态模量相位角 动态模量相位角 动态模量 相位角 0.1 Hz 2836 28.60 1278 27.71 4420 24.50 0.5 Hz 4673 27.12 2051 29.37 6304 20.44 1 Hz 5724 25.94 2562 29.52 7085 19.80 5 Hz 8628 22.09 4349 27.77 9294 18.78 10 Hz 10109 19.98 5410 26.30 10610 15.19 25 Hz 12124 17.51 6946 24.22 12106 14.99 表5. 不同加载频率下动态模量试验结果 由表5中数据可知，在相同试验温度下，试件加载频率越大其动态模量值越大，荷载频率的增加，相位角逐渐减少。结合图1(a)可知，在较低正弦荷载频率下(0.1~1 Hz左右范围内)，动态模量增长较快。相位角曲线在较低正弦荷载频率下切线斜率较大，变化较快。这是由于沥青混合料这种黏弹性材料在较大的正弦荷载频率作用下，轴向应变滞后于应力，试件被加载后，在力的作用下试件不会被立即压缩，卸载后也不会立即回弹，所以较大正弦荷载频率作用时发生的应变较小表现为较大的模量，同时增长速度也会降低。而在较低的正弦荷载频率作用下，应变滞后应力的效果减弱，所以表现为较小的模量及较快的变化速度。 3.3. 不同混合料组成级配对试验结果的影响 为了分析不同混合料组成级配动态模量和相位角的关系，本文选取了两种普遍使用的沥青混合料：沥青玛蹄子碎石混合料(SMA-13 + HMB)及AC密级配沥青混合料(AC-20 + HMB)在20℃下进行了测试。试验结果如图2(a)、(b)所示。 图2. 20℃下不同级配混合料动态模量试验结果对比 由图2(a)可知，AC密级配沥青混合料在各荷载频率下与沥青玛蹄子碎石混合料相比，动态模量相差不大。(b)中AC密级配沥青混合料相位角略低于沥青玛蹄子碎石混合料。这是由于试验温度较低与常温接近，在此试验温度下，沥青胶结料对沥青混凝土的劲度起关键的作用，混合料的矿料骨架的嵌挤作用影响不大。两种混合料所用沥青均为高模量改性沥青(HMB)，沥青胶结料性能相近，AC密级配沥青混合料中沥青胶浆所占体积百分比较大。从而结果表现为AC密级配沥青混合料动态模量略高于沥青玛蹄子碎石混合料。同时可以说明AC型级配的高模量沥青混合料具有现实使用意义，路面结构抗变形能力较普通AC型级配的沥青混合料有较大的提升。 3.4. 不同温度对试验结果的影响 在5℃和20℃试验温度下，(SMA-13 + HMB)沥青混合料分别在的动态模量和相位角变化如图3(a)、(b)所示。 图3. 不同温度下沥青混合料动态模量对比 由图3可知，沥青混合料的动态模量随着温度的升高而降低，随着温度的升高相位角有所升高。在10 Hz的荷载频率下，(a)中20℃混合料动态模量时5℃的2倍左右。原因是沥青混合料是一种具有粘弹性的复合材料，随着温度升高，沥青结合料逐渐变软，沥青黏度下降，黏性是此时沥青混合料主要的力学响应，导致沥青混合料对荷载反应滞后过程减弱，沥青混合料由弹性逐渐向塑性转化，抗变形能力降低，混合料动态模量降低。
system: [SYSTEM]

user: 4. 动态模量主曲线分析 针对组成设计相同的沥青混合料不同温度条件的动态模量，每个温度条件下进行试验去测定动态模量，试验量较大且过于繁琐。为了方便探究不同温度下和加载频率下的沥青混合料的动态模量，可利用现有动态模量的试验数据，采用时温置换原理，通过平移形成一条在参考温度下的曲线，该曲线称为主曲线 [ 12 ] [ 13 ] 。通过主曲线，可以得到在某一温度及加载频率下的动态模量值。 本文以5℃、20℃及40℃下沥青玛蹄子碎石混合料(SMA-13 + HMB)在不同加载频率下的动态模量数据，选用20℃作为参考温度，根据时间-温度置换原理，移位因子lgα(T)通过不同温度动态模量试验值下非线性最小二乘法拟合得到，最后求得西格摩德(Sigmoidal)函数。最终确定高模量沥青玛蹄子碎石混合料动态模量主曲线。拟合函数为西格摩德函数，其表达式为 lg | E ∗ | = δ + α 1 + e β + y ( lg ω r ) 利用该函数拟合试验所得的高模量沥青玛蹄子碎石混合料动态模量数据，依据上述方法求得 lg α ( T ) ， lg ω r 缩减频率根据公式 lg ω r = lg ω + lg α ( T ) 得到。最后得到参考温度下高模量沥青玛蹄子碎石混合料的动态模量主曲线方程为 log | E ∗ | = 2.4890 + 3.6025 1 + e 0.5870 − 0.6505 ( log ω r ) 参考温度下高模量沥青玛蹄子碎石混合料的动态模量主曲线如图4所示。 图4. 沥青混合料动态模量主曲线 由图4中动态模量主曲线与动态模量实测值对比可知，两者相关性较好。表明利用主曲线来预测不同温度和加载频率下沥青混合料动态模量值较为方便，试验数据更加合理精确。可精确快速得到在大范围内的试验温度和加载频率下的混合料动态模量。
system: [SYSTEM]

user: 5. 结论 1) 在相同的荷载频率下，动态模量随着温度的升高而降低，相位角随着温度升高而增大。在相同的试验温度下，相位角随着荷载频率的增大而升高，相位角随着荷载频率增大而降低。 2) 高模量改性沥青对沥青混合料动态模量的提升优于SBS改性沥青。在温度及荷载频率相同情况下，(SMA-13 + HMB)沥青混合料动态模量明显高于(SMA-13 + SBS改性沥青)沥青混合料。 3) 常温条件下，高模量改性沥青胶结料对沥青混凝土的动态模量起关键的作用。在20℃荷载频率相同情况下，(AC-20 + HMB)沥青混合料动态模量(SMA-13 + HMB)沥青混合料相差很小，源于其沥青结合料相同。 4) 高模量沥青混合料的动态模量主曲线与动态模量实测值有较好的相关性，可由主曲线预测不同试验温度和加载频率下的混合料动态模量。
system: [SYSTEM]

user: 摘 要 二战后军事手段在解决外交关系中逐渐减少使用，取而代之的是经济金融制裁，通过国际贸易、资金流动等影响相关国家产业链与相关领域。经济作为一国生存和发展的重要命脉，成为美国控制国际市场、维护其世界霸权主义的主要手段，而经济金融制裁则成为美国重要的对外政策。近年来，美国以维护国际安全为由对朝鲜发起多轮经济金融制裁，制裁的广度和深度大大超过以往，特别值得关注的是以联合国安理会的名义实施。因此，选择从近几年美国对朝鲜的经济金融制裁入手，分析美国对朝鲜经济金融制裁演进历程，由此对美国、朝鲜、第三方国家政治、经济等方面影响，最后结合我国与朝鲜差异，针对我国对外开放以及经济实力从对外开放格局、国际货币体系、人民币国际化等方面提出建议。 关键词 :美国，经济金融制裁，朝鲜 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 经济全球化给发展中国家和地区提供众多机遇，各国产业向规模化发展、资金流转增加，各经济体间依赖性加大，国际经济金融制裁次数增多。本文从近两年美国对朝鲜的经济金融制裁的鲜活案例着手，通过美国对朝鲜数轮经济金融制裁的制度演进、影响，寻求美国对他国制裁的一般手法和措施，从中得到一定的启示。目的在于，对未来美国可能对中国采取的金融制裁，我国能借鉴国际经验，快速反应，减少制裁带来的负面效应、降低我国国际贸易风险、加强本国经济稳定性和安全性。
system: [SYSTEM]

user: 2. 文献综述 经济制裁作为一种国际干预手段，5世纪后期就有雅典城邦对麦加拉城邦进行经济封锁的历史之鉴，二战结束后逐渐成为国家在进行国际制裁的有效举措。张曙光在《经济制裁研究》一书中提出经济制裁外交，即国家间的经济制裁关系构成现代国家外交的一种政治模式 [ 1 ] 。在他的另一篇文章中指出，对外经济制裁成为了一个介于立即采取军事行动和毫无作为之间的替代选择 [ 2 ] 。其影响范围不仅局限于国内金融体系，陈捷等(2019) [ 3 ] 认为国际金融制裁对离岸金融同样会产生显著的影响。 徐以升和马鑫(2014) [ 4 ] 认为经济全球化的今天，金融制裁被认为是最严厉的经济制裁形态，金融制裁通过限制被制裁方资金融通的方法，强行逼迫受制裁国停止与发起国外交政策或经济利益相悖的政策行为。随着美国对外经济金融制裁的增多，这方面的研究也逐渐增加，学者从不同的角度对这一现象进行了剖析。 案例分析视角。Mahmoud Jameel Jdeed (2013) [ 5 ] 对伊拉克和伊朗案例比较、采用博弈论模型，理论与实际结合重点研究国际经济制裁对被制裁国经济发展以及决策选择。Evsey Gurvich和Ilya Prilepskiy (2015) [ 6 ] 从美国对俄罗斯制裁的案例，评估制裁如何影响了对财务账户基本组成部分，运用模型详细分析制裁对俄罗斯资本流动等基本宏观经济指标以及国际收支构成部分的影响。ORLOVA， NATALIYA V (2016) [ 7 ] 在对俄罗斯案例分析中主要从银行层面入手，分析了美国对俄罗斯金融制裁的影响，提出银行面对金融制裁需要控制资本质量，否则当大量政府资源都用于补修银行漏洞，会导致长期经济增长滞缓。Lee， Jong-Woon和Gray， Kevin (2017) [ 8 ] 通过对中朝边境经济案例的考察，分析金融制裁对朝鲜的双边经济关系，提出制裁措施只是增加了规避经济活动的范围，进一步全面地认识金融制裁作为外交政策工具的不足之处。 法律法规视角。黄风(2012) [ 9 ] 首次从美国法律角度，归纳出美国实行国际金融制裁制度依据的理由和实现形式，向金融监管机关提出对金融制裁的研究与合规指导，保护我国巨额资产的安全。陈宇曈和成弋威(2017) [ 10 ] 对美国金融制裁类型、法律依据方面，联系案例提出国内银行操作过程中依据法律风险的防范措施。葛淼(2018) [ 11 ] 则提出美国实行单边金融制裁在并未违反国际强行法，中国更多的是需要提升本国实力、构建完善的反制裁机制。 运行机制视角。叶文辉(2015) [ 12 ] 通过美国金融制裁的运作体系总结出其制裁的主要手段包括，冻结被制裁国在美资产、限制或禁止与被制裁国间的金融往来、限制被制裁国在美市场上投融资、罚款与没收被制裁国资产、要求第三方制裁。刘建伟(2015) [ 13 ] 在对美国特有的金融制裁逻辑基础上进一步分析了其完善的金融制裁制度体系，包含了法律体系和组织运作制度，成为美国倚重的制裁方式。 其他视角。Carter， Barry E.和Farha， Ryan M. (2013) [ 14 ] 从经济制裁的久远历史，结合国际贸易扩张导致的金融制裁和经济制裁活动，论述了美国金融制裁在整个经济发展过程的演变意义。李伟和曾令勋(2013) [ 15 ] 在探讨国际话语权时提出，在国际市场掌握的权力大小决定了一国的国家利益和政策在国际层面的影响力、引导力和支配力。马鑫和许钊颖(2015) [ 16 ] 从手段、目标、制裁效果三个方面研究美国对俄罗斯的金融制裁，提出美欧对俄制裁的“温水煮青蛙”模式体现出现今大国博弈的新模式，对俄经济的挫伤和美国面临的国际“去美元化”的软约束会是一个漫长的战略博弈。
system: [SYSTEM]

user: 3. 美国对朝鲜制裁政策动因及措施 美国一直以人道主义和安全为由对其他国家采取单方面的制裁包括军事和金融制裁。对于朝鲜坚持核试验，美国依据《与敌国贸易法》、《全国紧急状态法》、《国际紧急经济权利法》、《爱国法》等金融制裁的法律体系，当美国的国家安全、外交政策和经济利益受到严重外来威胁时，美国有权针对有关外国、外国组织或者个人采取经济金融制裁措施(叶文辉，2015) [ 17 ] 。从2006年朝鲜开始核试验，美国一直以朝鲜核试验严重威胁美国及世界安全为由，单方面或推动联合国安理会对朝鲜采取数轮经济金融制裁，且强调会对其不断采取制裁措施直至朝鲜放弃核试验。 但对于美国多次颁布的制裁措施并没有能阻止朝鲜放弃核试验，2016年1月朝鲜成功的进行了第一次氢弹试验，再次引发美国对朝鲜的严重不满，同年2月奥巴马签署一份制裁朝鲜法案，再次拉开美国对朝鲜长期经济金融制裁的序幕。2016年3月落实安理会对朝鲜2270号决议制裁朝鲜法案，与美国对致力于对朝鲜政府持续施压的长期政策一致。2017年截止12月份美国对朝鲜共采取了13次制裁，旨在切断资金来源。 通过对美国制裁朝鲜的具体制裁手段整理分析可以看出，美国对朝鲜制裁途径包括三种：一是美国单方面对朝鲜颁布制裁法案或是政策；二是美国推动国际组织通常是联合国安理会，以维护国际安全为由通过制裁朝鲜决议；三是制裁与目标国家有重要贸易、金融合作的国家，达到间接制裁目标国的目的。法案或是政策内容涵盖限制贸易进出口、禁止目标国公民赴海外打工、冻结在美境内资产、禁止与被制裁国经济往来。从经济上封锁被制裁国，限制朝鲜重要的进出口货物，减少朝鲜国民经济收入。从金融上限制流动资金，禁止美国公司及个人向朝鲜投资、融资。
system: [SYSTEM]

user: 4. 美国对朝鲜制裁的相关影响 (一) 对朝鲜的影响 基于美朝间经济关联较少的实际情况，有学者提出美国对朝鲜经济金融制裁并不会动摇朝鲜强军治国政策，无论是利用对美国境内领土、公民权力的行使采取贸易限制、资产冻结、经济合作中断，还是通过联合国安理会或第三方国家经济援助的方式，就目前状况显示，从2006年起跨12年的制裁对朝鲜经济都未产生重创，这也是朝鲜至今都没有放弃核试验的重要原因之一。 首先，依据世界贸易组织数据(图1)，从2006年朝鲜首次开始核试验后美国采取针对性的一系列制裁措施，从整体上对朝鲜的影响不明显，相较于2000年，朝鲜进出口额整体呈现增长状态。从局部年份来看，会有进出口下降，而2008~2009年的下降不排除受到全球金融危机影响可能性，2009年之后便开始上升。但从2014年开始，朝鲜的进出口都呈现下降的趋势。韩国产业银行数据显示，朝鲜2017年1~6月出口总额为8.44亿美元，同比减少24%。随后，联合国安理会在2017年8月通过的2371号决议被称为“迄今为止最严厉”的对朝制裁，全面禁止朝鲜出口煤炭、铁、铁矿石、铅、铅矿石以及水产品，这些货物约占朝鲜每年出口额的1/3。结合2016年朝鲜“七大”中金正恩明确提出“经济强国”目标，表明朝鲜对国际社会持开放和接受的态度，结合已有数据，未来美国经济金融制裁对朝鲜贸易额的影响可能会加大。 图1. 朝鲜进出口变化趋势 其次，美国禁止其境内企业、个人与朝鲜任何经济往来，且禁止扩大已有的经济合作规模，但是朝鲜与美国经济合作很少，朝鲜对美国的依赖度较低。因此，美国中断与朝鲜合作对朝鲜影响有限。美国学者则认为，虽然美国对朝鲜制裁无法造成朝鲜经济明显波动，但从信号理论分析，对朝鲜不断采取措施可以向其他国家及朝鲜传达美国对朝鲜的不满。 最后，美国推动联合国安理会或通过制裁第三国家，切断对朝鲜的经济和技术支援。大韩贸易投资振兴公社称，2016年朝鲜对外贸易依赖排名前四的国家为中国、俄罗斯、印度和泰国，依赖度分别92.5%、1.2%、0.9%、0.8%。这也就不难理解美国通过媒体、采取制裁向中国施压的意图。美国总统特朗普近期在推特宣称中国不顾世界安全仍与朝鲜贸易，以2017年第一季度累积数据增长40%为由，片面的质疑中国执行联合国安理会法案的态度。对此，我国海关总署表示，自2017年3月份以来，我国自朝鲜月度进口值连续四个月大幅下降，以美元计价，3月同比下降36.5%，4月同比下降41.6%，5月同比下降31.6%，6月同比下降28.9%，一定程度上对朝鲜国民经济产生负面作用。 (二) 对第三方国家的影响 经济金融制裁能够产生效用的原因是国家间存在依赖性，依赖程度的高低影响一国解决外部制裁成本高低；同时不同国家脆弱性差异反映一国寻求替代贸易或投融资合作国能力不同(王玄，2011) [ 18 ] 。因此，二战后美国经济金融制裁出现辅助性制裁和针对性制裁的特点，针对性制裁主要制裁对象是被制裁国政府领导人和颁布不被国民接受政策的社会团体，而辅助性制裁是在针对性制裁基础上迫使更多国家参与制裁，降低被制裁国寻求国际可替代合作国机率，提高制裁效率(计宏亮，2008) [ 19 ] 。由于美国在执行单方面制裁时无法任意得到其他国家支持，美国国会投票或总统签署的相关法案中会执行“二级制裁”，对目标国以外的第三方国家进行制裁，迫使其他国家加入美国制裁行列，形成多边制裁的现实状态(颜剑英、熊伟，2005) [ 20 ] 。此外，美国积极推动联合国安理会通过制裁朝鲜决议，禁止各国与朝鲜贸易合作和提供技术、经济支持，使得制裁合理化。 中国作为朝鲜的临近国家，与朝鲜的传统友谊导致我国在对朝制裁过程中受到影响。对我国最直接的影响一方面是美国为切断朝鲜经济来源和外贸收入，将与朝鲜有贸易往来的中国公司列入“黑名单”；另一方面是冲击中朝经贸关系，进出口朝鲜货物的通关检验更加严格，增加了时间成本和贸易成本。除此以外，对于当前复杂严峻的半岛局势，中国严格遵守联合国的对朝决议，但反对单方面对朝鲜制裁，美国则以维护世界安全和人权为由，宣扬中国在面对朝鲜核试验问题立场不明确，不能支持其实施制裁朝鲜措施，用朝鲜核问题强行向中国政府施加舆论压力。特朗普政府认为中国可以对朝鲜能发挥更大影响力，继而将朝鲜核试验问题的责任推给中国。
system: [SYSTEM]

user: 摘 要 反应性是一个间接测量的参数，文章回顾总结了费曼方法进行反应性测量的基本原理及其物理意义。就该方法在反应堆临界时的发散，两次计数的时间间隔造成计数丢失效应，使用多个探测器测量进行数据综合，中子源释放相关的中子等新问题进行了分析讨论。对费曼方法应用时仍存在的问题提出了看法。 关键词 :反应性，费曼方法，加速器驱动次临界系统 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 测量反应堆反应性的方差均值比方法最先是由费曼提出的，所以又叫费曼方法 [ 1 ] 。即使名义上的零功率反应堆，也仍然有大量的中子存在。堆中子统计规律的一个重要特点就是堆内中子不是完全互相无关的。 有些中子是由同一裂变链产生的，它们有着共同的“祖先”，在时间上有着相关性。如图1所示，中子a与中子b、c相关，与中子d不相关。这种相关性的存在，使得中子探测器计数率的统计特性偏离了泊松分布。 图1. 中子裂变链示意图 将缓发中子看作是相互独立无关的，和一般的外中子源一样是服从泊松分布的本底，忽略缓发中子之间的相关性，则中子计数的方差与平均值之比为 [ 1 ] σ 2 m = 1 + ε D v ρ p 2 ( 1 − 1 − e − α T α T ) = 1 + Y (1) 其中， D v = v ( v − 1 ) ¯ v ¯ 2 ，v为每次裂变产生的中子数， ε 是探测器的效率，T为测量统计时间， α 为瞬发中子衰减常数， ρ p 为瞬发反应性。 对于完全随机不相关的中子，它们的统计特性服从泊松分布，其方差与平均值之比为1，即 Y = 0 。若中子探测器在某个时刻的计数可能与以前的计数是相关的，即某个时刻引起计数的中子与以前引起计数的中子是属于同一个裂变链，则 Y > 0 。所以Y代表的是堆内中子相关部分，即为堆内中子统计特性偏离泊松分布程度的测度。 对于不同的统计时间T测量出对应的Y，然后利用最小二乘法可以得到瞬发中子衰减常数 α 。根据 α = β − ρ l ，得到反应性 ρ 。其中l是瞬发中子平均寿命， β 是缓发中子份额。 对于以U 235 为燃料的热中子水堆，统计时间T一般在<10 ms的范围内选择，如10 ms、5 ms、1 ms、0.5 ms、0.1 ms等 [ 2 ] 。 当考虑缓发中子的衰减过程，分析缓发中子存在的互相关中子，则相关中子测度Y为 [ 3 ] Y = ε D v ∑ i = 1 7 2 A i α i H 0 ( α i ) ( 1 − 1 − e − α i T α i T ) (2) A i 、 α i 由下面的零功率传递函数决定， H 0 ( s ) = ∑ i = 1 7 A i α i + s 决定。文献 [ 4 ] 试验检测了缓发中子对中子相关测度Y的影响。 测量反应性的Feynman方法提出的时间应该说比较久远了，但是随着应用的不断深入，人们发现了一些新的问题，并着力进行了研究。
system: [SYSTEM]

user: 2. 临界发散的问题 当反应堆临界时，零功率传递函数的最后一个极点 α 7 为零，方差均值比表达式(2)会出现发散。Bennet [ 3 ] [ 5 ] [ 6 ] 提出如下方法来测量堆内中子相关程度，以避免临界时发散的问题。 〈 ( c k + 1 − c k ) 2 〉 2 〈 c k 〉 = 1 + W (3) 其中 W = ∑ i = 1 7 H 0 ( α i ) ( 1 − 3 2 + 1 2 e − 2 α i T − 2 e − α i T α i T ) ， 〈 ⋅ 〉 表示均值， c k 为第k个时间间隔T内的计数，W反映了堆内中子相关程度。 Wallerbos [ 7 ] 发现由于试验测量时间是有限的，实际试验中相关中子测度Y并不会发散，提出了修正后的有限测量时间模型 〈 Y ( N T ) 〉 = 2 ε D v ∑ i = 1 7 A i α i H 0 ( α i ) × [ f ( α i T ) − 1 N f ( α i N T ) ] (4) NT为测量时间， f ( T ) = 1 − 1 − e − α T α T 。如果测量时间无限， N → ∞ ，则 Y ( T ) 的期望值 lim N → ∞ 〈 Y ( N T ) 〉 = 2 ε D v ∑ i = 1 7 A i α i H 0 ( α i ) f ( α i T ) 与(2)式完全一致。当反应堆临界时，零功率传递函数的最后一个极点 α 7 为零，表达式(2)发散。而将函数 f ( T ) 进行泰勒展开取前两阶 lim x → 0 f ( x ) = 1 2 x − 1 6 x 2 代入有限测量时间模型(4)的零极点项，有 lim α 7 → 0 A 7 α 7 H 0 ( α 7 ) [ f ( α 7 T ) − 1 N f ( α 7 N T ) ] = 1 12 A 7 2 ( N − 1 ) T 2 可以看出反应堆临界时，方差均值比随着测量时间NT的增长而增加，除非测量时间无限，否则不会发散。Wallerbos在瑞士Paul Scherrer学院的HTR-PROTEUS试验装置上进行了试验验证 [ 8 ] 。测量统计时间NT在60 ms至6 s范围内， Y ( N T ) 试验值与理论计算值很接近，反应性的测量值与脉冲中子源法也相当吻合。
system: [SYSTEM]

user: 3. 计数时间间隔效应 两次中子计数之间不可避免存在一个时间间隔 Δ [ 9 ] ，如图2所示，造成计数丢失。 图2. 计数时间间隔示意图 如果 Δ 不可忽略，Yamane如下定义了这种情况下的均值和均方值 c ¯ = ∑ i = 1 N c i ， c 2 ¯ = ∑ i = 1 N c i 2 ¯ + 2 ∑ i < j c i c j ¯ 由此建立了考虑时间间隔 Δ 后计算相关中子测度Y的方法 [ 10 ] Y ( N T ) = 2 ε D v ∑ i = 1 7 A i α i H 0 ( α i ) [ 1 − 1 − e − α i T α i T + ( 1 − e − α i T ) 2 α i T     × ∑ j = 1 N N − j N exp ( − α i ( j Δ + j − 1 ) T ) ] (5) 符号含义同前。当 Δ = 0 时，Yamane验证了式(5)与式(2)保持一致。Wallerbos [ 11 ] 在HTR-PROTEUS试验装置上的试验表明，两次中子计数之间时间间隔 Δ = 2   ms 时， Y ( N T ) 试验值与理论计算值相当吻合。Kitamura [ 12 ] 在京都大学的临界装置KUCA上的试验表明，当 Δ 与每次计数时间T之比从0变化至100%的情况下，利用式(5)都能很好地实现反应性的测量。
system: [SYSTEM]

user: 4. 使用多个探测器的问题 如果同时使用多个探测器进行反应性测量，测量采集的数据如何进行综合呢？Yamane从主方程出发，引入几率母函数，提出了这种情况下计算相关中子测度Y的方法 [ 13 ] Y s y n ( T ) = ( 1 ∑ i = 1 n ε i D ( d i ) ) { ∑ i = 1 n ε i 2 D 2 ( d i ) × [ Q T ( 2 − d i T ) ( − d i T ) + A ( 1 − d i T ) e − α d i f ( T − d i ) ]     + 2 ∑ i = 1 n ∑ j = i + 1 n ε i ε j D ( d i ) D ( d j ) A f ( T ) } (6) n为探测器的个数， A = v ( v − 1 ) ¯ α 2 τ f 2 ， f ( T ) = 1 − 1 − e − α T α T ，Q为裂变率， τ f 为发生一次裂变的平均时间， D ( d i ) 是关于不灵敏时间 d i 的函数，其余符号含义同前。Yamane讨论了n个探测器的效率 ε i 和不灵敏时间 d i 都相同均为 ε 和d的情形，假设不灵敏时间d与测量时间T之比 d / T 很小， f ( T − d ) 几乎等于 f ( T ) ，则(6)简化为 Y s y n ( T ) = − 2 ε D ( d ) Q d + n ε D ( d ) A ( 1 − 1 − e − α T α T ) 不考虑不灵敏时间效应， d = 0 ， D ( d i ) = 1 ，进一步简化为 Y s y n ( T ) = n ε A ( 1 − 1 − e − α T α T ) (7) 与式(2)仍然保持一致。
system: [SYSTEM]

user: 5. 中子源释放相关中子 以上分析的各种情况都是假设使用传统的中子源，如Am-Be，Po-Be等，通过衰变射线在靶核上引起中子发射，每次事件只释放一个中子。每次事件释放的中子相互独立，堆内中子的相关性仅仅来源于同一裂变链上中子的相关性。 随着对加速器驱动次临界系统(Accelerator Driven Sub-critical System ADS)的研究越来越多，研究者尝试将Feynman方法应用于ADS，并发现了新的问题。 加速器驱动次临界系统包括中能质子加速器、外源中子产生器和次临界反应堆。加速器加速的质子轰击重金属靶(如铅)与靶核发生散裂反应，1个质子引起的散裂反应可产生几十个中子。用散裂中子作为外中子源作用于次临界反应堆上。 因此加速器驱动次临界系统中子源的特点是每次释放的中子数不只一个，增加了一个堆内中子数涨落的因素。对于自发裂变中子源 252 Cf也是这样。堆内中子的相关性不仅来源于同一裂变链上中子的相关性，而且来源于中子源同一次事件(ADS的散裂反应或 252 Cf的自发裂变)产生的中子之间的相关性。 中子源Am-Be，Po-Be等，一次事件释放一个中子，中子源释放的中子服从泊松分布。使用 252 Cf作为中子源的反应堆(有的反应堆是使用 252 Cf正比计数管， 252 Cf正比计数管既做中子探测器，又做中子源)， 252 Cf一次自发裂变释放1~8个中子，中子源释放的中子服从复合泊松分布。由加速器驱动的中子源，1个质子引起的散裂反应可产生几十个中子，中子源释放的中子同样服从复合泊松分布。 假设 p ( t 1 ， t 2 ) 是在 ( t 1 ， t 1 + d t 1 ) 和 ( t 2 ， t 2 + d t 2 ) 时间间隔内都探测到一个中子的几率，有 p ( t 1 ， t 2 ) = p c o r r ( t 1 ， t 2 ) + p ( t 1 ) p ( t 2 ) = p c ， f ( t 1 ， t 2 ) + p c ， s ( t 1 ， t 2 ) + p ( t 1 ) p (t2) p c ， f ( t 1 ， t 2 ) 表示由同一裂变链引起中子相关的几率， p c ， s ( t 1 ， t 2 ) 表示由中子源引起中子相关的几率， p ( t 1 ) 、 p ( t 2 ) 分别是在 ( t 1 ， t 1 + d t 1 ) 和 ( t 2 ， t 2 + d t 2 ) 时间间隔内各自独立探测到一个中子的几率。在不考虑缓发中子效应时，Pazsit得到了此时计算相关中子测度Y的方法 [ 14 ] Y ( T ) = ε D v 1 ρ 2 ( 1 + δ ) ( 1 − 1 − e − α T α T ) (8) 其中， δ = q ¯ D q v ¯ D v ， D q = q ( q − 1 ) ¯ q ¯ 2 ，q为中子源每次事件释放的中子数，其余符号同前。与(1)式比较，可以看出由于中子源引起的相关中子的存在，使得中子相关部分的测度 Y ( T ) 增大了，这对于反应性测量是有利的。
system: [SYSTEM]

user: 6. 小结 对利用Feynman方法进行反应性测量的半个多世纪来的文献资料分析总结表明，国内外研究人员对Feynman方法做了持续不懈的研究。针对理论分析和实验过程中碰到的一些新问题，努力提出了解决方法。并利用研究装置或采用蒙特卡罗方法检验处理过程中的假设是否合理、方法是否有效。 关于测量不灵敏时间问题的处理存在不小的分歧，不论采用基于Markov过程的方法 [ 15 ] [ 16 ] [ 17 ] [ 18 ] [ 19 ] 还是相关滤波技术的方法 [ 5 ] [ 20 ] [ 21 ] [ 22 ] ，分析过程和结论都十分复杂，这应该是式(6)中的不灵敏时间函数未能显式表示的原因 [ 13 ] 。而且由于计算机技术的快速发展，传统的多道分析仪可以使用高速数据采集卡来代替。
system: [SYSTEM]

user: 摘 要 以广东大亚湾核电站常规岛热力系统为仿真对象，基于SimStore软件包，根据满负荷状态下的质量守恒和能量守恒，运用矩阵法和定功率法，搭建仿真模型，进行典型负荷稳态仿真和70%~100%负荷升负荷仿真。仿真系统在100%、90%和70%负荷下稳定运行时，精度可以接受；在70%~100%的升负荷仿真中，主要设备的流量与压力的变化符合逻辑，为核电站常规岛热力系统模拟运行、仿真与优化，提高相关教育与实践培训的效率和开发海上核动力平台提供参考。 关键词 :常规岛，热力系统，仿真，SimStore Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 核电由于排放少，且能量巨大，是未来清洁能源的主流发展方向之一。然而，核电站的核泄漏事故风险使有关核电的大部分实验与教学无法在真实的核电站中进行。因此，仿真技术在核电以及相关领域起到不可替代的作用。在核电仿真技术领域，国外开展研究的时间相对较早，主要仿真软件有RELAP 5、ADMIRE、vPower等。目前，国内已有多名学者从事核电仿真研究。文献 [ 1 ] [ 2 ] [ 3 ] 运用RELAP 5进行压水堆核电站常规岛热力系统建模；文献 [ 4 ] 运用C++Builder 6.0搭建核电站常规岛仿真平台。 以大亚湾核电站常规岛热力系统为仿真对象，基于SimStore软件包，结合由矩阵法、定功率法和质量守恒定律改进得来的满负荷下二回路热力参数算法与结合范德瓦尔斯方程改进弗留格尔公式得来的非满负荷下二回路热力参数算法，搭建核电站常规岛热力系统仿真模型并在实践中进行改进，对典型负荷(100%、90%和70%负荷)进行稳态仿真和70%~100%负荷升负荷仿真，获得相关热工参数与变化趋势，与实际运行参数进行比对，为核电站常规岛热力系统模拟运行、仿真与优化，提高相关教育与实践培训的效率和开发海上核动力平台提供参考。改进后的模型部分截图如图1所示。 图1. 仿真模型局部图
system: [SYSTEM]

user: 2. 仿真模型热力参数计算方法 2.1. 满负荷下的热力参数计算方法 2.1.1. 效率计算法 任何热力系统或系统中的设备能够有效利用的能量占外界提供能量的比例，定义为系统或设备的效率。效率计算法以热力学第一定律为理论基础，利用系统和设备的热效率对热力系统和设备的能量利用情况进行定量表示 [ 5 ] [ 6 ] 。 2.1.2. 热平衡法及矩阵法 热平衡法可分为串联法与并联法。根据热力学第一定律，对抽汽回热系统中的每一级加热器列出相应热平衡方程，与其它热量与质量守恒方程联立求解各加热器的抽汽量。矩阵法通过建立各级回热加热器的热平衡方程式，求解一个包括各级抽汽量的线性方程组，进而获得热力系统各项参数，其特点是一次能计算大量未知参数，具有程序化、通用化、快速化等特点，是热平衡法的衍生物 [ 7 ] 。以下是矩阵法需要的方程和方程组： 汽轮机装置的功率方程： 当系统的结构、循环参数和抽汽参数确定时，理论上抽汽回热机组1 kg新蒸汽所做的功可表示为 N = ( h 0 − h n ) − ∑ i = 1 m α i ( h e ， i − h n ) (1) 式中： h 0 ， h e ， i ， h n ——新蒸汽焓、抽汽焓和凝汽焓，单位为kJ/kg；m——总抽汽数； α i ——第i级抽汽份额。 定义 h i d 为1 kg抽汽所引起的做功不足，即 h i d = h e ， i − h n (2) 方程的向量形式为 N = ( h 0 − h n ) − [ α i ] T [ h i d ] T (3) 式中： [ α i ] = [ α 1 ， ⋯ ， α i ， ⋯ ， α m ] T (4) [ h i d ] = [ h 1 − h n ， ⋯ ， h i − h n ， ⋯ ， h m − h n ] T (5) 给水加热器热平衡方程式： 如图2所示，系统中的给水加热器存在以下热平衡关系： 图2. 热力系统原理图 α 1 q 1 = τ 1 (6) α 1 γ 2 + α 2 q 2 = τ 2 (7) ( α 1 + α 2 ) γ 3 + α 3 q 3 = τ 3 (8) α 4 q 4 = ( 1 − α 3 − α 2 − α 1 ) τ 4 (9) α 4 γ 5 + α 5 q 5 = ( 1 − α 3 − α 2 − α 1 ) τ 5 (10) α 6 q 6 = ( 1 − α 5 − α 4 − α 3 − α 2 − α 1 ) τ 6 (11) α 7 q 7 = ( 1 − α 5 − α 4 − α 3 − α 2 − α 1 ) τ 7 (12) 联立以上各式，并将其写成标准的矩阵形式，可得 A [ α i ] = [ b i ] (13) 式中矩阵A是热力系统的结构矩阵； [ α i ] 是抽汽系数向量； [ b i ] 是特征向量。其中，矩阵A的每一行代表一个加热器的热平衡方程，每一列代表了每股抽汽所经历的放热过程；向量 [ b i ] 表示给水从凝汽器开 始经过多级加热器的升温过程。 由于需要预处理的参数较多，逐个设备计算耗时较长。在将源于热平衡法的矩阵法与源于效率计算法的定功率法适当结合的基础上，结合质量守恒定律，利用MATLAB进行编程，一次性快速求解满负荷下的常规岛仿真模型的热力参数。 2.2. 非满负荷下的热力参数计算方法 针对满负荷情况下核电站常规岛热力参数的计算方法不适用于非满负荷下的情况。在任何负荷下，核电站发电与其汽轮机所处工况密切相关，弗留格尔公式是汽轮机级组在变工况状态下的平衡关系式，如下所示： G 2 G 1 = p 2 1 2 − p 2 2 2 p 1 1 2 − p 1 2 2 T 0 2 T 0 1 (14) 上式中， G 1 、 G 2 分别为变工况前后机组通流量，单位为kg∙s −1 ； T 0 1 、 T 0 2 分别为变工况前后级组前蒸汽温度，单位为K； p 2 1 、 p 2 2 分别为变工况后的级前与级后压力，单位为Pa； p 1 1 、 p 1 2 分别为变工况前的 级前与级后压力，单位为Pa。 然而，该公式没有考虑到核电站汽轮机通流介质不是理想气体，导致公式精度差，需要结合范德瓦尔斯方程改造弗留格尔公式，获得筛选仿真数据的公式以避免过高或过低数据影响分析，如下所示： f ( p ， G ， T ， x ) i = G 2 i G 1 i ( p 2 1 + a p 2 1 v m 2 1 2 ) 2 − ( p 2 2 + a 2 v m 2 2 2 ) 2 ( p 1 1 + a p 1 1 v m 1 1 2 ) 2 − ( p 1 2 + a p 1 2 v m 1 2 2 ) 2 T 0 2 T 0 1 (15) f c = | f ( p ， G ， T ， x ) a f ( p ， G ， T ， x ) b − 1 | (16) 其中：x——湿蒸汽状态干度；a p ——相应状态范式参数； v m ——相应状态比体积；i、a、b、c——顺序与位置； f ( p ， G ， T ， x ) ——契合度； f c ——偏离度。 在仿真模型调试完毕后，先进行不同负荷下的稳态仿真。在利用Excel插件，结合上式计算得到每 组仿真数据的契合度 f ( p ， G ， T ， x ) i 与每个数据之间的偏离度 f c 后，通过MATLAB循环程序，建立 f c 的行向量或列向量，得到概率不低于90%的情况下，每组数据中汽轮机级组热力参数的偏离度小于设定值 θ 的相应时刻的测得的全部数据，求出均值作为仿真数据。经过实践，当设定 θ = 1.0 × 10 − 2 时，可保证达到要求。