Risk Evaluation Model of Blockchain Integration Based on CKKS Encryption Scheme
-
摘要:
针对区块链平台中存在的多方数据交互不可信以及隐私数据易泄露等问题,基于CKKS (Cheon-Kim-Kim-Song)全同态加密方案,提出了一种集成风险评价模型,把同态加密算法应用到风险评价中,将多种评价模型与同态加密结合起来. 首先,利用三角模糊综合评价方法确定各评价指标的权重,通过多种评价方法处理分布式数据库中的样本数据,获得相关节点对同一交易事件的风险评价结果;其次,利用公钥对评价结果进行加密并进行同态运算,获得密文综合评价结果,以避免风险评价过程中的数据泄露;再次,利用私钥对评价结果进行解码,获得明文综合评价结果;最后,选取5 000个中欧班列企业的样本数据作为案例,利用决策树模型、Adaboost模型、Bagging模型、ExtraTree极端随机数模型、GBDT (gradient boosting regression trees)模型、KNN(K-nearest neighbor)模型、随机森林模型、SVM (support vector machine)模型等最为常见的评价模型进行风险评价,并将经CKKS方案加密后的综合评价结果与明文直接计算的综合评价结果和经BFV (Brakerski-Fan-vercauteren)方案加密后的综合评价结果进行了对比. 结果表明:该集成风险评价模型具有普适性,对较为常见的评价模型均能适用;模型的综合评价结果误差率较小,与实际结果的误差率均在10−9以内;与BFV方案加密后的结果相比,经CKKS方案加密后的结果误差率小于前者的十万分之一,评价结果更为准确.
Abstract:Aiming at solving problems of distrust in the data interaction between multiple parties and leakage of privacy data, an integrated risk assessment model based on the Cheon-Kim-Kim-Song (CKKS) fully homomorphic encryption scheme is proposed by using homomorphic encryption algorithm in combination with multiple assessment models for risk assessment. Firstly, the triangular fuzzy comprehensive evaluation method is used to determine weights of evaluation indexs, and a variety of evaluation methods are used to process the sample data in the distributed database to obtain the risk evaluation results of related nodes for the same transaction event. Secondly, the public key is used to encrypt the evaluation results, and then homomorphic operation is performed to obtain the encrypted comprehensive evaluation results, so as to avoid data leakage in the process of risk assessment. Thirdly, the private key is used to decode the evaluation results to obtain decrypted comprehensive evaluation results. Finally, 5000 samples from China Railway Express companies are taken as cases to evaluate the risk with eight common evaluation models, including the decision tree model, Adaboost model, Bagging model, ExtraTree model, gradient boosting regression tree (GBDT) model, K-nearest neighbor (KNN) model, random forest model, and support vector machine (SVM) model. In addition, the comprehensive evaluation results encrypted by CKKS scheme are compared with the results directly calculated by plaintext and encrypted by BFV scheme. The findings show that: the integrated risk assessment model is universal and applicable to more common assessment models; the error rates of the comprehensive evaluation results obtained by the integrated risk assessment model are small, within 10−9 against the actual results; compared with the result encrypted by BFV scheme, the error rate of the result encrypted by CKKS scheme is less than 1/100000 of the former, and the evaluation result is more accurate.
-
Key words:
- Blockchain /
- data leakage /
- fully homomorphic encryption scheme /
- risk assessment /
- integrated model
-
车体作为地铁车辆关键承载零部件,在运行过程中承受来自转向架、车钩、车下设备以及乘客质量等一系列复杂载荷. 新造地铁车体由于轻量化需求,多采用铝合金型材整体焊接而成,一节车厢焊缝多达数千条. 以传统A型地铁为例,车体侧墙两边各开有5门4窗,且门窗尺寸远大于干线车体,大大增加了车身强度薄弱位置的数量.基于以上特点,地铁车体在承受复杂载荷时的疲劳问题越来越突出.
传统A型地铁铝合金车体长约为22 m、宽度约为3 m、高度约为3 m,整备质量约为25 t,超员载荷约为50 t. 整车的疲劳试验试验周期长、经费高、风险大. 在ANSYS软件中,子模型方法通过对局部关注区域的网格进行细化,并在区域边界处施加等效的边界条件,使得计算结果更精确,计算速度更快. 故可考虑截取车体受力最恶劣的部分,对该部分的边界施加等效边界条件后进行疲劳试验. 根据车体的受力状态以及运营维修经验,车体疲劳裂纹往往产生于牵枕缓部位[1],因此可以截取车体底架端部包含牵枕缓部位的部分进行子模型疲劳试验.
胡杰鑫等[2]利用Monte Carlo方法,对铝合金地铁车体枕梁进行了虚拟疲劳试验,并与实际疲劳试验结果进行对比,验证了该方法的有效性;王玉伟[3]基于实测载荷编制CRH2型车枕梁关键部位的应力谱,研究了运营里程与损伤之间的关系;李珊珊[4]利用子模型方法,对CRH2型车枕梁结构网格进行了细分,提出了枕梁的补强方案;臧伟锋等[5]针对飞机机身壁板设计了一种新型工装,该工装使得机身壁板的应力分布与实际全尺寸机身的应力分布相一致,并进行了内压载荷疲劳试验;四方股份联合航空工业飞机强度研究所对CRH380中间车车体全长的1/3结构进行了气密强度疲劳试验[6];Gutiérrez-Carvajal等[7]设置了多种载荷工况对铝合金车体枕梁进行了疲劳试验,利用累积损伤法评估了枕梁的疲劳寿命.
本文将设定边界条件下的静强度试验结果与有限元计算结果进行对比分析,验证了端部底架子模型(以下简称端部底架)有限元模型的正确性. 研究了端部底架边界条件模拟方法和试验加载方法,根据标准EN12663-1:2010 + A1:2014[8](以下简称EN 12663标准)规定的载荷进行了疲劳损伤计算,得到了端部底架与全尺寸车体关键位置损伤分布的规律,验证了以端部底架子模型疲劳试验代替全尺寸车体疲劳试验的可行性. 通过端部底架疲劳试验结果对全尺寸车体的疲劳寿命进行了评估.
1. 端部底架边界条件及载荷
1.1 边界条件
在全尺寸车体模型中,车体的底架边梁与车体的侧墙连接,底架的端梁与车体的端墙连接,底架的尾部与车体地板及底架边梁连接. 截取下来的端部底架长5 m,与端墙、侧墙、地板以及边梁的连接均被截断. 所设置的边界条件需保证端部底架关键部位的应力分布与应力大小与全尺寸车体相一致.
在有限元计算时,容易获取端部底架边界节点的力和弯矩,该力和弯矩与施加于端部底架的外部载荷使得底架处于力平衡的状态. 在疲劳试验过程中,底架车钩、空簧和中心销处的集中力可以通过作动器施加,由于底架边界在全尺寸车体中变形较小,边界处的分布力则以约束产生的反力施加.
根据端部底架与全尺寸车体的连接方式,在端部底架的头部两侧与车体侧墙连接处、枕梁两侧与车体侧墙连接处约束端部底架的垂向自由度. 在尾部地板处约束底架的纵向和垂向自由度. 横向约束在运营过程中由车体的惯性产生,在台架试验中无法复现,因此选取与牵枕缓焊缝应力无关的边梁侧面进行约束,如图1所示,图中:X轴所示方向为纵向;Y轴所示方向为横向;Z轴所示方向为垂向.
1.2 试验载荷计算
根据EN 12663标准中关于车体疲劳载荷的规定:地铁车体应根据定员载荷(m1 + m3)和超员载荷(m1 + m4)分别设计疲劳载荷工况,其中:m1为车体整备质量;m3为定员乘客质量;m4为超员乘客质量,并按其在服役寿命期间所占的比例进行疲劳试验. 其他类型轨道交通车辆车体根据其载荷特点,选取恒幅定员或恒幅超员疲劳载荷进行疲劳试验.
由于地铁车体定员载荷和超员载荷在服役期间所占比重难以准确统计,为得到较为保守的结果,认为车体在服役过程中一直承受最大载荷,车体的疲劳寿命最短,即以恒幅超员载荷(m1 + m4)进行数值计算和疲劳试验,m1 = 28.279 t,m4 = 22.440 t.
EN 12663标准中规定整车纵向和横向疲劳载荷为 ± 0.15g,垂向载荷为g ± 0.15g,但均为单个载荷分别作用1000万次的疲劳载荷. 对于组合载荷,EN 12663标准第6.8节推荐引用VDV 152-2016[9]标准(以下简称VDV 152标准).
VDV 152标准第5.3.2节中将X、Y、Z方向的载荷进行组合,建立了3个载荷工况,每个工况循环200万次,见表1.
表 1 标准之间疲劳载荷工况对比Table 1. Comparison of fatigue load cases between standards标准 工况 X Y Z VDV 152 1 ± 1.20 m/s2 g ± 1.90 m/s2 2 ± 1.60 m/s2 g ± 1.00 m/s2 3 ± 1.50 m/s2 ± 0.60 m/s2 g ± 1.00 m/s2 EN 12663 1 ± 0.15g 2 ± 0.15g 3 ± 0.15g VDV 152标准中载荷工况的组合与EN 12663标准中规定的3个方向载荷每个方向单独作用1000万次是等效的,即在疲劳试验中,可用VDV 152标准规定的疲劳载荷代替EN 12663标准中规定疲劳载荷[9].
对于X方向的载荷,VDV 152标准只在一个工况下循环了200万次,因此需根据等损伤原则,将X方向200万次载荷折算成EN 12663标准中与垂向、横向共同作用的1000万次载荷.
1.2.1 车钩载荷
车钩载荷主要源于列车牵引或制动时的不同步,对于A型地铁车辆,一般采用4动2拖编组,拖车为首尾车,最大车钩载荷发生在拖车和动车连接的车钩位置. 按启动加速度a = 1 m/s2计算,最大车钩载荷为
FXmax=4a(m1+m4)2=101.44kN. (1) 根据VDV 152标准第5.3.2.1节,最大车钩力仅在牵引(制动)不同步的瞬间产生,对于动力分散列车,车钩力的疲劳载荷为最大车钩力的25%. 因此,作用200万次车钩载荷的动态幅值为
${F_{X {\max} }}/{4}=25\; {\rm{KN}}$ .根据疲劳寿命-应力曲线幂函数表达式:
σmN=C. (2) 式中:σ为应力;m、C为材料常数,VDV 152标准第5.3.3.3节推荐m取4;N为循环次数.
按VDV 152标准,将0.15g作用200万次,根据式(2)折算成1000万次,则车钩载荷动态幅值
$F_{\text {coupler }}= 16.96\; \mathrm{kN}.$ 1.2.2 中心销牵引(制动)载荷
依据EN 12663标准按 ± 0.15g施加,中心销牵引(制动)载荷如式(3).
FpivotX=0.15g(m1+m4)2=37.3kN. (3) 按式(3)将200万次载荷折算成1000万次,载荷
$F_{{\rm{pivot}}X}=24.94\; \mathrm{kN}.$ 中心销牵引(制动)载荷的加载高度根据实车确定,本文中加载高度距枕梁下表面为586.5 mm.
1.2.3 中心销横向载荷
列车在曲线运行时会承受来自于转向架横向止挡的横向力FY,依据EN 12663标准按± 0.15g施加,如式(4).
FY=0.15g(m1+m4)2=37.3kN. (4) 该载荷由中心销和空气弹簧共同承担,在试验时仅施加中心销承受的横向载荷. 中心销横向载荷加载高度距枕梁下表面为297 mm. 根据空气弹簧参数,每一个转向架单个空气弹簧承担的横向载荷
${F_{{\rm{airsring}}Y}} = 3.78\;{\rm{kN}}.$ 因此中心销承担的横向载荷为
FpivotY=FY−2FairsringY=29.74kN. (5) 1.2.4 空簧垂向载荷
根据表1超员载荷确定空簧承受的垂向载荷. 依据EN 12663标准,按g± 0.15g施加,由空气弹簧承担.
空簧载荷均值为
FZm=g(m1+m4)4=124.39kN. (6) 空簧载荷幅值为
FZa=0.15g(m1+m4)4=18.66kN. (7) 2. 端部底架及整车有限元建模分析
2.1 端部底架及整车有限元建模
采用Hypermesh对端部底架进行网格划分,网格类型为Shell181. 底架所用材料为铝合金,材料密度为2.7×10−9 t/mm3,弹性模量为69 GPa,泊松比为0.3. 中心销所用材料为钢,材料密度为7.85×10−9 t/mm3,弹性模量为210 GPa,泊松比为0.3. 中心销与枕梁下表面通过梁单元(Beam188)模拟的螺栓连接. 端部底架有限元模型共划分单元213586个,节点197352个. 底架的约束方式如图2所示,FZ为Z方向的载荷.通过刚性元在车钩座、中心销和空簧位置施加集中力载荷.
全尺寸车体建模方式与端部底架相同,共划分单元1956597个,节点1779596个,如图3所示.
车体在4个空簧位置约束垂向,在中心销横向止挡位置约束横向. 通过刚性元在车钩座、中心销空簧位置施加纵向集中力载荷,垂向载荷及横向载荷以重力场形式施加. 乘客质量通过柔性元均布在地板上,空调及车下设备等考虑其对整车刚度的加强,通过刚性元连接到设备安装座,车门及车窗通过刚性元连接到安装位置.
2.2 模型正确性验证
试验过程中,利用应变片对端部底架关键部位进行应力测量,测点布置见图4,S1~S9为应变片测点,R1~R4为应变花测点,应变花用相当应力进行评定.
当载荷加到最大时,读取各个位置的应力值如图5所示.
由图5可知:测点S1、S2、S3、S4、S5、S9、R1、R2、R3、R4的实测应力与计算应力相对误差较小,测点S3误差最小为0.64%;测点S6、S7、S8误差较大,其中测点S8误差最大为122.2%.
测点S8位于枕内加强筋的焊趾处,该处应力梯度较大,且打磨后的几何轮廓不平整,难以准确确定焊趾位置,因此应变片位置的细微差别会导致测量结果较大的变化;测点S6与测点S3为对称点、测点S7与测点S1为对称点,其仿真应力值误差分别为8.2%和12.2%,且测点S1和测点S3的实测应力值和仿真应力值相差较小,因此测点S6和测点S7误差较大的原因可能是试验中的约束条件没有仿真中的约束条件理想,约束工装本身与车体之间存在间隙导致.
本文中端部底架所需评估的部位主要位于枕梁内侧加强框与周边零部件的焊缝,无法直接粘贴应变片. 测点S6、S7、S8与所评估部位无直接联系,且测点部位在实际运营过程中状态良好. 除测点S6、S7、S8外,其余测点实测应力与计算应力误差较小,因此,认为端部底架有限元模型能够真实反映实际模型的应力状态.
3. 端部底架及整车疲劳寿命分析
3.1 疲劳理论
基于Miner线性累积损伤理论进行疲劳计算,如式(8).
D=∑niNi, (8) 式中:D为损伤值,当D<1时,认为材料没有疲劳失效的风险.;ni为第i级应力谱下循环的次数;Ni为第i级应力谱下的疲劳寿命.
在ANSYS下计算单位静载荷下焊缝单元的名义应力,在nCode Designlife中施加恒幅疲劳载荷谱.
端部底架同时承受来自车钩座、中心销横向、中心销纵向和空簧垂向载荷,在试验过程中其单元的主应力方向不断发生变化,属于典型的多轴疲劳. 因此采用带绝对值的最大主应力进行疲劳分析.
在疲劳试验中,端部底架承受空簧载荷均值产生的应力,需进行平均应力修正,本文中选取FKM指导书[10]中规定的方法进行平均应力修正. 列车车体在实际服役过程中,结构承载的载荷循环次数远远超过材料疲劳极限处的循环次数,大量小幅值载荷仍能对车体造成损伤,因此将S-N曲线的第二段延伸至应力幅为0处,考虑小载荷循环的影响[11].
3.2 S-N曲线
本文根据Eurocode 9[12]标准建立端部底架铝合金材料不同焊缝类型的S-N曲线.
$\Delta \sigma $ 代表循环应力范围. 根据端部底架焊缝类型选用Eurocode 9中7.2.3号焊缝、7.4.3号焊缝、9.4号焊缝和9.6号焊缝. 焊缝类型如图6所示.各焊缝对应的S-N曲线如图7所示. 由图7可知:各焊缝S-N曲线的拐点位于500万次循环处;S-N曲线的第一段斜率k1 = −1/3.4,第二段斜率k2 = −1/5.4.
3.3 疲劳计算结果
根据有限元模型与图6焊缝类型的对应关系,在nCode中将图7各焊缝的S-N曲线分别赋予给对应的焊缝单元,按1.2节计算的载荷进行疲劳计算. 得到端部底架和全尺寸车体的损伤云图如图8,对比如表2.
表 2 端部底架与全尺寸车体损伤对比Table 2. Damage comparison of the end underframe and the full-scale carbody端部底架损伤 整车损伤 损伤部位 0.5028 0.1910 地板下表面与枕内加强筋底板焊缝直角拐角处 0.1796 0.0962 地板下表面与枕内加强筋底板135°拐角附近区域 0.0857 0.0591 外侧枕内加强筋与枕梁下盖板焊缝 从表2可以看出:端部底架损伤最大的3个部位与全尺寸车体损伤最大的3个部位位置一致,且同一位置端部底架的损伤均大于全尺寸车体损伤. 因此通过端部底架疲劳试验验证整车的疲劳强度可以获得较为保守的结果.
3.4 端部底架疲劳试验
端部底架包括了车体运营过程中受力最恶劣的牵枕缓结构,采用接口工装和MTS液压作动器进行疲劳加载,如图9.
端部底架的垂向和横向运动由头部及枕梁两侧边梁的“C”字形工装提供约束,纵向运动通过固定尾部地板面提供约束,如图10所示. 按1.2节计算的载荷进行1000万次循环的疲劳试验,对图8(a)所示的区域进行渗透探伤,未发现裂纹,因此由端部底架的疲劳试验结果认为整车的疲劳寿命满足EN 12663标准的要求.
4. 结 论
1) 全尺寸车体的疲劳试验对试验设备、试验成本均有较高要求,用小尺寸子模型代替整车进行疲劳试验能够有效地缩短研发周期,降低试验成本. 本文采用端部底架子模型代替全尺寸车体的疲劳损伤评估结果偏于保守,从数值仿真和台架试验的角度验证了方法对车体疲劳寿命进行评估是可行的.
2) 使用子模型进行疲劳试验时,需保证子模型的边界远离可能的裂纹潜在部位,且保证子模型结构潜在危险点与全尺寸车体结构保持一致. 边界条件的设置可根据边界处刚度的方向特性,选取刚度最大的方向进行约束,释放他方向的约束. 通过有限元仿真设置的边界条件与试验进行对比,获得正确的有限元数值分析模型,为车体的寿命评估提供依据.
3) 结合名义应力法S-N曲线和线性损伤理论,分别计算了端部底架与全尺寸车体的疲劳损伤. 发现该方法依赖于网格质量和S-N曲线的选取,针对车体和其焊缝特点,在计算中均进行了考虑,因此采用的方法适合车体的寿命评估.
4) 采用端部底架子模型的疲劳试验与数值计算,并与全尺寸车体结构的对比分析,表明该车体能够满足EN 12663标准规定的疲劳寿命要求.
致谢:铝合金车体疲劳寿命评估及试验方法研究(2682019CX45).
-
表 1 0.1~0.9标度的含义
Table 1. Meaning of scale 0.1~0.9
标度 对应的三角模糊数 含义 0.1 (0.1,0.1,0.2) 指标 i 相对于指标 j 极端不重要 0.3 (0.2,0.3,0.4) 指标 i 相对于指标 j 明显不重要 0.5 (0.4,0.5,0.6) 指标 i 与指标 j 同样重要 0.7 (0.6,0.7,0.8) 指标 i 相对于指标 j 明显重要 0.9 (0.8,0.9,0.9) 指标 i 相对于指标 j 极端重要 表 2 定性指标量化评分表
Table 2. Quantitative scoring of qualitative indices
定性指标 评价指标分档/分 [8, 10] [4, 8) [0, 4) 交易履约情况(X8) 好 中 差 核心企业的对外担保状况(X9) 几乎无 少量 较多 供应链关系的强度(X10) 高 中 低 表 3 三角模糊打分表
Table 3. Triangular fuzzy scoring results
k B1 B2 B3 B4 B5 1 (0.1,0.2,0.3) (0,0.1,0.3) (0.1,0.3,0.5) (0.5,0.2,0.3) (0.6,0.8,0.9) 2 (0.2,0.3,0.3) (0.4,0.7,0.8) (0.3,0.2,0.2) (0.7,0.8,0.9) (0.1,0.2,0.4) 3 (0.4,0.5,0.3) (0.1,0.3,0.5) (04,0.5,0.7) (0.1,0.4,0.6) (0.1,0.6,0.7) 4 (0.5,0.7,0.3) (0.3,0.6,0.9) (0.4,0.6,0.8) (0.3,0.5,0.7) (0.1,0.3,0.4) 5 (0.6,0.4,0.3) (0.2,0.4,0.5) (0.5,0.8,0.9) (0.5,0.6,0.8) (0.2,0.4,0.5) wm (0.15, 0.21,0.32) (0.13,0.22,0.35) (0.11,0.18,0.28) (0.18,0.25,0.39) (0.09,0.14,0.24) 表 4 指标赋权结果
Table 4. Index weighting results
指标名称 变量 权重 指标名称 变量 权重 企业规模 X1 0.0246 核心企业对外担保情况 X9 0.1512 销售利润率 X2 0.0833 供应链关系强度 X10 0.0014 速动比率 X3 0.0255 赊销周期 X11 0.1146 存货周转率 X4 0.0117 产品可替代性 X12 0.0884 资产负债率 X5 0.0140 权益乘数 X13 0.0176 交易量 X6 0.0748 质押物变现能力 X14 0.0799 交易金额 X7 0.0363 总资产周转率 X15 0.1329 交易履约情况 X8 0.1437 表 5 CKKS方案和BFV方案同态加密结果及误差率
Table 5. Homomorphic encryption results and error rates of CKKS and BFV schemes
编号 KNN 模型 +
随机森林模型CKKS 全同态加密
评价结果CKKS 方案
误差率BFV 全同态加密
评价结果BFV 方案
误差率1 48.6117737805 48.6117737789 3.34 × 10−11 48.6094744655 4.73 × 10−5 2 38.4173642857 38.4173642916 1.53 × 10−10 38.3934853887 6.22 × 10−4 3 66.6095043478 66.6095043506 4.09 × 10−11 66.6003539367 1.37 × 10−4 4 40.7277690992 40.7277690999 1.66 × 10−11 40.7098999864 4.39 × 10−4 5 36.0318162194 36.0318162184 2.58 × 10−11 36.0170565909 4.10 × 10−4 6 39.2235467290 39.2235467305 3.98 × 10−11 39.1413627605 2.10 × 10−3 7 45.2551751634 45.2551751619 3.20 × 10−11 45.2384881180 3.69 × 10−4 8 85.5722701987 85.5722701993 7.54 × 10−12 85.5655826327 7.82 × 10−5 9 59.7378161329 59.7378161324 8.14 × 10−12 59.7254652272 2.07 × 10−4 10 49.4582034252 49.4582034281 5.74 × 10−11 49.4489389677 1.87 × 10−4 11 50.6262683040 50.6262683016 4.62 × 10−11 50.6237741777 4.93 × 10−5 12 59.0855247219 59.0855247221 4.30 × 10−12 59.0248410817 1.03 × 10−3 13 50.8218075188 50.8218075195 1.39 × 10−11 50.8154754434 1.25 × 10−4 14 51.5208750000 51.5208750014 2.77 × 10−11 51.5204698244 7.86 × 10−6 15 40.6476813842 40.6476813860 4.43 × 10−11 40.6440481786 8.94 × 10−5 16 41.8885428177 41.8885428182 1.16 × 10−11 41.8792671831 2.21 × 10−4 17 37.0070685441 37.0070685404 1.01 × 10−10 36.9933666065 3.70 × 10−4 18 66.6304171240 66.6304171283 6.33 × 10−11 66.6209557226 1.42 × 10−4 19 14.8742302481 14.8742302472 6.03 × 10−11 14.8490706101 1.69 × 10−3 20 38.9383642857 38.9383642857 3.34 × 10−11 38.8372207385 2.60 × 10−5 -
[1] ELGHAISH F A K, ABRISHAMI S, HOSSEINI M R. Integrated project delivery with blockchain: an automated financial system[J]. Automation in Construction, 2020, 114(1): 209-224. [2] LIU H, ZHANG Y, YANG T. Blockchain-enabled security in electric vehicles cloud and edge computing[J]. IEEE Network, 2018, 32(3): 78-83. doi: 10.1109/MNET.2018.1700344 [3] 邰雪,孙宏斌,郭庆来. 能源互联网中基于区块链的电力交易和阻塞管理方法[J]. 电网技术,2016,40(12): 3630-3638.TAI Xue, SUN Hongbin, GUO Qinglai. Electricity transactions and congestion management based on blockchain in energy internet[J]. Power System Technology, 2016, 40(12): 3630-3638. [4] 张宁,王毅,康重庆,等. 能源互联网中的区块链技术:研究框架与典型应用初探[J]. 中国电机工程学报,2016,36(15): 4011-4023.ZHANG Ning, WANG Yi, KANG Chongqing, et al. Blockchain technique in the energy Internet: preliminary research framework and typical applications[J]. Proceedings of the CSEE, 2016, 36(15): 4011-4023. [5] 袁勇,王飞跃. 区块链技术发展现状与展望[J]. 自动化学报,2016,42(4): 481-494.YUAN Yong, WANG Feiyue. Blockchain: the state of the art and future trends[J]. Acta Automatica Sinica, 2016, 42(4): 481-494. [6] OMRAN Y, HENKE M, HEINES R, et al. Blockchain-driven supply chain finance: towards a conceptual framework from a buyer perspective[C]//2017: 26th Annual Conference of the International Purchasing and Supply Education and Research Association. Budapest: [s.n.], 2017: 1-15. [7] 徐忠,邹传伟. 区块链能做什么、不能做什么?[J]. 金融研究,2018(11): 1-16.XU Zhong ZOU Chuanwei. What can blockchain do and cannot do?[J]. Journal of Financial Research, 2018(11): 1-16. [8] 龚强,班铭媛,张一林. 区块链、企业数字化与供应链金融创新[J]. 管理世界,2021,37(2): 3,22-34. doi: 10.3969/j.issn.1002-5502.2021.02.004GONG Qiang, BAN Mingyuan, ZHANG Yilin. Blockchain, enterprise digitalization and supply chain finance innovation[J]. Journal of Management World, 2021, 37(2): 3,22-34. doi: 10.3969/j.issn.1002-5502.2021.02.004 [9] RIVEST RL, ADLEMAN L M, DERTOUZOS M L. On databanks and privacy homomorphisms[J]. Foundations of Secure Computation, 1978, 76(4): 169-179. [10] GOLDWASSER S, MICALI S. Probabilistic encryption[J]. Journal of Computer and System Sciences, 1984, 28(2): 270-299. doi: 10.1016/0022-0000(84)90070-9 [11] GENTRY C. Fully homomorphic encryption using ideal lattices[C]//In Proceedings of the Forty-First Annual ACM Symposium on Theory of computing, Association for Computing Machinery. New York: ACM Press, 2009: 169-178. [12] GENTRY C. A fully homomorphic encryption scheme[M]. Ann Arbor: [s.n.], 2009. [13] CHEON J H , KIM A , KIM M , et al. Homomorphic encryption for arithmetic of approximate numbers[C]//International Conference on the Theory and Application of Cryptology and Information Security. Hong Kong: Springer, 2017: 409-437. [14] 郑尚文,刘尧,周潭平,等. 优化的基于错误学习问题的CKKS方案[J]. 计算机应用,2021,41(6): 1723-1728.ZHENG Shangwen, LIU Yao, ZHOU Tanping, et al. Optimized CKKS scheme based on learning with errors problem[J]. Journal of Computer Applications, 2021, 41(6): 1723-1728. [15] NAKASUMI M. Information sharing for supply chain management based on block chain technology[C]//2017 IEEE 19th Conference on Business Informatics. Thessaloniki: IEEE, 140-149. [16] 刘彦松,夏琦,李柱,等. 基于区块链的链上数据安全共享体系研究[J]. 大数据,2020,6(5): 92-105.LIU Yansong, XIA Qi, LI Zhu, et al. Research on secure data sharing system based on blockchain[J]. Big Data Research, 2020, 6(5): 92-105. [17] 钱萍,吴蒙. 同态加密隐私保护数据挖掘方法综述[J]. 计算机应用研究,2011,28(5): 1614-1617,1622. doi: 10.3969/j.issn.1001-3695.2011.05.004QIAN Ping, WU Meng. Survey of privacy preserving data mining methods based on homomorphic encryption[J]. Application Research of Computers, 2011, 28(5): 1614-1617,1622. doi: 10.3969/j.issn.1001-3695.2011.05.004 [18] DU M X, CHEN Q J, XIAO J, et al. Supply chain finance innovation using blockchain[J]. IEEE Transactions on Engineering Management, 2020, 67(4): 1045-1058. doi: 10.1109/TEM.2020.2971858 [19] 王竹泉,宋晓缤,王苑琢. 我国实体经济短期金融风险的评价与研判——存量与流量兼顾的短期财务风险综合评估与预警[J]. 管理世界,2020,36(10): 156-170,216. doi: 10.3969/j.issn.1002-5502.2020.10.012WANG Zhuquan, SONG Xiaobin, WANG Yuanzhuo. Objective evaluation and rational judgment of short-term financial risk in China’s real economy: comprehensive assessment and early warning of short-term financial risk considering stock and flow[J]. Management World, 2020, 36(10): 156-170,216. doi: 10.3969/j.issn.1002-5502.2020.10.012 [20] 廖礼坤,张炜. 导入期风险企业的风险综合评价[J]. 西南交通大学学报,2004,39(5): 590-594. doi: 10.3969/j.issn.0258-2724.2004.05.008LIAO Likun, ZHANG Wei. Synthetical evaluation of risks of venture business during start-up period[J]. Journal of Southwest Jiaotong University, 2004, 39(5): 590-594. doi: 10.3969/j.issn.0258-2724.2004.05.008 [21] 龙云飞. 基于熵值法的中小企业供应链融资信用风险评价[J]. 统计与决策,2013(13): 177-179. [22] 程昔武,丁忠明. 高校负债融资风险及其评价方法研究[J]. 财贸研究,2009,20(6): 131-138,152. doi: 10.3969/j.issn.1001-6260.2009.06.022CHENG Xiwu, DING Zhongming. Research on risk of debt financing to colleges & universities and it’s evaluating[J]. Finance and Trade Research, 2009, 20(6): 131-138,152. doi: 10.3969/j.issn.1001-6260.2009.06.022 [23] 曹清玮,戴丽芳,孙琪,等. 社会网络环境下基于分布式信任的在线评价方法[J]. 控制与决策,2020,35(7): 1697-1702.CAO Qingwei, DAI Lifang, SUN Qi, et al. A distributed trust based online evaluation under social network[J]. Control and Decision, 2020, 35(7): 1697-1702. [24] 仇文革,李俊松,胡兰,等. 基于WebGIS的地下工程安全风险管理系统[J]. 西南交通大学学报,2011,46(6): 953-959,965. doi: 10.3969/j.issn.0258-2724.2011.06.011QIU Wenge, LI Junsong, HU Lan, et al. WebGIS-based safety risk management system of underground engineering[J]. Journal of Southwest Jiaotong University, 2011, 46(6): 953-959,965. doi: 10.3969/j.issn.0258-2724.2011.06.011 [25] 胡兰,胡培. 基于概率论-逻辑学的隧道各方关系与风险研究[J]. 西南交通大学学报,2013,48(6): 1122-1128. doi: 10.3969/j.issn.0258-2724.2013.06.024HU Lan, HU Pei. Relation and risk in owner and contractors of tunnel projects based on probability and logic theories[J]. Journal of Southwest Jiaotong University, 2013, 48(6): 1122-1128. doi: 10.3969/j.issn.0258-2724.2013.06.024 [26] 何正友,冯玎,林圣,等. 高速铁路牵引供电系统安全风险评估研究综述[J]. 西南交通大学学报,2016,51(3): 418-429. doi: 10.3969/j.issn.0258-2724.2016.03.002HE Zhengyou, FENG Ding, LIN Sheng, et al. Research on security risk assessment for traction power supply system of high-speed railway[J]. Journal of Southwest Jiaotong University, 2016, 51(3): 418-429. doi: 10.3969/j.issn.0258-2724.2016.03.002 [27] 吴波前,蔡伯根,陆德彪,等. 基于GNSS/INS的列车定位风险评估方法[J]. 西南交通大学学报,2020,55(6): 1191-1198. doi: 10.3969/j.issn.0258-2724.20190981WU Boqian, CAI Baigen, LU Debiao, et al. GNSS/INS based risk assessment in train localization[J]. Journal of Southwest Jiaotong University, 2020, 55(6): 1191-1198. doi: 10.3969/j.issn.0258-2724.20190981 [28] 刘曙阳. CI 系统开发技术[M]. 北京: 国防工业出版社, 1997. [29] 刘颖,张丽娟,韩亚男,等. 基于粒子群协同优化算法的供应链金融信用风险评价模型[J]. 吉林大学学报(理学版),2018,56(1): 119-125.LIU Ying, ZHANG Lijuan, HAN Yanan, et al. Financial credit risk evaluation model of supply chain finance based on particle swarm cooperative optimization algorithm[J]. Journal of Jilin University (Science Edition), 2018, 56(1): 119-125. -