CN104916337B - A kind of data processing method of control rod calculus value measurement experiment - Google Patents

A kind of data processing method of control rod calculus value measurement experiment Download PDF

Info

Publication number
CN104916337B
CN104916337B CN201410095030.5A CN201410095030A CN104916337B CN 104916337 B CN104916337 B CN 104916337B CN 201410095030 A CN201410095030 A CN 201410095030A CN 104916337 B CN104916337 B CN 104916337B
Authority
CN
China
Prior art keywords
control rod
value
reactivity
measurement
rod
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410095030.5A
Other languages
Chinese (zh)
Other versions
CN104916337A (en
Inventor
杨晓强
李载鹏
姚进国
李文双
王建瑜
李友谊
王汗
董超
李宝库
叶刘锁
孙暖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Nuclear Power Corp
Original Assignee
Jiangsu Nuclear Power Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Nuclear Power Corp filed Critical Jiangsu Nuclear Power Corp
Priority to CN201410095030.5A priority Critical patent/CN104916337B/en
Publication of CN104916337A publication Critical patent/CN104916337A/en
Application granted granted Critical
Publication of CN104916337B publication Critical patent/CN104916337B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • G21C17/104—Measuring reactivity
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Plasma & Fusion (AREA)
  • General Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The invention belongs to a kind of reactor start-up experimental technique field, and in particular to a kind of data processing method of control rod calculus value measurement experiment.The method includes the steps of：Step 1: to the continuous low discharge injected clear water of reactor core or boric acid solution；Step 2: mobile control rod, record data point；Step 3: to data point linear fit；Step 4: obtain changes of reactivity；Step 5: obtain the differential value and integral worth of reactivity；Step 6: repeat step two to five；Step 7: design rod position is mapped with designing reactive differential value；Step 8: the figure of step 5 is fitted；Step 9: measurement of comparison value and design load.This method, which is reduced, even ignores influence of the three-dimensional effect to result of the test, reduces reactivity measurement value error, and precision improves about 5%；By carrying out six fittings of unitary to design load, the error that conventional method is brought to measured value linear fit is reduced, while reduce standard deviation.

Description

A kind of data processing method of control rod calculus value measurement experiment
Technical field
The invention belongs to a kind of reactor physics firing test technical field, and in particular to a kind of control rod calculus value Measure the data processing method of experiment.
Background technology
According to the requirement of nuclear safety codes, the reactor core after reloading must carry out neutronics firing test, to verify reactor core The accuracy and security of design, the value measurement of control rod calculus is important one therein.VVER-1000/428 types react Heap control rod is divided into 10 groups, wherein the 8th, 9,10 group is regulating rod group, the 10th group is main regulating rod group.Experiment is by adjusting boron method Measured for the 9th, 10 group of control rod, can obtain the rod position changes delta Hi of each control rod moving period, and correspondingly Reactor core changes of reactivity Δ ρ i, by Δ ρ i/ Δs Hi obtain rod differential be worth, integral worth is ∑ Δ ρ i.Interpolation obtain with Reactivity measurement value under the identical rod position of design load, evaluate whether to meet according to formula (measured value/design load -1) * 100% Acceptance criteria.This traditional data processing method comes with some shortcomings, and is mainly manifested in following two aspect：
Reactivity measurement value bigger error.Experiment process and data processing are to be based on Point reactor nutron kinetics model, and in Sub- flux is three-dimensional spatial distribution in reactor core, and the change of control rod position not only changes the total netron-flux density of reactor core, Change local distribution, but currently used reactivity meter can not reflect partial power's change of reactor core.Using control rod underthrust All data points of all data points recorded afterwards or short time after, if carrying out described point interception reactivity, it can not reduce The influence of three-dimensional effect, cause reactivity measurement value bigger error.
Measured value causes bigger error with design load control methods.In order to be contrasted with design load, generally by measured value according to Rod position corresponding to design load carries out linear fit, then evaluates whether that satisfaction is tested according to (measured value/design load -1) * 100% Receive criterion.Because measured value inherently has larger measured deviation, caused using the method to measured value linear fit total Bigger error.
The content of the invention
It is an object of the invention to provide the data processing method that a kind of control rod calculus is worth measurement experiment, reduce data The error that processing procedure is brought, improve test accuracy.
In order to solve the above problems, a kind of data processing method of control rod calculus of the invention value measurement experiment, The method includes the steps of：
Step 1: to the continuous low discharge injected clear water of reactor core or boric acid solution；
Step 2: when homophony rod group control rod is from about 85% core height underthrust, inserted by the reactivity meter minute book cycle The change of reactor core reactivity before and after rod, record reactivity；
Step 3: the point that reactivity meter records in the t times after selection control rod movement, is intended using least square method Close, obtain changes of reactivity fit line, a upper cycle data is fitted with same method, obtains cycle changes of reactivity Fit line；
Step 4: using the reactive fit line obtained by step 3, obtain reactivity meter before control rod movement record it is anti- Answering property variation delta ρ 'iThe changes of reactivity amount Δ ρ " that reactivity meter records after being moved with control rodi, each cycle rod value Δ ρiComputational methods be Δ ρi=(Δ ρ 'i+Δρ″i)/2；
Step 5: pass through Δ ρi/ Δ Hi is obtained in HMeasurementLocate rod differential value, pass through Σ Δs ρiIntegral worth is obtained, wherein The depth that Δ Hi moves for this cycle control rod, HMeasurementFor the rod bit depth in this cycle, with HMeasurementFor abscissa, the control rod of measurement Differential value records the point for ordinate；
Step 6: repeat step two arrives step 5, the discrete point of multiple control rod differential value measured values is obtained；
Step 7: to design rod position HDesignFor abscissa, control rod differential Value design value is that ordinate does control rod differential Value design value changes figure；
Step 8: being fitted to the control rod differential Value design value of step 7, corresponding relational expression is obtained；
Step 9: according to control rod differential Value design value relational expression, by measuring stick position HMeasurementBring formula into, measured Rod position HMeasurementUnder rod differential Value design value, obtain rod differential Comparison of Value with measurement, then according to (measured value/design load- 1) * 100% evaluate whether to meet acceptance criteria.
In the step 2, in the movement of control rod, reactive need is controlled between ± 0.05 β eff.
In described step three, t >=200s.
In described step eight, six fittings of unitary are carried out to control rod differential Value design value.
The advantageous effects of this method are：By the analysis and optimization to test data point selection, can reduce even Ignore influence of the three-dimensional effect to result of the test, reduce the reactivity measurement value error of conventional data analysis method, precision improves About 5%.By carrying out six fittings of unitary to design load, the error that conventional method is brought to measured value linear fit is reduced, Standard deviation is reduced to 0.012pcm/cm by 0.054pcm/cm.
Brief description of the drawings
Fig. 1 is control rod reactivity worth data processing schematic diagram；
Fig. 2 is the data point that reactivity meter records in experiment；
Fig. 3 is reactive differential Value design value and control rod away from reactor core bottom position graph of a relation；
Fig. 4 is that reactive differential is worth measured value with control rod away from reactor core bottom position graph of a relation.
Embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.
This method comprises the following steps：
Step 1: to the continuous low discharge injected clear water of reactor core or boric acid solution；
Step 2: when homophony rod group control rod is from about 85% core height underthrust, inserted by the reactivity meter minute book cycle The change of reactor core reactivity before and after rod, reactivity control is between ± 0.05 β eff, record reactivity, as shown in Figure 2；
Step 3: the data point that t >=200s reactivity meters record after selection control rod movement, is carried out using least square method Fitting, obtains changes of reactivity solid line, a upper cycle data is fitted with same method, obtains cycle changes of reactivity Solid line；
Step 4: as shown in figure 1, using the reactive fit line obtained by step 3, reactivity before control rod movement is obtained The changes of reactivity amount Δ ρ ' of instrument recordiThe changes of reactivity amount Δ ρ " that reactivity meter records after being moved with control rodi, Mei Yizhou Phase rod value Δ ρiComputational methods be Δ ρi=(Δ ρ 'i+Δρ″i)/2；
Step 5: pass through Δ ρi/ Δ Hi is obtained in HMeasurementLocate rod differential value, pass through ∑ Δ ρiIntegral worth is obtained, wherein The depth that Δ Hi moves for this cycle control rod, HMeasurementFor the rod bit depth in this cycle, with HMeasurementFor abscissa, the control rod of measurement Differential value records the point for ordinate,
Step 6: repeat step two arrives step 5, the discrete point of multiple control rod differential value measured values, such as Fig. 4 are obtained It is shown；
Step 7: as shown in figure 3, to design rod position HDesignFor abscissa, control rod differential Value design value is done for ordinate Control rod differential Value design value changes figure；
Step 8: carrying out six fittings of unitary to Fig. 3 design load, corresponding relational expression is obtained；
Step 9: the formula obtained according to Fig. 3 curve matchings, by measuring stick position HMeasurementBring formula into, obtain measuring stick position HMeasurementUnder rod differential Value design value.Rod differential Comparison of Value is obtained with measurement, then according to (measured value/design load -1) * 100% evaluates whether to meet acceptance criteria.

Claims (3)

1. a kind of data processing method of control rod calculus value measurement experiment, the method includes the steps of：
Step 1: to the continuous low discharge injected clear water of reactor core or boric acid solution；
Step 2: when homophony rod group control rod is from about 85% core height underthrust, before reactivity meter minute book cycle plunger The change of reactor core reactivity afterwards, record reactivity；
Step 3: the point that reactivity meter records in the t times after selection control rod movement, is fitted using least square method, obtained To changes of reactivity fit line, a upper cycle data is fitted with same method, obtains cycle changes of reactivity plan Zygonema；
Step 4: using the reactive fit line obtained by step 3, the reactivity that reactivity meter records before control rod moves is obtained Variation delta ρi' changes of reactivity amount Δ the ρ recorded with reactivity meter after control rod movementi", each cycle rod value Δ ρi's Computational methods are Δ ρi=(Δ ρi′+Δρi″)/2；
Step 5: pass through Δ ρi/ Δ Hi is obtained in HMeasurementLocate rod differential value, pass through ∑ Δ ρiObtain integral worth, wherein Δ Hi For the depth of this cycle control rod movement, HMeasurementFor the rod bit depth in this cycle, with HMeasurementFor abscissa, the control rod differential of measurement It is worth and records the point for ordinate；
Step 6: repeat step two arrives step 5, the discrete point of multiple control rod differential value measured values is obtained；
Step 7: to design rod position HDesignFor abscissa, control rod differential Value design value is that ordinate does control rod differential value Design load variation diagram；
Step 8: being fitted to the control rod differential Value design value of step 7, corresponding relational expression is obtained；
Step 9: according to control rod differential Value design value relational expression, by measuring stick position HMeasurementBring formula into, obtain measuring stick position HMeasurementUnder rod differential Value design value, rod differential Comparison of Value is obtained with measurement, then according to (measured value/design load -1) * 100% evaluates whether to meet acceptance criteria.
2. a kind of data processing method of control rod calculus value measurement experiment according to claim 1, its feature exist In：In the step 2, in the movement of control rod, reactive need is controlled between ± 0.05 β eff.
3. a kind of data processing method of control rod calculus value measurement experiment according to claim 1, its feature exist In：In described step three, t >=200s.
CN201410095030.5A 2014-03-14 2014-03-14 A kind of data processing method of control rod calculus value measurement experiment Active CN104916337B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410095030.5A CN104916337B (en) 2014-03-14 2014-03-14 A kind of data processing method of control rod calculus value measurement experiment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410095030.5A CN104916337B (en) 2014-03-14 2014-03-14 A kind of data processing method of control rod calculus value measurement experiment

Publications (2)

Publication Number Publication Date
CN104916337A CN104916337A (en) 2015-09-16
CN104916337B true CN104916337B (en) 2017-11-24

Family

ID=54085340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410095030.5A Active CN104916337B (en) 2014-03-14 2014-03-14 A kind of data processing method of control rod calculus value measurement experiment

Country Status (1)

Country Link
CN (1) CN104916337B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105895174B (en) * 2016-05-06 2017-07-04 中国核动力研究设计院 A kind of presurized water reactor control rod worth computational methods
CN106373624B (en) * 2016-10-20 2018-06-22 中核核电运行管理有限公司 A kind of dynamic carves stick controlling test stick and is inserted into starting point automatic judging method
CN109166639A (en) * 2018-08-28 2019-01-08 中国核动力研究设计院 Reactor control rod measurement data preparation method and differential value capture method
CN109979625B (en) * 2019-03-28 2022-02-18 江苏核电有限公司 Test method for measuring symmetry of reactor core proliferation characteristics
CN110010256B (en) * 2019-04-29 2021-07-20 岭澳核电有限公司 Method and device for measuring control rod value of nuclear power station in subcritical state

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1147968A (en) * 1966-05-16 1969-04-10 Gen Electric Improvements in method of controlling nuclear reactor
US4251321A (en) * 1967-12-15 1981-02-17 General Electric Company Nuclear reactor utilizing plutonium
US4642216A (en) * 1984-03-12 1987-02-10 Westinghouse Electric Corp. Control rod cluster arrangement
US4877575A (en) * 1988-01-19 1989-10-31 Westinghouse Electric Corp. Core reactivity validation computer and method
CN102714063A (en) * 2009-11-06 2012-10-03 希尔莱特有限责任公司 Systems and methods for controlling reactivity in a nuclear fission reactor
CN103268728A (en) * 2013-04-27 2013-08-28 国家电网公司 Method for constructing power system dynamic simulation pressurized water reactor control system model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59184885A (en) * 1983-04-06 1984-10-20 株式会社日立製作所 Fuel assembly

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1147968A (en) * 1966-05-16 1969-04-10 Gen Electric Improvements in method of controlling nuclear reactor
US4251321A (en) * 1967-12-15 1981-02-17 General Electric Company Nuclear reactor utilizing plutonium
US4642216A (en) * 1984-03-12 1987-02-10 Westinghouse Electric Corp. Control rod cluster arrangement
US4877575A (en) * 1988-01-19 1989-10-31 Westinghouse Electric Corp. Core reactivity validation computer and method
CN102714063A (en) * 2009-11-06 2012-10-03 希尔莱特有限责任公司 Systems and methods for controlling reactivity in a nuclear fission reactor
CN103268728A (en) * 2013-04-27 2013-08-28 国家电网公司 Method for constructing power system dynamic simulation pressurized water reactor control system model

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SHB-5临界装置铀水栅堆芯硼微分价值和总后备反应性的测量;黄文楼 等;《核动力工程》;19971031;第18卷(第5期);392-418页 *
压水堆控制棒价值误差分析;付学峰 等;《核科学与工程》;20130630;第33卷(第2期);第152-158页 *
反应堆物理试验用便携式数字反应性仪实堆考验试验;黄礼渊;《核电子学与探测技术》;20100430;第30卷(第4期);472-475页 *
插棒法空间效应修正因子计算中MC方法的应用;赵柱民 等;《核动力工程》;20051231;第26卷(第6期);531-534页 *
秦山核电厂反应堆控制棒价值和硼价值测量试验;蒋树庆;《核动力工程》;19930228;第14卷(第1期);第73-78页 *

Also Published As

Publication number Publication date
CN104916337A (en) 2015-09-16

Similar Documents

Publication Publication Date Title
CN104916337B (en) A kind of data processing method of control rod calculus value measurement experiment
CN101105986B (en) Reactor reactivity measuring method
CN107910083B (en) A kind of nuclear power plant's measurement of rod worth method
CN107122545B (en) Method for accurately calculating space-time neutron distribution in nuclear reactor
CN105788666A (en) Nuclear reactor critical process control method
CN105895174A (en) Method for calculating value of control rod for pressurized water reactor
CN105425174A (en) Method and apparatus for obtaining leakage magnetic field of core type transformer
CN109273119B (en) Method for optimizing the position of a neutron detector during the measurement of high reactivity on a critical device
Kaviani et al. Integral form of the control rod calibration curve in the new core configuration of HWZPR using rod insertion method
CN109192341B (en) Large reactivity measurement method based on three-dimensional space-time dynamics
CN105808910A (en) Irradiation deformation calculation method of zirconium-niobium pressure pipe in CANDU reactor
CN104898155A (en) Method for measuring undercritical degree of undercritical system based on average neutron energy
CN102800373B (en) A kind of method of fission foil detector scale reactor absolute power
Salam et al. Measurement of control rod reactivity and shut down margin of 3 MW TRIGA Mark-II research reactor using analogue and digital I&C system
CN108828651B (en) Active neutron analysis method for uranium plutonium content in waste cladding
CN112668165B (en) Reactor core design system and application
CN204731395U (en) A kind of multistation magnetization degree measurement mechanism
Liang et al. Application of dynamic rod worth measurement in a VVER
CN109324070A (en) A kind of passive neutron assay method of uranium plutonium content in cladding waste
Dupont et al. Development of critical experiments to benchmark moderator temperature reactivity worth
Borland Evaluation of end-of-cycle moderator temperature coefficient measurements for Davis-Besse Unit 1
CN104049225A (en) Automatic magnetic powder magnetization degree measuring device
Vasiliev et al. Validation studies of computational scheme for high-fidelity fluence estimations of the Swiss BWRs
Dudley et al. The reactor core neutronic model for the pebble bed modular reactor
Borodkin et al. Application of ex-vessel neutron dosimetry combined with in-core measurements for correction of neutron source used for RPV fluence calculations

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant