An Algorithm for Short-Circuit Current Interval in Distribution Networks with Inverter Type Distributed Generation Based on Affine Arithmetic

During faults in a distribution network,the output power of a distributed generation(DG)may be uncertain.Moreover,the output currents of distributed power sources are also affected by the output power,resulting in uncertainties in the calculation of the short-circuit current at the time of a fault.Additionally,the impacts of such uncertainties around short-circuit currents will increase with the increase of distributed power sources.Thus,it is very important to develop a method for calculating the short-circuit current while considering the uncertainties in a distribution network.In this study,an affine arithmetic algorithm for calculating short-circuit current intervals in distribution networks with distributed power sources while considering power fluctuations is presented.The proposed algorithm includes two stages.In the first stage,normal operations are considered to establish a conservative interval affine optimization model of injection currents in distributed power sources.Constrained by the fluctuation range of distributed generation power at the moment of fault occurrence,the model can then be used to solve for the fluctuation range of injected current amplitudes in distributed power sources.The second stage is implemented after a malfunction occurs.In this stage,an affine optimization model is first established.This model is developed to characterizes the short-circuit current interval of a transmission line,and is constrained by the fluctuation range of the injected current amplitude of DG during normal operations.Finally,the range of the short-circuit current amplitudes of distribution network lines after a short-circuit fault occurs is predicted.The algorithm proposed in this article obtains an interval range containing accurate results through interval operation.Compared with traditional point value calculation methods,interval calculation methods can provide more reliable analysis and calculation results.The range of short-circuit current amplitude obtained by this algorithm is slightly larger than those obtained using the Monte Carlo algorithm and the Latin hypercube sampling algorithm.Therefore,the proposed algorithm has good suitability and does not require iterative calculations,resulting in a significant improvement in computational speed compared to the Monte Carlo algorithm and the Latin hypercube sampling algorithm.Furthermore,the proposed algorithm can provide more reliable analysis and calculation results,improving the safety and stability of power systems.

Practical Equivalent Method of Power Network Based on PSASP Short-Circuit Calculation Program

When using the ATP-EMTP software to simulate the electromagnetic transient of local line or nearby line for a large area power grid, it is necessary to divide the area grid into internal network and external network so as to make the external network simple for the restriction of the simulation calculation scale of software. Based on the existing equivalence method, in this paper we proposed a practical equivalent method using the short-circuit calculation program of PSASP. In this way, we completed the equivalence of northwest power grid operating at the large mode of summer 2016. By comparing the power flow and the characteristics of short-circuit in the system before and after the equivalence, this equivalent method is proved to be correct.

面向高比例新能源电网短路计算的机电-电磁融合电源模型

为解决新型电力系统计算与分析中传统电源模型难以继续适用的问题,提出一种面向高比例新能源电网短路计算的机电-电磁融合电源模型。一方面,考虑新能源电源通用控制策略作用下的外部输出特性,采用单一压控电流源等效电源侧海量新能源电源的整体输出;另一方面,根据短路电流故障的时间划分尺度,保留传统发电机工频基波意义下的外部输出特性,并采用3个大小不一且衰减速度不同的电压源等效电源侧所有发电机的整体输出。进而采用3个电压源和1个压控电流源并联的形式作为融合电源模型的结构,并采用内点法对融合电源模型参数进行在线辨识。通过短路仿真数据和电网实测故障录波数据验证了融合电源模型的适用性。

双制式列车接地系统的车-地联合牵引供电计算

为研究双制式列车在不同供电制式区段下的车体环流和车体—轴端电位分布,针对某型双制式列车构建列车接地系统的链式电路模型,并提出车-地一体化的交、直流区段联合牵引供电计算方法;同时,建立列车接地保护电阻优化模型,并分析2种列车接地系统配置方案.以国内某条双制式线路为算例进行仿真验证,研究结果表明:方案1相较方案2,列车的车体—轴端电位最大值降低36.58%~41.04%,车体环流最大值降低18.49%~22.97%;并且在方案1中,头尾车保护电阻设置为20 mΩ,可使列车的车体—轴端电位最大值为1.15 V,车体电流最大值为50.30 A,达到抑制车体—轴端电位的最优效果.该车-地一体化的交、直流区段联合牵引供电计算方法可适用于单一供电制式或者多供电制式列车接地系统的分析.

含主动限流控制的MMC-HVDC电网直流短路故障电流解析计算

传统的MMC-HVDC电网直流短路故障电流计算一般采用等效RLC电路方法。但是该方法没有考虑主动限流控制(ACLC)对等效电路RLC参数的影响,因此无法直接用于含ACLC的MMC-HVDC电网的直流短路故障电流的解析计算问题。该文通过将ACLC的控制效果等效为虚拟阻抗产生的压降,建立了体现ACLC的控制影响的等效RLC电路,以包含电流变化率限流环节的典型ACLC算法为例,推导了直流双极短路故障条件下ACLC的控制参数与虚拟阻抗的映射关系以及相应的直流故障电流解析表达式。最后,以基于含ACLC的半桥型四端MMC-HVDC电网为例,建立PSCAD/EMTDC仿真模型,对采用ACLC的直流双极短路故障电流解析计算结果进行了验证。结果表明:故障后10ms内直流短路故障电流的PSCAD仿真结果与所提出的含ACLC的等效RLC电路解析计算结果的平均误差在8.29%以内,最大误差在10.99%。因此,使用所提出的方法具有工程实用性。

Synthesis, Crystal Structure and Theoretical Calculation of [2-(2,6-Dioxacyclohexyl)-5-methoxylphenols]_2

The title compound [2-（2,6-dioxacyclohexyl）-5-methoxylphenols]2 was synthesized by the reaction of 2-hydroxyl-3-methoxylbenzaldehyde and 1,3-propanediol in the presence of DMF-DMS adduct and characterized by IR spectrum, UV-Vis spectrum and X-ray crystallography. The title compound belongs to monoclinic, space group P21 with a = 9.8967（10）, b = 8.2240（9）, c = 13.3654（14） A^°, β= 90.016（2）°, C11H14O4, Mr = 210.22, V = 1087.8（2）A^°3, Z = 4, Dc = 1.284 Mg/m^3, F（000） = 448,μ = 0.098 mm^-1, the final R = 0.0300 and wR = 0.0761 for 2070 observed reflections with I 〉 2σ（I）. The molecules are connected via intermolecular O-H…O hydrogen bonds into a 2D network structure. Moreover, theoretical investigations of the title compound with HF/6-31G（d） method were performed, and its stability, frontier molecular orbital composition and Mulliken charge distribution were also discussed. The compound is a bis-molecule. The two molecules stay together and could not be separated. Two crystallographically independent molecules exist in an asymmetric unit. The bond lengths and bond angles of the two molecules are slightly different with each other.

Structural and Vibrational Study of 2-（Quinolin-8-yloxy）-Acetic Acid based on FT-IR-Raman Spectroscopy and DFT Calculations

We have prepared the 2-（quinolin-8-yloxy）-acetic acid and characterized it by infrared and Raman spectroscopies in the solid phase. The Density Functional Theory （DFT） method, together with the 6-31G^＊ and 6-311＋＋ G^＊＊ basis sets, show that three stable molecules, for the anhydrous and monohydrated compounds were theoretically determined in the gas phase, and that probably the two more stable conformations are present in the solid phase of the monohydrated compound. The harmonic vibrational wavenumbers for the optimized geometries were calculated at B3LYP/6-31G^＊and B3LYP/6-311＋＋G^＊＊ levels. For a complete assignment of all the observed bands in the vibrational spectra the DFT calculations were combined with Pulay＇s scaled quantum mechanical force field （SQMFF） methodology in order to fit the theoretical Wavenumber values to the experimental ones. The characteristics of the electronic delocalization of all structures of both forms were performed by using natural bond orbital （NBO）, while the corresponding topological properties of electronic charge density are analysed by employing Bader＇s atoms in molecules theory （AIM）.

Synthesis,Structure and DFT Calculations of a Novel Copper(Ⅱ)Complex Based on Tert-butyl 2-[N-(tert-butyloxycarbonylmethyl)-2-picolyamino]acetate

A novel copper（Ⅱ） complex with tert-butyl 2-[N-（tert-butyloxycarbonylmethyl）-2-picolyamino]acetate（ampy） was synthesized and structurally characterized by elemental analysis, FT-IR spectrum, electrospray ionization-mass spectrometry, UV-vis spectrum and single-crystal X-ray diffraction, respectively. A mononuclear copper（II） complex with ampy, [Cu（ampy）Cl2]（1）, was formed irrespective of the metal-to-ligand molar ratios（[Cu2＋]：[ampy] = 0.5：1, 1：1, and 2：1） as a single product. Complex 1 crystallizes in the orthorhombic system, space group Pbca with a = 12.343（2）, b = 18.928（3）, c = 20.058（4） A, V = 4686.1（14） A3, Z = 8, Dc = 1.3349（4） g/cm3, F（000） = 1920, S = 1.016, R = 0.0693 and w R = 0.1721 for 3151 observed reflections with I 〉 2σ（I）. X-ray diffraction analysis reveals that the central copper（II） ion is bound by pyridyl N, tertiary amine N and carbonyl O atoms of the quadridentate ampy as well as two Cl anions, constructing a slightly distorted octahedral geometry. Complex 1 further constructs a stable 3D supramolecular architecture by intermolecular C–H…Cl hydrogen bonds. In addition, the molecular geometry was calculated by density functional theory（DFT/B3LYP） method with the basis sets（6-31＋G（d,p） for H, C, N, O and Cl atoms, and LANL2 DZ for Cu atom, respectively）. The calculated results show that the optimized geometrical parameters are in good agreement with the experimental data. Natural bond orbital（NBO） analysis and frontier molecular orbitals（FMOs） analysis were investigated at the same level.

Calculation of Radii and Density of ^(7-19)B Isotopes Using Effective Skyrme Force

Abstract In this study, the proton and neutron densities, charge densities, rms nuclear charge radii, rms nuclear mass radii, rms nuclear proton, neutron radii, and neutron skin thickness are calculated by using Harfree-Fock method with an effective nucleon-nucleon Skyrme interactions with SⅠ, SⅡ, SⅣ, T3, SKM, and SKM^＊ parameters. These nuclear properties for the neutron-rich isotopes of B （Boron） are presented. The calculated results are compared with the experimental and theoretical results of other researchers.

SELECTING CLUSTER MODEL IN Sn - BASED SOLDER ALLOY DESIGN WITH DV - X_α CALCULATION METHOD

Applying calculation method in alloy design should be an important tendency due to its characters of inexpensive cost, high efficiency and prediction. DOS calculations of AuSn, AsSn and SbSn Sn- based alloys have been investigated by employing DV - Xa method, in which different cluster models were adopted to calculate electron structure.It is proved that some regulations must be taken into ac- count in order to carry out alloy design calculation successfully,which are described in this paper in detail.

CALCULATION OF PHASE DIAGRAM FOR PSEUDOBINARY GaAs-InAs

ＣＡＬＣＵＬＡＴＩＯＮＯＦＰＨＡＳＥＤＩＡＧＲＡＭＦＯＲＰＳＥＵＤＯＢＩＮＡＲＹＧａＡｓ－ＩｎＡｓ￥Ｑｉａｏ，Ｈｕａｎ；Ｓｈｅｎ，Ｊｉａｎｙｕｎ；Ｌｉ，Ｇｕｏｘｕｎ（ＧｅｎｅｒａｌＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｆｏｒＮｏｎｆｅｒｒｏｕｓＭｅｔａｌ...

Synthesis, Structure and Quantum Mechanical Calculations of Methyl 2-(5-((Quinolin-8-yloxy)-methyl)-1,3,4-oxadiazol-2-ylthio)-acetate

The title compound was synthesized by the base catalyzed reaction of 5-（（quinolin- 8-yloxy）methyl）-1,3,4-oxadiazole-2（3H）-thione with methyl chloroacetate. The structure was supported by the spectroscopic data and unambiguously confirmed by single-crystal X-ray diffraction studies. It crystallizes from a methanol solution in the triclinic space group Pi with unit cell dimensions a = 7.4509（9）, b = 10.2389（12）, c = 12.2299（15）A, a = 74.771（2）, β = 77.956（2）, 7 = 69.263（2）°, V = 834.98（17） A3 and Z = 2. In order to gain some valuable insights into the molecular structure, the quantum mechanical calculations were performed using both HF and time-dependent density functional theory at the B3LYP/6-31G（d,p） level. The molecular geometry from X-ray determination of the title compound in the ground state has been compared using the Hartree-Fock （HF） and density functional theory （DFT） with the 6-31G（d） basis set. The calculated results show that the DFT and HF can well reproduce the structure of the title compound. The energetic behavior of the title compound was examined using the B3LYP method with the 6-31G（d） basis set. The harmonic vibrational frequencies calculated have been compared with the experimental FTIR and FT-Raman spectra. The restricted Hartree-Fock and density functional theory-based nuclear magnetic resonance （NMR） calculation procedure was also performed, and it was used for assigning the 13C and 1H NMR chemical shifts of the title compound. Moreover, molecular electrostatic potential and thermodynamic parameters of the title compound were investigated by theoretical calculations.

A 330-500 GHz Zero-Biased Broadband Tripler Based on Terahertz Monolithic Integrated Circuits

A 330-500 GHz zero-biased broadband monolithic integrated tripler is reported. The measured results show that the maximum efficiency and the maximum output power are 2% and 194μW at 348 GHz. The saturation characteristic test shows that the output i dB compression point is about -8.5 dBm at 334 GHz and the maximum efficiency is obtained at the point, which is slightly below the 1 dB compression point. Compared with the conventional hybrid integrated circuit, a major advantage of the monolithic integrated circuit is the significant improvement of reliability and consistency. In this work, a terahertz monolithic frequency multiplier at this band is designed and fabricated.

计及零直流电压控制的混合型MMC-HVDC输电系统短路电流计算方法

混合型模块化多电平换流器通过零压控制限制直流侧故障电流,为准确分析计算故障后直流侧电流,考虑不同阶段控制系统作用建立其等值模型,并设计故障电流求解计算方法。在零压控制前,构建混合型MMC的串联RLC等值模型,利用受控电流源表征交流馈入能量对故障电流的影响;在零压控制后,分析上下桥臂电压动态调整对故障电流的影响,构建相应混合型MMC等值模型提高计算精度。采用状态空间方程计算2个阶段的直流短路电流,通过电容电压置零来计及零压控制后子模块投切引起直流侧电压复零的特性。基于MATLAB/Simulink仿真平台搭建了双极混合型MMC直流输电系统,在不同工况下通过与详细电磁暂态仿真结果对比,验证了所提短路电流计算方法的准确性。

Process Analysis of Double Circuit Transmission Line Trip-out on the Same Tower Due to Lightning

In recent years,several failures of double circuit transmission line on the same tower due to lightning were happened in Beijing power grid.Although it can be reclosed successful,the lightning strike caused a grave threat to the power grid security.The cause of the accident and the accident process were studied for the sake of further understanding of the impact of lightning on power grid.As an example,110 kV double circuit transmission line(Xilong-line) was analyzed.At first,the system topology was given.Through the analysis on relay protection actions and the fault recorder data,over voltage on the insulator strings was calculated.Based on the analysis and the calculation,accident cause and the process were presented respectively.Secondly,it comes to the conclusion that the lightning failure was caused by counterattack.The wave of the lightning over voltage would spread to the not grounded neutral point of the transformers,and make the neutral protective gap breakdown,then cause freewheeling with the frequency of 50 Hz.As results of the relay protection,the double circuit transmission line all tripped out.Finally,the causes of the accident were proposed that included terrain features,large corner towers,strong thunderstorm weather and poor grounded contact of the tower.

Heat Production Rate Calculation of Thermal Manikin with a Compensation Method

Thermal manikin plays important roles in simulating thethermal state of human bodies to facilitate the evaluationof thermal comfort properties of various clothing assem-blies.Based on the heat balance principle and analysis ofrelated factors,the heat production rate of manikin isrecommended as an efficient evaluation index.Whereas,its inside heat production which occurs as a result of theexistence of temperature difference between its insidepart and outside surface,should not be ignored.Through a series of theoretical analysis and calculations,a compensative equation is deduced in this paper.

基于改进K-均值算法的双馈风场故障等值建模方法

新能源场站故障建模与分析方法是风电、光伏集中并网地区的故障特性和保护研究的基础,场站的多机等值建模多采用聚类算法,传统算法基于机组间运行状态的相似程度进行分群,然而将运行状态相近的机组分为一类并不等价于等值后模型精度更高。因此首先分析双馈风机故障电流的影响因素并作为分群指标,提出了一种改进的K-均值算法,算法以机组的故障电流等值误差作为样本间距,将同类机组的等值电气量作为算法的中心坐标;并结合了风电机组的等值误差与聚类算法的轮廓系数,整定了算法的分群数,构建了双馈风场的多机等值故障稳态模型,以提高网络短路计算的效率和精度。仿真结果验证了分群等值方法的合理性,所提方法相比单机等值模型与基于传统K-均值算法的多机等值模型精度均有一定的提升。

基于电路-概率理论的对神经电刺激不稳定性的定性分析

神经电刺激通过影响中枢或外周神经系统,进而治疗某疾病。在实际应用中,神经电刺激响应强度的不稳定性较为棘手,一般认为是由于电刺激扰动了神经轴突的膜电位,造成了神经电刺激强度的不稳定性。然而,由于神经电刺激缺乏宏观可计算模型,其膜电位扰动在电刺激中具体产生何种影响,长期以来未得到有效结论。该研究基于电路-概率理论,对神经电刺激不稳定性进行定性分析,有效地分析了膜电位的扰动对电刺激的影响。实验结果表明,在体实验数据和定性仿真的电流幅值-不稳定性曲线具有高度一致性,说明电路-概率理论可能对电刺激产生的膜电位扰动有较为准确的解释。该研究对神经电刺激的实际应用具有指导性意义。

计及风光相位特性和机-网间元件的短路电流计算方法

以风电和光伏为代表的新能源机组大量接入电网使得短路电流的精确计算校核变得困难。针对最新研究中构建新能源并网点短路电流映射模型时,未考虑相位映射及未计及箱变、集电线路、主变及送出线路等机端到并网点间电气元件(机-网间元件)影响的理论缺陷,文中提出一种计及新能源相位特性和机-网间元件影响的改进工程化电网短路电流计算方法。首先,结合国家运行规程对新能源机端电压-短路电流相位映射进行有补充意义的理论建模。其次,对该模型进行工程化处理,并提出一种将所得映射逆推至并网点的迭代计算方法。继而,将工程化后的并网点短路电流幅相映射模型应用到现有局部迭代计算中,得到改进的计算方法。最后,在PSCAD仿真软件中搭建新能源并网模型,验证了新能源机端电压-短路电流理论幅相映射的准确性、机端映射模型逆推至并网点的迭代计算理论有效性;在此基础上,在IEEE 39节点系统中对所得映射模型进行实验验证。结果表明所提改进计算方法能在一定程度上提升短路电流计算精度。

不同稳态工况对同步调相机单相短路后暂态温升的影响

同步调相机可在系统故障时为特高压直流换流站提供暂态无功支持,以保持系统电压稳定,防止直流换相失败,但调相机的高无功输出也会导致其温升增加,运行能力受限。系统故障前同步调相机的初始工况是影响其暂态温升和运行能力的一个重要因素,为研究系统低电压时调相机暂态温升与初始工况之间的关系,该文将调相机的电磁场、温度场与电网模型耦合,计算调相机不同稳态工况下的定转子损耗和温度分布;以线路单相短路故障引起的系统低电压为例,研究调相机初始工况对暂态温升的影响。在考虑各结构件最大容许温度的条件下,揭示调相机单相短路持续时间随初始工况的变化规律。研究结果为提升调相机在系统低电压下的暂态运行能力提供技术支持。