CN104318972A - Design method of X-ray flattening filter - Google Patents

Design method of X-ray flattening filter Download PDF

Info

Publication number
CN104318972A
CN104318972A CN201410610085.5A CN201410610085A CN104318972A CN 104318972 A CN104318972 A CN 104318972A CN 201410610085 A CN201410610085 A CN 201410610085A CN 104318972 A CN104318972 A CN 104318972A
Authority
CN
China
Prior art keywords
expansion device
ray
frustum
different levels
thickness
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.)
Granted
Application number
CN201410610085.5A
Other languages
Chinese (zh)
Other versions
CN104318972B (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.)
CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co Ltd
Original Assignee
CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co Ltd
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 CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co Ltd filed Critical CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co Ltd
Priority to CN201410610085.5A priority Critical patent/CN104318972B/en
Publication of CN104318972A publication Critical patent/CN104318972A/en
Application granted granted Critical
Publication of CN104318972B publication Critical patent/CN104318972B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21KTECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K1/00Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
    • G21K1/02Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21KTECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K1/00Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
    • G21K1/02Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators
    • G21K1/025Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diaphragms, collimators using multiple collimators, e.g. Bucky screens; other devices for eliminating undesired or dispersed radiation
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21KTECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K2201/00Arrangements for handling radiation or particles

Abstract

The invention relates to a design method of an X-ray flattening filter. A traditional design method of the flattening filter has the defects of being long in time consumed, high in cost and poor in flattening effect. The design of the multi-stage conical flattening filter is realized through the Monte Carlo method, a data fitting method and the like. The design method can realize a good flattening effect, and the whole process is clean and easy to realize. Meanwhile, the design cost is low and the design period is short.

Description

A kind of X ray expansion device method for designing
Technical field
The present invention relates to a kind of X ray expansion device method for designing, particularly for cone-beam x-ray, realize the method for designing of multi-stage taper expansion device based on the method such as Monte Carlo, data fitting.
Background technology
Traditional expansion device method for designing empirically designs the taper expansion device that is no more than three grades, and by revising expansion device structure repeatedly, test of many times reaches expansion effect after measuring.This method length consuming time, cost is high, and expansion effect is poor.
Also someone proposes a kind of expansion device method for designing.The method adopts the patten's design one generated step by step to have the expansion device of multistage frustum structure.Every one-level frustum all will through repeatedly adjusting, consuming time longer.In addition, due to the impact of ray scattering, often newly increase one-level frustum and all can frustum before be impacted.Therefore often newly-increased one-level frustum all needs to modify to frustum before, further increases the time of design, whole design need the number of times of amendment about to square being directly proportional of frustum quantity n.Along with the requirement of flatteness and the increase of frustum quantity, the global design time also will increase substantially.The design process of this method of formation is step by step also very complicated simultaneously.
Summary of the invention
In view of this, the object of the present invention is to provide a kind of X ray expansion device method for designing, the method realizes the design of multi-stage taper expansion device based on the method such as monte carlo method and data fitting.
For achieving the above object, the invention provides following technical scheme:
A kind of X ray expansion device method for designing, comprises the following steps:
1) X ray forms cone-beam x-ray after front collimation device, by place, expansion device bottom surface plane spatially angle be divided into n place, n is corresponding with expansion device frustum number, records each corresponding space angle θ 1, θ 2..., θ n;
2) be the base of m at expansion device base place plane thickness setting, base and radiogenic distance are L; According to parameters such as x-ray source beam energy, target structures, monte carlo method is utilized to calculate X ray on expansion device base plane by step 1) the roentgen dose X A of the n position that divides and expansion device center 0, A 1... .., A n, A 0for the roentgen dose X at center place;
3) monte carlo method is utilized to calculate the power spectrum of X ray on this position, n+1 place:
A i = Σ j = 1 k A ji , i = 0,1,2 , . . . . . n
4) design an one-level frustum, thickness is d 0, make the dosage of center ray after the decay of this frustum with the roentgen dose X A of edge nequal, A 0'=A n
5) wide beam modifying factor B is calculated
A 0 ′ = A n = B Σ j = 1 k A j 0 e - μ j d 0
Wherein, A j0centered by ray energy spectrum component; μ jfor the attenuation coefficient of correspondence; d 0for frustum thickness, solve above-mentioned equation and can obtain wide beam modifying factor B;
6) the attenuation length d of each frustum is calculated i
A i ′ = A n = B Σ j = 1 k A ji e - μ j d i
D ifor the unique unknown number in formula, solve the attenuation length d that above-mentioned equation can obtain frustum at different levels i;
7) according to the attenuation length d of frustum at different levels iand the space angle θ of correspondence icalculate the radius r of frustum at different levels iwith thickness l i;
8) according to the parameter r of frustum at different levels iand l iset up expansion device model, utilize monte carlo method to calculate the uniformity coefficient of ray after expansion device; Adjust the attenuation length d of frustum at different levels i, the flatteness that analog computation is obtained reaches requirement;
9) process expansion device according to the later parameter of adjustment, expansion device is arranged on the position set and measures; If actual
Flatteness reaches requirement and does not then modify; If actual flatteness does not reach requirement, on this basis expansion device is slightly made an amendment, make the actual flatteness of ray reach requirement.
Further, the material that expansion device is used is metal material.
Further, the material that expansion device is used is copper or stainless steel.
Further, step 2) in X-ray energy scope from 50keV to 200MeV.
Further, step 7) in the radius r of frustum at different levels iwith thickness l iand the attenuation length d of frustum at different levels iand the Space Angle θ of correspondence ithere is following recurrence relation in degree:
θ 0=0
l i=d i-1 cosθ i-1-d i cosθ i,i=1,2,....,n
r i = ( L - d 0 + Σ j = 1 i l j ) tan θ i , i = 1,2 , . . . . , n .
Further, step 1) to 8) in the foundation of roentgen dose X and power spectrum and expansion device model all at Geant4, or EGS, or other carry out based in the software of monte carlo method.
Beneficial effect of the present invention is: X ray expansion device method for designing provided by the invention can reach good expansion effect, and whole process flow is clear, is easy to realize; Meanwhile, design cost is lower, and the cycle is shorter.
Accompanying drawing explanation
In order to make object of the present invention, technical scheme and beneficial effect clearly, the invention provides following accompanying drawing and being described:
Fig. 1 is expansion device schematic diagram;
Fig. 2 is 6MeV medical accelerator expansion device structural section figure in embodiment;
Fig. 3 is the schematic flow sheet of the method for the invention.
Embodiment
Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described in detail.
Fig. 1 is expansion device schematic diagram, and Fig. 3 is the schematic flow sheet of the method for the invention.
Be designed to example with 6MeV accelerator for electron therapy expansion device be specifically described the method for the invention, Fig. 2 is 6MeV medical accelerator expansion device structural section figure in the present embodiment, the technical requirement of this figure is: 1, cone boss portion size tolerance is by the processing of IT8 level, and all the other dimensional tolerences press the processing of IT10 level; 2, material is brass.
Specifically comprise the following steps:
1) the distance L=128mm of expansion device base distance radiographic source target spot, ray subtended angle is 14 °.The expansion device that one has 10 grades of frustum structures is set, the space angle θ that each frustum is corresponding ibe respectively 1 °, 2 °, 3 °, 4 °, 5 °, 6 °, 8 °, 10 °, 12 °, 14 °.Expansion device Material selec-tion copper.
2) monte carlo method is utilized to calculate the dosage A of 6MeV ray in corresponding position 0, A 1... .., A 10, A 0for the roentgen dose X at center place.Calculate respective spectral distribution simultaneously.
A i = Σ j = 1 60 A ji , i = 0,1,2 , . . . . . 10
3) designing a thickness is d 0one-level frustum, make the dosage of the ray of center through frustum decay after equal with the roentgen dose X of edge.
A 0′=A 10 d 0=20mm
4) wide beam modifying factor B is calculated
A 0 ′ = A 10 = B Σ j = 1 60 A j 0 e - μ j d 0
5) attenuation length of frustum at different levels is calculated
A i ′ = A 10 = B Σ j = 1 60 A ji e - μ j d i , i = 1,2 , . . . , 10
6) according to the attenuation length d of frustum at different levels iand the space angle θ of correspondence icalculate the radius r of frustum at different levels iwith thickness l i.
7) according to the parameter r of frustum at different levels iand l iset up expansion device model, utilize monte carlo method to calculate the uniformity coefficient of ray after expansion device.Adjust the attenuation length d of frustum at different levels i, the flatteness that analog computation is obtained reaches requirement.
8) process expansion device according to the later parameter of adjustment, expansion device is arranged on the position set and measures.The basis of measurement result is once finely tuned expansion device structure, makes actual flatteness be less than 6%, as shown in Figure 3, reach designing requirement.
What finally illustrate is, above preferred embodiment is only in order to illustrate technical scheme of the present invention and unrestricted, although by above preferred embodiment to invention has been detailed description, but those skilled in the art are to be understood that, various change can be made to it in the form and details, and not depart from claims of the present invention limited range.

Claims (6)

1. an X ray expansion device method for designing, is characterized in that: comprise the following steps:
1) X ray forms cone-beam x-ray after front collimation device, by place, expansion device bottom surface plane spatially angle be divided into n place, n is corresponding with expansion device frustum number, records each corresponding space angle θ 1, θ 2..., θ n;
2) be the base of m at expansion device base place plane thickness setting, base and radiogenic distance are L; According to parameters such as x-ray source beam energy, target structures, monte carlo method is utilized to calculate X ray on expansion device base plane by step 1) the roentgen dose X A of the n position that divides and expansion device center 0, A 1... .., A n, A 0for the roentgen dose X at center place;
3) monte carlo method is utilized to calculate the power spectrum of X ray on this position, n+1 place:
A i = Σ j = 1 k A ji , i = 0,1,2 . . . . . n
4) design an one-level frustum, thickness is d 0, make the dosage of center ray after the decay of this frustum with the roentgen dose X A of edge nequal, A 0'=A n
5) wide beam modifying factor B is calculated
A 0 ′ = A n = B Σ j = 1 k A j 0 e - μ j d 0
Wherein, A j0centered by ray energy spectrum component; μ jfor the attenuation coefficient of correspondence; d 0for frustum thickness, solve above-mentioned equation and can obtain wide beam modifying factor B;
6) the attenuation length d of each frustum is calculated i
A i ′ = A n = B Σ j = 1 k A ji e - μ j d i
D ifor the unique unknown number in formula, solve the attenuation length d that above-mentioned equation can obtain frustum at different levels i;
7) according to the attenuation length d of frustum at different levels iand the space angle θ of correspondence icalculate the radius r of frustum at different levels iwith thickness l i;
8) according to the parameter r of frustum at different levels iand l iset up expansion device model, utilize monte carlo method to calculate the uniformity coefficient of ray after expansion device; Adjust the attenuation length d of frustum at different levels i, the flatteness that analog computation is obtained reaches requirement;
9) process expansion device according to the later parameter of adjustment, expansion device is arranged on the position set and measures; If actual flatteness reaches requirement, do not modify; If actual flatteness does not reach requirement, on this basis expansion device is slightly made an amendment, make the actual flatteness of ray reach requirement.
2. a kind of X ray expansion device method for designing according to claim 1, is characterized in that: described expansion device material used is metal material.
3. a kind of X ray expansion device method for designing according to claim 2, is characterized in that: described expansion device material used is copper or stainless steel.
4. a kind of X ray expansion device method for designing according to claim 1, is characterized in that: step 2) in X-ray energy scope from 50keV to 200MeV.
5. a kind of X ray expansion device method for designing according to claim 1, is characterized in that: step 7) in the radius r of frustum at different levels iwith thickness l iand the attenuation length d of frustum at different levels iand the Space Angle θ of correspondence ithere is following recurrence relation in degree:
θ 0=0
l i=d i-1cosθ i-1-d i cosθ i,i=1,2,....,n
r i = ( L - d 0 + Σ j = 1 i l j ) tan θ i , i = 1,2 , . . . . , n .
6. a kind of X ray expansion device method for designing according to claim 1, it is characterized in that: step 1) to 8) in the foundation of roentgen dose X and power spectrum and expansion device model all at Geant4, or EGS, or other carry out based in the software of monte carlo method.
CN201410610085.5A 2014-10-29 2014-10-29 A kind of X-ray expansion device method for designing Active CN104318972B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410610085.5A CN104318972B (en) 2014-10-29 2014-10-29 A kind of X-ray expansion device method for designing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410610085.5A CN104318972B (en) 2014-10-29 2014-10-29 A kind of X-ray expansion device method for designing

Publications (2)

Publication Number Publication Date
CN104318972A true CN104318972A (en) 2015-01-28
CN104318972B CN104318972B (en) 2017-07-18

Family

ID=52374191

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410610085.5A Active CN104318972B (en) 2014-10-29 2014-10-29 A kind of X-ray expansion device method for designing

Country Status (1)

Country Link
CN (1) CN104318972B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105608282A (en) * 2015-12-31 2016-05-25 重庆真测科技股份有限公司 Design method of ray filter

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6345114B1 (en) * 1995-06-14 2002-02-05 Wisconsin Alumni Research Foundation Method and apparatus for calibration of radiation therapy equipment and verification of radiation treatment
CN1346291A (en) * 1999-03-08 2002-04-24 加利福尼亚大学董事会 Correlated histogram representation of monte carlo derived medical accelerator photon-output phase space
JP2004041292A (en) * 2002-07-09 2004-02-12 Mitsubishi Electric Corp Dose calculation method
EP2157548A2 (en) * 2008-08-21 2010-02-24 Siemens Medical Solutions USA, Inc. System to estimate X-ray scatter
EP2248551A1 (en) * 2009-05-05 2010-11-10 7Sigma N.V. Method for the verification of a radiotherapy treatment apparatus
CN203787094U (en) * 2014-01-22 2014-08-20 上海联影医疗科技有限公司 Split type primary collimator

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6345114B1 (en) * 1995-06-14 2002-02-05 Wisconsin Alumni Research Foundation Method and apparatus for calibration of radiation therapy equipment and verification of radiation treatment
CN1346291A (en) * 1999-03-08 2002-04-24 加利福尼亚大学董事会 Correlated histogram representation of monte carlo derived medical accelerator photon-output phase space
JP2004041292A (en) * 2002-07-09 2004-02-12 Mitsubishi Electric Corp Dose calculation method
EP2157548A2 (en) * 2008-08-21 2010-02-24 Siemens Medical Solutions USA, Inc. System to estimate X-ray scatter
EP2248551A1 (en) * 2009-05-05 2010-11-10 7Sigma N.V. Method for the verification of a radiotherapy treatment apparatus
CN203787094U (en) * 2014-01-22 2014-08-20 上海联影医疗科技有限公司 Split type primary collimator

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
《物理学实用手册》编写组: "《物理学实用手册》", 31 January 1991, 福建人民出版社 *
K.VAN LAERE, W.MONDELAERS: "DESIGN OF FIELD FLATTENING FILTERS FOR A HIGH-POWER BREMSSTRAHLUNG CONVERTER BY FULL MONTE CARLO SIMULATION", 《RADIATION PHYSICS AND CHEMISTRY》 *
R.MOHAN, C.CHUI: "Energy and angular distributions of photons from medical linear accelerators", 《MEDICAL PHYSICS》 *
代志力 刘咸生: "一种基于MCNP设计X射线均整片的方法", 《全国射线数字成像与CT新技术研讨会论文集》 *
李泉凤 严慧勇: "用Monte Carlo方法设计X射线均整器", 《清华大学学报(自然科学版)》 *
陈宏芳 等: "《粒子探测技术》", 30 June 2009, 中国科学技术大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105608282A (en) * 2015-12-31 2016-05-25 重庆真测科技股份有限公司 Design method of ray filter
CN105608282B (en) * 2015-12-31 2018-09-11 重庆真测科技股份有限公司 A kind of ray filtration design method

Also Published As

Publication number Publication date
CN104318972B (en) 2017-07-18

Similar Documents

Publication Publication Date Title
Jabbari et al. Monte carlo simulation of siemens oncor linear accelerator with beamnrc and dosxyZnrc code
MX2022010978A (en) Shielding facility and method of making thereof.
Fulkerson et al. Dosimetric characterization and output verification for conical brachytherapy surface applicators. Part I. Electronic brachytherapy source
Carlsson et al. Point kernels and superposition methods for scatter dose calculations in brachytherapy
Van Der Heyden et al. Modelling of the focal spot intensity distribution and the off-focal spot radiation in kilovoltage x-ray tubes for imaging
CN104599735B (en) Gamma ray collimator for reference gamma radiation field
US10272268B2 (en) Method for estimating the dose administered by an external radiotherapy system
Goodall et al. Clinical implementation of photon beam flatness measurements to verify beam quality
CN104318972A (en) Design method of X-ray flattening filter
Harty et al. Absolute x‐ray dosimetry on a synchrotron medical beam line with a graphite calorimeter
Standaert et al. UCL irradiation facilities status
Kupfer et al. Commissioning of a PTW 34070 large‐area plane‐parallel ionization chamber for small field megavoltage photon dosimetry
Faddegon et al. Experimental depth dose curves of a 67.5 MeV proton beam for benchmarking and validation of Monte Carlo simulation
CN105699403B (en) Arrangement method for detecting vertical transillumination by directional X-ray
Tomic et al. Verification of cell irradiation dose deposition using a radiochromic film
EP3722747A3 (en) Volume measurement method, device, system and computer-readable storage medium
CN204314473U (en) The device of quick acquisition flushing-free checking film dosimetry response curve
US8742326B2 (en) Determination of beam parameters for asymmetrical photon beams
Vilsmeier et al. Investigation of the Effect of Beam Space-charge on Electron Trajectories in Ionization Profile Monitors
JP2014121635A (en) Particle beam medical treatment device and irradiation dose setting method
CN204067298U (en) A kind of radiation detector ionization chamber
Van Laere et al. Design of field flattening filters for a high-power bremsstrahlung converter by full Monte Carlo simulation
Halls Preliminary Investigation of Apparent Mass Loss in Objects Due to Image Blur using X-ray Radiography.
Lee et al. Measurements of the Beam Energy and Beam Profile of 100 MeV Proton Linac at KOMAC
US20160178795A1 (en) Method for measuring x-ray energy of an accelerator in an inspection system

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
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A Design Method for X-ray Equalizer

Effective date of registration: 20230605

Granted publication date: 20170718

Pledgee: Chongqing Bank Co.,Ltd. Major Branch

Pledgor: CHONGQING ZHENCE SCIENCE AND TECHNOLOGY Co.,Ltd.

Registration number: Y2023500000039

PE01 Entry into force of the registration of the contract for pledge of patent right