首页 | 官方网站   微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Values of effective dose (E) were calculated for the entire range of incident directions of broad parallel photon beams for selected photon energies using the Monte Carlo N-Particle (MCNP) transport code with a hermaphroditic phantom. The calculated results are presented in terms of conversion coefficients transforming air kerma to effective dose. This study also compared the numerical values of E and H(E) over the entire range of incident beam directions. E was always less than H(E) considering all beam directions and photon energies, but the differences were not significant except when a photon beam approaches some specific directions (overhead and underfoot). This result suggests that the current H(E) values can be directly interpreted as E or, at least, as a conservative value of E without knowing the details of irradiation geometries. Finally, based on the distributions of H(E) and E over the beam directions, this study proposes ideal angular response factors for personal dosimeters that can be used to improve the angular response properties of personal dosimeters for off-normal incident photons.  相似文献   

2.
The scientific background of the standard procedure for determination of the mean electron energy at the phantom surface (E0) from the half-value depth (R50) has been studied. The influence of energy, angular spread and range straggling on the shape of the depth dose distribution and the R50 and Rp ranges is described using the simple Gaussian range straggling model. The relation between the R50 and Rp ranges is derived in terms of the variance of the range straggling distribution. By describing the mean energy imparted by the electrons both as a surface integral over the incident energy fluence and as a volume integral over the associated absorbed dose distribution, the relation between E0 and different range concepts, such as R50 and the maximum dose and the surface dose related mean energy deposition ranges, Rm and R0, is analysed. In particular the influence of multiple electron scatter and phantom generated bremsstrahlung on R50 is derived. A simple analytical expression is derived for the ratio of the incident electron energy to the half-value depth. Also, an analytical expression is derived for the maximum energy deposition in monoenergetic plane-parallel electron beams in water for energies between 2 and 50 MeV. Simple linear relations describing the relative absorbed dose and mass ionization at the depth of the practical range deposited by the bremsstrahlung photons generated in the phantom are derived as a function of the incident electron energy. With these relations and a measurement of the extrapolated photon background at Rp, the treatment head generated bremsstrahlung distribution can be determined. The identification of this photon contamination allows an accurate calculation of the absorbed dose in electron beams with a high bremsstrahlung contamination by accounting for the difference in stopping power ratios between a clean electron beam and the photon contamination. The absorbed dose determined using ionization chambers in heavily photon contaminated (10%) electron beams may be too low--by as much as 1.5%--without correction.  相似文献   

3.
In photon beam convolution, the distribution of energy deposition about a primary photon interaction site due to charged particles set in motion at that site is represented by the primary kernel. Energy deposited due to scattered photons, bremsstrahlung, and annihilation photons is represented by the scatter kernel. As the energy deposited in each kernel voxel is normalized to the energy imparted at the interaction site, it is known as a fractional energy distribution. In terma-based convolution, where kernels are normalized to total energy imparted at the interaction site and are convolved with the terma in the dose calculation process, the sum of fractional energies contained in the primary kernel is equal to the ratio of collision kerma (Kc) to terma (T) corresponding to the energy spectrum used to generate the kernel. Since the ratio of collision kerma to terma increases with depth as the beam hardens, the integral fractional energy in a primary kernel formed for the spectrum at the surface is less than the ratio Kc/T at depth. This causes primary dose to be increasingly underestimated with depth and scatter dose to be increasingly overestimated. Single polyenergetic convolution (using polyenergetic primary and scatter kernels formed using a polyenergetic primary photon spectrum) is thus not as rigorous as if a separate convolution is performed for each energy component. The ratio of true primary dose to single polyenergetic primary dose increases almost linearly with depth and is almost equal to the Kc/T ratio. Primary and scatter dose are calculated correctly if a single polyenergetic convolution is performed in terms of Kc (for primary) and T-Kc (for scatter), where the kernels are weighted sums of monoenergetic kernels normalized to Kc and T-Kc. With this method, it is ensured that total primary energy deposited due to primary photon interactions in a unit mass at a point is equal to Kc at that point.  相似文献   

4.
To investigate validity of scatter correction by the TEW method in 201Tl imaging, we performed an experimental study using the gamma camera with the capability to perform the TEW method and a plate source with a defect. Images were acquired with the triple energy window which is recommended by the gamma camera manufacturer. The result of the energy spectrum showed that backscattered photons were included within the lower sub-energy window and main energy window, and the spectral shapes in the upper half region of the photopeak (70 keV) were not changed greatly by the source shape and the thickness of scattering materials. The scatter fraction calculated using energy spectra and, visual observation and the contrast values measured at the defect using planar images also showed that substantial primary photons were included in the upper sub-energy window. In TEW method (for scatter correction), two sub-energy windows are expected to be defined on the part of energy region in which total counts mainly consist of scattered photons. Therefore, it is necessary to investigate the use of the upper sub-energy window on scatter correction by the TEW method in 201Tl imaging.  相似文献   

5.
PURPOSE: To determine the dose perturbation effects at the tissue-metal implant interfaces in head and neck cancer patients treated with 6 MV and 10 MV photon beams. METHODS AND MATERIALS: Phantom measurements were performed to investigate the magnitude of dose perturbation to the tissue adjacent to the titanium alloy implants with (100 mu and 500 mu thick) and without hydroxylapatite (HA) coating. Radiographic and radiochromic films were placed at the upper (and lower) surface of circular metal discs (diameter x thickness: 15 x 3.2, 48 x 3.2, 48 x 3.8 mm2) in a solid water phantom and were exposed perpendicular to radiation beams. The dosimeters were scanned with automatic film scanners. Using a thin-window parallel-plate ion chamber, dose perturbation were measured for a 48 x 3.2 mm2 disc. RESULTS: At the upper surface of the tissue-dental implant interface, the radiographic data indicate that for 15 x 3.2 mm2 uncoated, as well as 100 mu coated discs, dose perturbation is about +22.5% and +20.0% using 6 MV and 10 MV photon beams, respectively. For 48 x 3.2 mm2 discs, these values basically remain the same. However, for 48 x 3.8 mm2 discs, these values increase slightly to about +23.0% and +20.5% for 6 MV and 10 MV beams, respectively. For 48 x 3.2 mm2 discs with 500 mu coating, dose enhancement is slightly lower than that obtained for uncoated and 100 mu coated discs for each beam energy studied. At the lower interface for 15 x 3.2 mm2 and 48 x 3.2 mm2 uncoated and 100 mu coated discs, dose reduction is similar and is about -13.5% and -9.5% for 6 MV and 10 MV beams, respectively. For 48 x 3.8 mm2 discs, dose reduction is about -14.5% and -10.0% for 6 MV and 10 MV beams, respectively. For 48 x 3.2 mm2 discs with 500 mu coating, the dose reduction were slightly higher than those for uncoated and 100 mu coated discs. CONCLUSIONS: For the beam energies studied, dose enhancement is slightly larger for the lower energy beam. The results of dose perturbation were similar for 100 mu coated and uncoated discs. These results were slightly lower for the 500 mu coated discs but are not clinically significant. The dosimetry results obtained from radiochromic films were similar to the ones obtained from radiographic film. The dose enhancement results obtained from ion chamber dosimetry are higher than those obtained from film dosimetry. The ion chamber data represent the data at "true" tissue-titanium interface, whereas the ones obtained from film dosimetry represent the data at film-titanium interface.  相似文献   

6.
The phantom scatter correction factor Sp of megavoltage photon beams can be accurately described using a three-Gaussian fit. The model leads to six parameters, with which Sp(r) is described as a smooth function of the field radius r for beam qualities in the range from 60Co up to 25 MV. The parameters allow Sp values to be calculated at intermediate beam energies and for any field shape. Calculated Sp(X, Y) values for rectangular fields (X, Y) can be subsequently used as reference values to compare with measured Sp(X, Y) values, for example when appraising a new beam.  相似文献   

7.
Successful radiotherapy requires accurate dosimetry for treatment verification. Existing dosimeters such as ion chambers, TLD, and diodes have drawbacks such as relatively long measurement time and poor spatial resolution. These disadvantages become serious problems for dynamic-wedged beams. Thus the clinical use of dynamic wedges requires an improved dosimetry method. X-ray film may serve this purpose. However, x-ray film is not clinically accepted as a dosimeter for photon beams, because it overresponds to photons with energies below about 400 keV. This paper presents and develops a method which was initially proposed by Burch to improve the dose response of x-ray film in a phantom. The method is based on placing high-atomic number foils next to the film. The foils are used as filters to preferentially remove low-energy photons. The optimal film and filter configuration in a phantom was determined using a mathematical scheme derived in this study and a Monte Carlo technique (ITS code). The optimal configuration thus determined is as follows: the filter-to-film distance of 6 mm and the filter thickness of 0.15 mm for percent depth-dose measurement; the distance of 1 cm and the thickness of 0.25 mm for off-axis (dose) ratio measurement. The configuration was then tested with photon beams from a 4 MV linac. The test result indicates that the in-phantom dose distribution based on the optimal configuration agrees well with those measured by ion chambers.  相似文献   

8.
A new method for producing very narrow and intense 50 MV bremsstrahlung beams with a half-width as low as 35 mm at a distance of 1 m from the target is presented. Such a beam is well suited for intensity modulation using scanned photon beams. An algorithm has been developed to minimize the width of the bremsstrahlung beam generated in a multilayer target by varying the individual layer thicknesses and atomic numbers under given constraints on the total target thickness and the mean energy of the transmitted electrons. Under such constraints the narrowest possible bremsstrahlung beam is obtained with a target composed of layers of monotonically increasing atomic number starting with the lowest possible value at the entrance side where the electrons impinge. It is also shown that the narrowest photon beam profile is associated with the highest possible forward photon yield. To be able to use the optimized target clinically it is desirable to be able to collect and stop all the electrons that are transmitted through the target. The electrons are most efficiently collected if they are kept close together, i.e. by minimizing the multiple scatter of the electrons and consequently the half-width of the generated bremsstrahlung beam. This is achieved by a thin low-atomic-number target. A dedicated electron stopper has been developed and integrated with the purging magnet. When the electron stopper is combined with a purging magnet, a primary photon collimator and a multileaf collimator, almost all of the transmitted electrons and their associated bremsstrahlung contamination can effectively be collected. The narrow photon beams from thin low-atomic-number targets have the additional advantage of producing the hardest and most penetrative photon spectrum possible, which is ideal for treating large deep-seated tumours.  相似文献   

9.
This paper presents two methods for absorbed dose determination with ionization chambers at short distance from 60Co and 192Ir brachytherapy sources. The methods are modifications of the Bragg-Gray and large cavity principles given in the IAEA code of practice for high- and medium-energy photon beams. A non-uniformity correction factor to account for the non-uniform electron fluence in the air cavity is introduced into the methods. The absorbed dose rates were determined from ionization chamber measurements at distances between 1.5 and 5.0 cm from the brachytherapy sources. The agreement between the two methods is excellent in 60Co brachytherapy dosimetry. For 192Ir dosimetry, the difference is less than 2.5% at all distances. In absorbed dose rate calculations with the 60Co brachytherapy source, the ratios between calculated and experimentally determined absorbed dose rates are 0.987 and 0.994 depending on the method used for absorption and scatter correction. In 192Ir dosimetry, the large cavity principle gives almost identical values to those which can be obtained with the AAPM recommendations. Using the chambers according to the Bragg-Gray principle in 192Ir dosimetry, the agreement with AAPM calculated absorbed dose rates is within 2.5% at all distances. The uncertainty, expressed as one standard deviation, in the experimentally determined absorbed dose is estimated to be between 3 and 4%.  相似文献   

10.
In routine dosimetry we assume separability of the collimator (Sc) and phantom (Sp) scatter components that together comprise the total scatter factor (Sc,p). In practice, the addition of blocking also affects the photon fluence attributable to the treatment head and flattening filter in a complicated way. The reduced aperture blocks out some of the head scatter contribution, while the block and tray add back secondary scatter. In the following we present techniques for directly measuring the aperture effect on Sc in air or in a full-scatter phantom. The change in Sc is found to be a scaleable quantity that can be modelled as a simple linear fit to the ratio of projected open-to-blocked equivalent square fields. Measurements have been made for 6, 18 and 24 MV photon beams on one Varian 2500 and two Varian 2100c accelerators. Results indicate a progressive loss of collimator scatter contribution with increased field blocking that is amplified with increasing energy. Block and tray scatter only contribute significantly to Sc for large fields and treatment distances of 80 cm or less. Application of these corrections in monitor unit calculations is presented.  相似文献   

11.
PURPOSE: To evaluate dose variations at bone/titanium interfaces in an experimental model designed to simulate postoperative radiotherapy in patients with mandibular reconstructions using a titanium hollow-screw osseointegrating reconstruction plate (THORP) system. MATERIALS AND METHODS: The model consisted of a 25 x 25 x 10 mm3 block of fresh bovine femoral diaphysis, to the surface of which a segment of THORP system reconstruction plate was fixed by means of a solid titanium screw 4 mm in diameter and 10 mm in length. Using specially designed thermoluminescent dosimeters (TLD) 2 mm in diameter and 0.13 mm in thickness, dose measurements were carried out at four distances from the screw axis (0.1, 0.3, 0.6, and 1 mm). 60Co and 6-MV photon beams were used at incidences both perpendicular and parallel ("axial") to the screw axis. RESULTS: For 6-MV X-ray beams incident perpendicular to the screw axis, the maximum dose enhancement (due to backscatter) and the maximum dose reduction (due to attenuation) at the bone/titanium interface were 5% (+/- 2%) and 6% (+/- 2%), respectively. The corresponding values for 60Co beams were 6% (+/- 5%) and 10% (+/- 5%). For the axial incidences, a maximum dose enhancement of 5-7% was noted for both 6-MV X-rays and 60Co for beams incident on the surface containing the THORP plate segment, whereas beams incident on the opposite surface induced only a very small dose enhancement (2-3%). CONCLUSION: Using a new experimental model, TLD measurements showed only marginally significant dose variations at bone/titanium interfaces around THORP screws, all measured values being very close to the uncertainty limits (+/- 5%) associated with the method. For both 60Co and 6-MV beams, dose variations appeared smaller for axial than for perpendicular incidences. Because photon beams used in head and neck cancer treatment are most often directed parallel to the screw axes, these results suggest that failures of prosthetic osseointegration are unlikely to be explained by an overdosage at the bone/titanium interface.  相似文献   

12.
PURPOSE: To measure the effect of silicon diode detectors used for in vivo dosimetry on beam characteristics and determine whether this effect is clinically significant. METHODS AND MATERIALS: Commercially available photon and electron diodes were placed on the central axis of photon and electron beams. The beam characteristics were measured for 6- and 10-MV photon and 6-20-MeV electron energies from a Varian Clinac 1800 medical linear accelerator. Water was used for the medium, and measurements were made for various clinically common field sizes and depths. RESULTS: Beam attenuations along the central axis were 10 and 7.5% for 6- and 10-MV photons, respectively. Electron beam dose reductions were between 13 and 25% for 20-6-MeV electrons. Photon beam flatness varied up to 7% at different depths, but the symmetry was not affected much. Electron beam flatness and symmetry were significantly changed to as much as 18 and 6%, respectively. CONCLUSION: Use of diode detectors on central axis of photon and electron beams for in vivo dosimetry causes significant attenuation and alteration of the beam characteristics. The percentage of the volume affected is significant (e.g., 23% of the volume in a 4 x 4 field gets 10% less dose for a 6-MV photon beam), especially if these diodes are used for in vivo dosimetry on the central axis every day for every treatment, as is done in some clinics. Other beam parameters such as penumbra and skin dose are also affected. It is therefore recommended that the diodes be used only as needed.  相似文献   

13.
The accuracy of quantitation of radionuclide distributions in human tissue with the scintillation camera is decreased by attenuation and scatter of photons. If scatter correction is applied satisfactorily, narrow beam attenuation can be applied. In this article, a scatter correction technique, the channel ratio (CR) method, is introduced. The CR scatter correction method is proposed for quantitation of the radionuclide distribution in organs. The improvement in the geometrical resolution was measured and examples of clinical images are presented. In this method, the change in the ratio of counts from two symmetrical adjacent energy windows straddling the energy photopeak was used to eliminate the contribution of scattered photons during imaging with 99mTc. The theory and methods for the empirical affirmation are described. To apply the CR scatter correction method, two constants, the ratio of primary photons G and the ratio of scattered photons H in the same windows, were determined. Different sized sources in varying depths of water were imaged. When the source activities were quantified after scatter correction with the CR method, the measurements ranged from 96%-108% in comparison to the reference value in 100 mm water. The scatter fraction increased from 0.20 in 10 mm water to 1.44 in 200 mm water. The geometrical resolution expressed as full width at tenth maximum in 150 mm water improved by 30.4% and was restored to the value of the geometrical resolution in air. The CR scatter correction method is a simple method to correct for scatter in order to facilitate accurate quantitation of the radionuclide distribution during imaging with a scintillation camera.  相似文献   

14.
Accuracy of dose delivery at low monitor unit setting is studied for a dual photon energy linear accelerator. Dose delivered per MU is found to be constant for both the photon beams for MU settings above 30. For lower MUs there is definite deviation from the calibrated value and the error is found to be increasing as fewer MUs are set for dose delivery. This dose/MU ratio at low MU setting is found to be dose-rate dependent, showing an increasing trend with pulse repetition frequency (PRF). Also, the dosimetric ratio is observed to be mode dependent; its value for an 18 MV beam is almost double that observed in the case of a 6 MV beam at very low MU setting. The magnitude of this error should be determined for each energy so that appropriate corrections can be applied if very low MUs are to be used.  相似文献   

15.
Backscatter factors for orthovoltage x-ray beams are fitted with the expression BSF = 1 + ms/(s + n), where s is field size and m and n are adjustable parameters. The accuracy is better than 0.5%. Good results were also obtained for Cs-137 gamma rays at depth 0.12 cm and Co-60 gamma rays at depth 0.5 cm using the same expression, which failed, however, for peak scatter factors at high x-ray energies.  相似文献   

16.
A realistic photon beam model based on Monte Carlo simulation of clinical linear accelerators was implemented in a convolution/superposition dose calculation algorithm. A primary and an extra-focal sources were used in this beam model to represent the direct photons from the target and the scattered photons from other head structures, respectively. The effect of the finite size of the extra-focal source was modeled by a convolution of the source fluence distribution with the collimator aperture function. Relative photon output in air (Sc) and in phantom (Scp) were computed using the convolution method with this new photon beam model. Our results showed that in a 10 MV photon beam, the Sc, Sp (phantom scatter factor), and Scp factors increased by 11%, 10%, and 22%, respectively, as the field size changed from 3 x 3 cm2 to 40 x 40 cm2. The variation of the Sc factor was contributed mostly by an increase of the extra-focal radiation with field size. The radiation backscattered into the monitor chamber inside the accelerator head affected the Sc by about 2% in the same field range. The output factors in elongated fields, asymmetric fields, and blocked fields were also investigated in this study. Our results showed that if the effect of the backscattered radiation was taken into account, output factors in these treatment fields can be predicted accurately by our convolution algorithm using the dual source photon beam model.  相似文献   

17.
A strategy for Compton scatter correction in brain SPECT images was proposed recently. It assumes that two radioisotopes are used and that a significant portion of photons of one radioisotope (for example, Tc99m) spills over into the low energy acquisition window of the other radioisotope (for example, Tl201). We are extending this approach to cases of several radioisotopes with mutual, multiple and significant photon spillover. In the example above, one may correct not only the Tl201 image but also the Tc99m image corrupted by the Compton scatter originating from the small component of high energy Tl201 photons. The proposed extension is applicable to other anatomical domains (cardiac imaging).  相似文献   

18.
Not until the last decade has flexible intensity modulated three-dimensional dose delivery techniques with photon beams become a clinical reality, first in the form of heavy metal transmission blocks and other beam compensators, then in dynamic and segmented multileaf collimation, and most recently by scanning high-energy narrow electron and photon beams. The merits of various treatment unit and bremsstrahlung target designs for high-energy photon therapy are investigated theoretically for two clinically relevant target sites, a cervix and a larynx cancer both in late stages. With an optimized bremsstrahlung target it is possible to generate photon beams with a half-width of about 3 cm at a source to axis distance (SAD) of 100 cm and an initial electron energy of 50 MeV. By making a more compact treatment head and shortening the SAD, it is possible to reduce the half-width even further to about 2 cm at a SAD of 70 cm and still have sufficient clearance between the collimator head and the patient. One advantage of a reduced SAD is that the divergence of the beam for a given field size on the patient is increased, and thus the exit dose is lowered by as much as 1%/cm of the patient cross section. A second advantage of a reduced SAD is that the electron beam on the patient surface will be only about 8 mm wide and very suitable for precision spot beam scanning. It may also be possible to reduce the beamwidth further by increasing the electron energy up to about 60 MeV to get a photon beam of around 15 mm half-width and an electron beam as narrow as 5 mm. The compact machine will be more efficient and easy to work with, due to the small gantry and the reduced isocentric height. For a given target volume and optimally selected static multileaf collimator, it is no surprise that the narrowest possible scanned elementary bremsstrahlung beam generates the best possible treatment outcome. In fact, by delivering a few static field segments with individually optimized scan patterns, it is possible to combine the advantage of being able to fine tune the fluence distribution by the scanning system with the steeper dose gradients that can be delivered by a few static multileaf collimator segments. It is demonstrated that in most cases a few collimator segments are sufficient and often a single segment per beam portal may suffice when narrow scanned photon beams are employed, and they can be delivered sequentially with a negligible time delay. A further advantage is the increase of therapeutically useful photons and improved patient protection, since the pencil beam is only scanned where the leaf collimator is open. Consequently, some of the problems associated with dynamic multileaf collimation such as the tongue and groove and edge leakage effects are significantly reduced. Fast scanning beam techniques combined with good treatment verification systems allow interesting future possibilities to counteract patient and internal organ motions in real time.  相似文献   

19.
Recently the compensator has been shown to be an in expensive and reliable dose delivery device for photon beam intensity-modulated radiation therapy (IMRT). The goal of IMRT compensator design is to produce an optimized primary fluence profile at the patient's surface obtained from the optimization procedure. In this paper some of the problems associated with IMRT compensator design, specifically the beam perturbations caused by the compensator, are discussed. A simple formula is derived to calculate the optimal compensator thickness profile from an optimized primary fluence profile. The change of characteristics of a 6 MV beam caused by the introduction of cerrobend compensators in the beam is investigated using OMEGA Monte Carlo codes. It is found that the compensator significantly changes the energy spectrum and the mean energy of the primary photons at the patient's surface. However, beam hardening does not have as significant an effect on the percent depth dose as it does on the energy spectrum. We conclude that in most situations the beam hardening effect can be ignored during compensator design and dose calculation. The influence of the compensator on the contaminant electron buildup dose is found to be small and independent of the compensator thickness of interest. Therefore, it can be ignored in the compensator design and included as a correction into the final dose distribution. The scattered photons from the compensator are found to have no effect on the surface dose. These photons produce a uniform low fluence distribution at the patient's surface, which is independent of compensator shape. This is also true for very large fields and extremely asymmetric and nonuniform compensator thickness profiles. Compared to the primary photons, the scattered photons have much larger angular spread and similar energy spectrum at the patient's surface. These characteristics allow the compensator thickness profile and the dose distribution to be calculated from the optimized fluence profile of primary photons, without considering the scattered photons.  相似文献   

20.
We investigated the impact of air cavities in head and neck cancer patients treated by photon beams based on clinical set-ups. The phantom for investigation was constructed with a cubic air cavity of 4 x 4 x 4 cm3 located at the centre of a 30 x 30 x 16 cm3 solid water slab. The cavity cube was used to resemble an extreme case for the nasal cavity. Apart from measuring the dose profiles and central axis percentage depth dose distribution, the dose values in 0.25 x 0.25 x 0.25 cm3 voxels at regions around the air cavity were obtained by Monte Carlo simulations. A mean dose value was taken over the voxels of interest at each depth for evaluation. Single-field results were added to study parallel opposed field effects. For 10 x 10 cm2 parallel opposed fields at 4, 6 and 8 MV, the mean dose at regions near the lateral interfaces of the cavity cube were decreased by 1 to 2% due to the lack of lateral scatter, while the mean dose near the proximal and distal interfaces was increased by 2 to 4% due to the greater transmission through air. Secondary build-up effects at points immediately beyond the air cavity cube are negligible using field sizes greater than 4 x 4 cm2. For most head and neck treatment, the field sizes are usually 6 x 6 cm2 or greater, and most cavity volumes are smaller than our chosen dimensions. Therefore, the influence of closed air cavities on photon interface doses is not significant in clinical treatment set-ups.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司    京ICP备09084417号-23

京公网安备 11010802026262号