Several studies have shown that European police crash reports provide different detail degrees of work zone crash-related data. In this sense, the present study aims to verify the possibility of identifying significant risk factors involved in the occurrence of road work zone crashes with casualties, based on the official data usually available, through a descriptive, binary logistic, and probit regression statistical analysis. To accomplish the analysis, a total of 2597 police-reports related to 1767 Portuguese work zone crashes that occurred during the 2013–2015 period were considered and binary logistic and probit regression models were estimated by the main type of crash, contributing factor, and driver age group. Fifteen explanatory variables, selected based on the literature review and crash data provided in police crash reports, were considered in the analysis. The results obtained for the estimated coefficients and goodness-of-fit test values were found very similar for both link functions (logit and probit) and it was possible to identify risk factors. The modeling results pointed to excessive speed, disregard for vertical signs, luminosity, intersections, and motorcycle and heavy vehicle involvement as the most significant risk factors. Given the results, it is possible to conclude that binary logistic regression can be used in the statistical analysis of the available police official work zone crash data to identify and get some insight into the risk factors involved in work zone crashes. Data analysis also revealed the need to promote adequate and complete crash report filling by police officers. While police crash reports are not revised and standardized to incorporate more detailed work zone crash information, this approach can be used to support a more efficient road operation decision making and the review of some aspects related to work zone layout design.