Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 RESEARCH ARTICLE www.ijera.com OPEN ACCESS Case Study of the Chaq-Chaq Dam Failure: Parameter Estimation and Evaluation of Dam Breach Prediction Models Dr. KawaZedanAbdulrahman Abstract On 4th of February, 2006 at about 10:00 pm.Chaq-Chaq dam failed due to overtopping. The fall of 131.2 mm of rain over a 24-hour period was recorded at Sulaimani metrological gage station, which is located about 7.5Kmsouth-east of the dam. As a result, the reservoir level rose, the dam has been overtopped and finally breached near the spillway at the right abutment. Fortunately no human lives loss nor important structure destruction were reporteddue to the dam failure. The aim of this paper is to estimate the flood hydrograph passing through Chaq-Chaq dam breach using measured breach geometry as input to unsteady option of HEC RAS 4.1.0 and calibrating the breach formation time to obtain the measured maximum water surface at ChaqChaq Bridge (1.36 km downstream of dam axis). In addition the recent breach prediction models were evaluated to check their accuracy in predicting the breach geometry, breach formation time and peak breach discharge. I. Introduction Chaq-Chaq dam is located about 2 km NE of Sulaimani city (Iraq). Fig. 1 shows a satellite image of the area between Chaq-Chaq dam and Chaq-Chaq Bridge.Chaq-Chaq dam is a zoned earth dam of central clay core and gravelly shell as shown in Fig.2. Chaq-Chaq dam was designed and built by engineers of little experience in the field of dam design and construction. As a recognized design problem,one of the major mistakes was the building of the spillway beside the dam in the same valley not as a separate structure. The spillway wall has been made vertical. Compaction of an embankment near a vertical wall is notrecommended in constructing embankmentdams because this procedure will produce a weak bond at the interface of the wall and the embankment(FEMA, 2005).In addition; the required compaction for the materials close to the vertical wall will not be gained. This weak-compacted portion will be weaker compare to the other well-compacted portions of the dam. Therefore, the dam breached close to the spillway rather than other locations. In order to check the accuracy of existing breach prediction models in predicting the breach geometry, breach formation time and peak breach discharge; a bathymetric survey after the dam failure has been carried out to obtain the breach geometry. Extensive interviews with the surrounding habitants, owners of the tourism cabinets, and directorate of security have been done to gather information about the breach formation time and the highest water level at Chaq-Chaq Bridge.The breach formation time and www.ijera.com the highest water level at Chaq-Chaq Bridge is used as input to calibrate the HEC RAS 4.1 (Brunner, 2010 a,b) to achieve the maximum flood discharge passing through the dam breach as it will be presented as the followings. Figure 1: Satellite image showing the area between Chaq-Chaq dam and Chaq-ChaqBridge. 109|P a g e Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 Spillway www.ijera.com Dam Body Core Flow Direction Shell Figure 2: Photo of Chaq- Chaq dam after failure. II. Breach GeometryData Bathymetric survey has been carried out to obtain the breach geometries;Table 1 shows the geometry parameters of Chaq-Chaq dam and its breach. Table 1: Geometry parameters of Chaq-Chaq dam. Parameter Height Top Upstream Downstream Breach Breach Breach Dam of dam width slope (v:h) slope (v:h) Bottom Average Top crest H (m) (m) width width width level (m) (m) (m) (masl) Value 14.5 9 1:3 1:2 29.6 38 46 780 III. Breach hydraulic data Due to insufficient spillway capacity ChaqChaq dam was overtopped and then failed. According to a local witness (who was the formal responsible of the dam and his house was located about 100 m far from the dam) the maximum depth of water above the dam crest was between 0.5 − 0.6 m. So, he was also estimated the breach formation time to be between 1 to 1.5 hours. In addition the maximum water level due to the dam failure flood at Chaq-Chaq bridge which is located about 1.36 km downstream of Para meter Value the dam has been decided based on eyewitness accounts. The maximum water level at the bridge was around 759.4-759.5maslas corresponded to 20-30 cm below the lower cord of the bridge. There was a security team at the bridge to prevent peoples from passing the bridge because there was a potential of bridge failure due to the flood before and during the dam failure. The flood extent at the bridge was seen by the security team. Table 2 shows some of the hydraulic parameters of Chaq-Chaq dam and ChaqChaq Bridge. Table 2: Hydraulic parameters of the Chaq-Chaq dam and bridge. Depth Breach Reservoir Reservoir Reservoir Spillway Spillwa Minimum of formati storage at storage at storage at length y crest stream bed overt on NPL El. El. 780.0 El. 780.6 (m) level level at the oppin time 777.5 MCM MCM (masl) bridge g (m) (hr) MCM (masl) 0.6 1-1.5 www.ijera.com 1.4 2.344 2.55 15 777.5 754.4 High cord level of the bridge (masl) low cord level of the bridge (masl) 761 759.7 110|P a g e Maxim um water level at the bridge (masl) 759.6 Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 IV. Upstream and downstream cross sections data One-dimensional dam breach hydraulic modelof HEC-RAS is used frequentlyto predict the flood inundation area due to a dam breachflood through the downstream valley. It was found thatHEC-RAS performed well, with relatively good agreement between predicted and measured water levels(Yochum etal.,2008) and(Gee, 2010). HEC-RAS modeling system is a public domain model developed by the US Army Corp of Engineers (Brunner, 2010 a,b). It performs onedimensional (1D) steady and unsteady flow simulations on a full network of natural or man-made open channels. Additionally, it has the ability to model storage areas and dam break problems as well as bridges and culverts hydraulics. In order to model the flooding in the stream valley using HEC‐RAS; cross sections data are required. In this study a topographic map of 1m interval is obtained in AutoCAD format from the municipality of Sulaimanya.Then, the river reach in the Chaq-Chaq system extending over a length of 4.15 km from upstream end of the reservoir to the downstream portion of the damis considered for analysis. The cross sections data of the river reach aredeveloped by AutoCAD Civil 3D 2013, by using this software the main channel as well as right and left overbank have beennoted and coded in the hydraulic model. A number of 21 cross sections at the upstream of the dam are used to model the reservoir area and19 cross sections were developed at the downstream portion. Extra cross sections were added by interpolation at a maximum distance of 75 m. The values of Manning’s roughness coefficient were entered directly into the cross section editor to describe the channel and overbanks. These values were determined by visual inspection and satellite imagesbased on guidance fromChow (1959). The Manning’s roughness coefficient values were set at 0.028 for main channel and the two overbanks. These values have been assumed because the stream reach under study is clean with stones and high flow stages are expected during the dam break analysis (Parhi etal, 2012). www.ijera.com www.ijera.com V. HEC-RAS Model The unsteady option of HEC-RASrequires the breach geometry and breach formation time as input in order to model a dam breach flood.The breach geometry is readily available from the bathymetric survey but breach formation time is still a matter of uncertainty (1-1.5 hrs). Breach formation time is the most sensitive parameter in developing a hydraulic model for dam break problems and breach hydrograph development. Therefore, in this study it is attempted to calibrate the breach formation time through simulation of breach flood using HEC-RAS 4.1 unsteady model. For calibration of Breach formation time value; the observedWSE at the downstream bridge has been considered. A weir coefficient of 1.1 was used in this analysis; the trigger time of breach is set such that it corresponds to the time of peak of a developed inflow hydrograph as it will be explained in the next paragraph. At that time the water surface elevation was equal to 780.57 m which is close to the observed water surface elevation ( 780.5 − 780.6 m). This equality in the simulated and the observed WSE proves that the developed inflow hydrograph is accurate and that there was under-estimate for the inflow hydrograph in the design of Chaq-Chaq Dam. VI. Boundary conditions The upstream boundary condition is modeled using the flood hydrograph corresponding to the measured 131.2 mm rainfall depth during 24 hrs on a 151 𝑘𝑚2 of catchment area. The flood hydrograph is developed from contributing catchments using NRCSunit hydrograph (UH) method. The NRCS dimensionless UH is a synthetic unit hydrograph in which the discharge is expressed by the ratio of discharge to peak discharge and the time by the ratio of time to the time of rise of the unit hydrograph(Chow etal., 1988). Fig.3 shows the developed inflow flood hydrograph. Details on how to develop flood hydrograph using NRCS UH can be found in McCuin (2005). The downstream boundary condition is set to normal depth and an approximate water surface slope is assumed for the friction slope. 111|P a g e Flood (m3/sec) Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 www.ijera.com 200.00 180.00 160.00 140.00 120.00 100.00 80.00 60.00 40.00 20.00 0.00 0 10 20 30 40 50 60 70 80 90 Time (hrs) Figure 3: Inflow flood hydrograph. 1.60hrs.Fig. 4 provides a plot of modeled water surface profiles at different times of the simulation and Fig. 5 shows the outflow flood hydrograph through the dam breach. Each hydrograph was routed through the downstream reach to produce different water surface elevations at the downstream bridge; the results of the model at the bridge location are shown in Table 3. The percentage of errorsbetween the predicted water surface elevation and theobservedwater surface elevation at the downstream bridge are depicted in the Table 4. VII. Initial conditions The WSE upstream of the dam is set to 780 m which is the crest elevation of the dam; while WSE at the downstream reach is set such that 2 m depth of water is existing. VIII. Hydraulic model result Using of surveyed dam breach geometry combined with standard approaches for developing the upstream hydrograph boundary conditiona HEC RAS model was developed to generate different breach hydrograph corresponding to different breach formation times, namely 1.25 hrs, 1.50hrsand Chaq Chaq Ch1 Legend WS 04FEB2006 2340 WS 04FEB2006 2200 780 Elevation (m) Ground 770 0 First S ection at upstream 750 Dam Axis Bridge section 760 1000 2000 3000 4000 Main Channel Distance (m) Figure 4: water surface profiles at different times of simulation corresponding to 1.6 hrsof breach formation time. www.ijera.com 112|P a g e Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 www.ijera.com Plan: 1.1 weir 1.5 River: Chaq Chaq Reach: Ch1 RS: 19.5 1200 Legend Flow (m3/s) 1000 Flow Beginning of dam failure @time of 2200 800 600 Peak outflow@ time of 2340 400 200 0 1200 1800 04Feb2006 2400 0600 05Feb2006 Time Figure 5: Outflow flood hydrograph due to Chaq-Chaq dam breach. Table 3: Results of the model at the bridge location at different breach formation time. Breach River Time Q Min W.S. Crit E.G. E.G. VelC formati Sta of peak Total Ch Elev W.S. Elev Slope hnl. on time (hrs) El (m) (m) (m) (m) (m/s) (𝒎𝟑 /s (hr) (m) ) 1.25 1.5 1.6 Just upstrea m of bridge Just upstrea m of bridge Just upstrea m of bridge Flow Area (𝒎𝟐 ) Top Width (m) Froud e# Chl 4FEB200 62325 929.5 754.4 760.60 758.4 760.8 0.0005 17 2.37 492.4 287.2 0.32 04FEB20 06 2335 919.1 754.4 759.92 758.4 760.2 0.0009 71 2.94 392.1 208.1 0.43 04FEB20 06 2340 915 754.4 759.55 758.4 760 0.0012 81 3.24 352.9 165.3 0.48 Table 4: Departures of estimated and observed water surface elevations at Chaq-Chaq Bridge corresponding to different BFT. Breach formation time Estimated WSE at ObservedWSE at ChaqDifference between (hr) Chaq-Chaq bridge Chaq bridge (m) estimated and measured using HEC RAS (m) WSE (m) 1.25 760.60 759.50 1.10 1.5 759.92 759.50 0.42 1.6 759.55 759.50 0.05 A comparison of the predictedWSE with the observed WSE at the bridge indicates that a breach formation time of 1.60 hrs may be considered the most accurate value, with a differenceof 0.05 m in WSE.The modeling indicates a peak breach flood discharge of979.2 m3 s and this value attenuates at the bridge to 915.4 m3 s. IX. Existing Empirical Breach Prediction Models www.ijera.com Simulation of dam breach floods is essential to characterize and identify hazards due to hypothetical dam failures. Hydraulic modelssuch as HEC-RASare often used for the analysis of downstream impacts resulting from potential damfailures. Estimation of the dam breach parameters, such as formation time, width and side slopes,has usually done external to the hydraulic model. If input breach parameters cannot be predicted with sufficient accuracy, more conservative parameters and associated increased costs may be 113|P a g e Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 required (Wahl, 1997).This paper aims to check the reliability of the existing breach prediction methodologies in estimating the breach parameters of Chaq-Chaq dam. Four important breach parameters namelytop width, average width, breach formation time and peak discharge pass through the beach are estimated by thefollowing breach prediction models Froehlich (1995, 2008), Xu and Zhang (2009) and Pierce etal. (2010)and the results are compared to the measured values (breach geometries) and HEC-RAS www.ijera.com output values (breach formation time and peak discharge). Froehlich (1995) model was selected based on the results obtained by(Wahl, 2004) which showed that this model is more accurate than other existing prediction models up to the time the paper was published.Froehlich (1995) as cited in (Wahl, 2004), developed the following formulas based on 75, 34 and 31 case studies for Bavg , Tf and Q p ; respectively: Bavg =0.1803× k o × Vw 0.32 × Hb0.19 … … … … … … … … … … . . … … … … … … … (1) Tf = 0.00254 × Vw0.53 × Hb−0.9 … … … … … … … … … … . . … … … … … … … … … (2) 1.24 Q p = 0.607Vw 0.295 Hw … … … … … … … … … … … . … … … … … … … … … … . . (3) Where K o = constant = 1.4 if there is overtopping and 1.0 if else, Z=1.4 if there is overtopping, otherwise Z=1.0, Vw = volume of reservoir at the time of failure, hb =height of breach, Bavg = average width, Tf = breach formation time and Q p = peak discharge. Froehlich (2008) developed the following formulas based on 74, 23 case studies for Bavg , and t f ; respectively: Bavg = 0.27K o Vw 0.32 Hb 0.04 … … … … … . . … … … … … … . … … … . . … … … … . . (4) Tf = 0.0175 Vw gHb 2 … … … … … … … … … … … … … … … . … … … … … … … … . (5) Where K o = constant = 1.3 if there is overtopping and 1.0 if else, Z=1.0 if there is overtopping, if not Z=0.7. Xu and Zhang (2009) proved that his model is more accurate than other models. This model was developed using 182 case studies to estimateBt ,Bavg , Tf and Q p ; respectively: Bt Hd = 1.062 Hb Hr 0.092 1 0.508 Vw 3 Hw eB 1 … … … … … … … … … … … … … … … … . (6) With Bt = top width of the breach, Hd = dam height, Hr = 15m , Hw = height of water at the time of failure,B1 = b3 + b4 + b5 , in which b3 = 0.061, 0.088, and −0.089 for dams with core-walls, concrete faced dams, and homogeneous or zoned-fill dams, b4 = 0.299 and −0.239 for overtopping and seepage erosion or piping, b5 = 0.411, −0.062, and−0.289 for high, medium, and low dam erodibility Bave Hd = 0.787 Hb Hr 0.133 1 0.652 Vw 3 Hw eB 2 … … … … … . … … … … … … … … … . (7) with B2 = b3 + b4 + b5 , in which b3 = −0.041, 0.026, and − 0.226 for dams with core-walls, concrete faced dams, and homogeneous or zoned-fill dams, respectively, 𝑏4 = 0.149 𝑎𝑛𝑑 − 0.389 for overtopping and seepage erosion/piping, respectively, 𝑏5 = 0.291, −0.14, and − 0.391 for high, medium, and low dam erodibility, respectively Tf Hd = 0.304 Tr Hr 0.707 1 1.228 Vw 3 Hw eB 3 … … … … … … … … … … … … . … … … . (8) withTr = 1 hr., B3 = b3 + b4 + b5 , in which b3 = −0.327, −0.674, and − 0.189 for dams with core-walls, concrete faced dams, and homogeneous/ zoned-fill dams, respectively, b4 = −0.579 and − 0.611 for overtopping and seepage erosion/piping, respectively, b5 =−1.205, −0.564, and 0.579 for high, medium, and low dam erodibility, respectively. Qp gVw 5/3 Hd = 0.175 Hr 0.199 1 Vw 3 Hw −1.274 eB 4 … … … … . … … … … … . … … … (9) withB4 = b3 + b4 + b5 , in which b3 = −0.503, −0.591, and − 0.649 for dams with core-walls, concrete faced dams, and homogeneous or zoned-fill dams, respectively, b4 = −0.705 and − 1.039 for overtopping and www.ijera.com 114|P a g e Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 www.ijera.com seepage erosion/piping, respectively, b5 = −0.007, −0.375, and − 1.362 for high, medium, and low dam erodibility, respectively. Pierce (2010) showed that his developed multiple-regression model using 87 case studies is more accurate than the Froehlich (1995) in predicting peak-discharge through an embankment dam breach. Q p = 0.038Vw 0.475 Hd 1.09 … … … … … … … … … … … … … … … . . . … … … … … (10) X. Comparison with considered empirical breach prediction models Applying the above equations to Chaq-Chaq dam failure yields the results shown in Table 5. Table 5: Results of empirical models applied to Chaq-Chaq dam failure. Breach Observed HEC RAS Froehlich Froehlich Xu and Pierce parameter value Prediction (1995) (2008) Zhang(2009) (2010) Prediction Prediction 𝐏𝐫𝐞𝐝𝐢𝐜𝐭𝐢𝐨𝐧𝐛 𝑩𝒕 (m) 𝑩𝒂𝒗𝒆 (m) 𝑻𝒇 (hr) 46 38 1-1.5 N.A. 𝒎𝟑 N.A. N.A. 1.6 979.2 𝑸𝒑 ( ) 𝑺 1.13 N.A. Side slope Z a. Obtained by using values of Z and 𝐵𝑎𝑣𝑒 b. Medium dam erodibility is assumed. 57.25𝑎 47.1 0.57 1364 51𝑎 43.8 0.62 N.A. 54.5 38.4 1.17 1274 N.A. N.A. N.A. 809 1.4 1 N.A. N.A. Table 6: Percentage of errors between predicted and measured values. Breach Froehlich Froehlich (2008) Xu and Zhang Pierce parameter (1995) Prediction (2009) (2010) Prediction Prediction Prediction 𝑩𝒕 (m) 𝑩𝒂𝒗𝒆 (m) 𝑻𝒇 (hr)* 𝑸𝒑 ( 𝒎𝟑 𝑺 )* 24 23.9 -64.4 39.2 10.8 15.2 -61.2 N.A. 18.5 1.0 -26.8 30.1 N.A. N.A. N.A. -17.3 *HEC RAS results are considered as measured values Generally, all the models over-predict the breach top width andthe averagewidth.This trend of the models to over-predict the breach size may be attributed to the fact that they are developed based on the assumption of breach forms in a shape of trapezoid, while Chaq-Chaq breach has a vertical side near the spillway which may be considered as an odd case. However,Xu and Zhang (2009) predicts the average breach width more accurate than others, where the percentage of the error between the predicted and the measured values is 1% as shown in table (6). While the predicted breach top width using Froehlich (2008) appears to be better than others with an error of 10.8% and Xu and Zhang comes in the second order with an error of 18.5%. All the used models under-predict the breach formation time, includingXu and Zhang (2009) who was the best where it gives an error of -26.8%. Froehlich (1995 and 2008) errors are -65.4% and 62.4%; respectively. www.ijera.com The predictedpeak flood discharge using the considered empirical models shows that most of these equations tend to over-predictthe value of this parameter; except Pierce (2010) which yields a value lower than that indicated by HEC RAS model. Pierce (2010)yieldsapeak flood discharge with an error of 17.3%, Xu and Zhang (2009)estimates the peak discharge with an error of 30.1% andFroehlich (1995)estimates the peak discharge with an error of 39.2%. XI. Conclusions Simulation of dam breach floods is essential to characterize and identify hazards due to hypothetical dam failures. Hydraulic models such as HEC-RAS are often used for the analysis of downstream impacts resulting from potential dam failures. Estimation of the dam breach parameters, such as formation time, width and side slopes, has usually done external to the hydraulic model. 115|P a g e Dr. K Abdulrahman Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 5( Version 1), May 2014, pp.109-116 Due to uncertainty in determining the exact value of the breach formation time; different values of breach formation time have been coded into the HEC RAS 4.1 to calibrate its valueby using the observed highest water level at Chaq-Chaq Bridge.In this context a breach formation time of 1.6 hrs was achieved. Themaximum flood discharge passing through the dam breach for the corresponding breach formation time was 979.2 𝑚3 𝑠for the corresponding breach formation times. The most competitive and recent breach prediction models were examined to determine the most accurate onein predicting the breach parameters.In this context; Xu and Zhang (2009) performs better than other in predicting the average breach width and the breach formation time. Froehlich (2008) predicts the top breach width more accurate than other models andXu and Zhang(2009) is in the second order. The peak flood discharge passing the breach of the dam is under-estimated by pierce (2010) with an error of 17.3%, while Xu and Zhang (2009) over-estimates the peak discharge with an error of 30.1%. As a conclusion Xu and Zhang (2009) can be considered as the most accurate breach prediction model because it was the best in predicting the breach width and the breach formation time. [6] [7] [8] [9] [10] [11] [12] [13] XII. Acknowledgments The writer acknowledges the support from the municipality of Sulaimani especially the GIS department (Shahlaa A. F., Azad A. H., and Mohammed H.). Thanks also go to Dr. RizgarS. and Dr. NihadB.for their valuable notations. [14] www.ijera.com Chow, V. T., Maidment, D. R., Mays, L. W. (1988). Applied hydrology. McGRAWHILL. FEMA. (2005). Conduits through embankment dams. Federal Emergency Management Agency. Froehlich, D. C. (2008). Embankment dam breach parameters and their uncertainties. Journal of Hydraulic Engineering, 134(12), 1708-1721. Gee, D. M. (2010). Dam breach modeling with HEC-RAS using embankment erosion process models. World Environmental and Water Resources Congress (pp. 1347-1356). Davis: ASCE. McCuin, R. H. (2005). Hydrologic analysis and design. New Jersey: Pearson. Parhi, P. K., Sankhua, R. N., Roy, G. P. (2012). Calibration of Channel Roughness for Mahanadi River, (India) Using HECRAS Model. Journal of Water Resource and Protection, 4, 847-850. Wahl, T. L. (2004). Uncertainty of predictions of embankment dam. Journal of Hydraulic Engineering, 130(5), 389-397. Xu, Y., Zhang, L. M. . (2009). Breaching parameters for earth and rockfill dams. Journal of Geotechnical and Geoenvironmental Engineering, 135(12), 1957-1970. Yochum, S. E., Goertz L. A., Jones, P. H. (2008). Case study of the Big Bay dam failure: accuracy and comparison of breach predictions. Journal of Hydraulic Engineering, 134(9), 1285-1293. Bibliography [1] [2] [3] [4] [5] Pierce M. W.; Thornton C. I.; Abt S. R. (2010). Predicting Peak Outflow from Breached Embankment Dams. Journal of Hydrologic Engineering, 15(5), 339-349. Atallah, T. A. (2002). A review on dams and breach parameters estimation. Virginia: MSc. Thesis . Brunner, G. W. (2010 a). HEC-RAS river analysis system. User’s Manual, version4.1, Hydrologic Engineering Center. Institute for Water Resources, U. S. Army Corps of Engineers, Davis, Calif. Brunner, G. W. (2010 b). HEC-RAS river analysis system. Hydraulic Refference Manual version 4.1, Hydrologic Engineering Center. Institute for Water Resources, U.S. Army Corps of Engineers, Davis, Calif. Chow, V. T. (1959). Open channel hydraulics. New York: McGraw-Hill Book Company. www.ijera.com 116|P a g e
© Copyright 2024 ExpyDoc