We study XY and dimerized XX spin-1/2 chains with random exchange couplings by analytical and numerical methods and scaling considerations. We extend previous investigations to dynamical properties, to surface quantities and operator profiles, and give a detailed analysis of the Griffiths phase. We present a phenomenological scaling theory of average quantities based on the scaling properties of rare regions, in which the distribution of the couplings follows a surviving random walk character. Using this theory we have obtained the complete set of critical decay exponents of the random XY and XX models, both in the volume and at the surface. The scaling results are confronted with numerical calculations based on a mapping to free fermions, which then lead to an exact correspondence with directed walks. The numerically calculated critical operator profiles on large finite systems (L<=512) are found to follow conformal predictions with the decay exponents of the phenomenological scaling theory. Dynamical correlations in the critical state are in average logarithmically slow and their distribution show multi-scaling character. In the Griffiths phase, which is an extended part of the off-critical region average autocorrelations have a power-law form with a non-universal decay exponent, which is analytically calculated. We note on extensions of our work to the random antiferromagnetic XXZ chain and to higher dimensions.